NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



NASA 

National Aeronautics and 
Space Administration 


Lyndon B. Johnson Space Center 

Houston. Texas 77058 



NASA CR- 


JSC-16239 

r • '979 


Space and Life Sciences Directorate 
Earth Observations Division 



AS-BUILT DESIGN SPECIFICATION OF 
SOFTWARE FOR CLUSTER-BASED PROPORTION ESTIMATION 

Job Order 76-662 


il 8 2z 10045) as-built design specification 
of software for cluster-b asfd proportion 
ESTIMATION (Lockheed Electronics Co.) 17a d 
HC A08/MF AC 1 CSCL OSp 


Prepared by 


o 3/4 


NSO-16395 


Unclas 

00045 


Lockheed Electronics Company, Inc. 
Systems and Services Division 
Houston , Texas 


Contract NAS 9-15800 


October 1979 


LEC-13934 


1. Report No. 2. Government Accession No. 

JSC-16239 

3. Recipient's Catalog No. 

4. Title end Subtitle 

As-Built Design Specification of Software for 
Cluster-Based Proportion Estimation 

5. Report Date 

October 1979 

6. Performing Organization Code 

7. Authorfi) 

Nim-yau Chu 

Lockheed Electronics Company, Inc?. 

8. Performing Organiiation Raport No. 

LEC-13934 

10. Work Unit No. 

9. Performing Organization Name and Address 
Lockheed Electronics Company, Inc. 
Systems and Services Division 
Houston, Texas 77058 

11. Contract or Grant No. 
NAS 9-15800 

13. Type of Report and Period Covered 

Design specification 

12. Sponsoring Agency Name and Addreta 1 

National Aeronautics and Space Administration 
Lyndon B. Johnson Space Center i 

Houston, Texas 77058 
J. M. Sulester, Technical Monitor 

14. Sponsoring Agency Code 

lb. Supplementary Notes 
16. Abetract 


Nine dot allocation and labeling programs, two conversion programs, and 
a utility subroutines package have been developed for use in testing, 
evaluating, and predicting the performance of Procedure 1 replacements. 
These dot allocation and labeling schemes can be tested in combination 
with cluster maps produced by such clustering algorithms as CLASSY, 
AMOEBA, and ISOCLS to determine if unsupervised clustering followed by 
such proportion estimation strategies will efficiently provide accurate 
proportion estimates. 


17. Key Words (Suggested by Authorft)) 
Bayesian dot allocation 
Dot labeling 
Majority rule 
Proportion estimation 
Proportional dot allocation 

IB. Distribution Statement 

Unclassified-Unlimited 

L 

19 Security Ossif (of this report) 
Unclassified 

20. Security Cletiif. (of this page) 
Unclassified 

21. No. of Paget 

174 

22. Price* 


JSC form 1424 IRfv 781 


For tale by the National Technical Information Service, Springfield, Virginia 22161 


NASA — JSC 











AS-BUILT DESIGN SPECIFICATION OF 
SOFTWARE FOR CLUSTER-BASED PROPORTION ESTIMATION 

Job Order 76-662 


Prepared by 
Nim-yau Chu 


Approved by 



John I. Morrow, Supervisor 
Scientific Applications Software Section 


Lockheed Electronics Company, Inc. 
Systems and Services Division 
Houston, Texas 


Contract NAS 9-15800 
for 

Eartn Observations Division 
Space and Life Sciences Directorate 
National Aeronautics and Space Administration 
Lyndon B. Johnson Space Center 
Houston, Texas 


October 1979 


LEC -13934 


CONTENTS 


Section Page 

1. PURPOSE 1-1 

2. SCOPE 2-1 

2.1 CONVERSION PROGRAMS 2-1 

2.2 DOT ALLOCATION AND LABELING PROGRAMS 2-3 

2.3 UTILITY SUBROUTINES PACKAGE 2-3 

2.4 ORGANIZATION OF DOCUMENTATION 2-5 

2.5 APPLICABLE DOCUMENTS 2-5 

3. A81: TO STRIP A GROUND TRUTH MAP 3-1 

x 

4. A86: TO STRIP A CLUSTER MAP 4-1 

5. A82: PROPORTIONAL DOT ALLOCATION 5-1 

C. A83 : PROPORTIONAL DOT ALLOCATION, 

MAJORITY RULE LABELING 6-1 

7. A84: BAYESIAN DOT ALLOCATION (UNIFORM PRIOR) 7-1 

8. A85: BAYESIAN DOT ALLOCATION (NO PRIOR) 8-1 

9. A87 : BAYESIAN DOT ALLOCATION (QUADRATIC PRIOR). ... 9-1 

10. A89: BAYESIAN DOT ALLOCATION 

(MODIFIED QUADRATIC PRIOR) 10-1 

11. A91: BAYESIAN DOT ALLOCATION (ADAPTIVE PRIOR) .... 11-1 

12. A88: BAYESIAN DOT ALLOCATION, 

MAJORITY RULE LABELING 12-1 

13. A90 : NON-BAYESIAN SEQUENTIAL DOT ALLOCATION, 

MAJORITY RULE LABELING 13-1 

14. A06: THE UTILITY PACKAGE 14-1 

15. EXAMPLE OF OPERATING PROCEDURE 15-1 


v 



Appendix 




Page 

A. 

SUBPROGRAMS 

OF 

A81 : 

TO STRIP A GROUND TRUTH MAP . . . 

A-l 

B. 

SUBPROGRAMS 

OF 

A82 : 

PROPORTIONAL DOT ALLOCATION . . . 

B-l 

C. 

SUBPROGRAMS 

OF 

A83: 

PROPORTIONAL DOT ALLOCATION, 
MAJORITY RULE LABELING 

C-l 

D. 

SUBPROGRAMS 

OF 

A84: 

BAYESIAN DOT ALLOCATION 
(UNIFORM PRIOR) 

D-l 

E. 

SUBPROGRAMS 

OF 

A85: 

BAYESIAN DOT ALLOCATION 

(NO PRIOR) 

E-l 

F. 

SUBPROGRAMS 

OF 

A87: 

BAYESIAN DOT ALLOCATION 
(QUADRATIC PRIOR) 

F-l 

G. 

SUBPROGRAMS 

OF 

A89: 

BAYESIAN DOT ALLOCATION 
(MODIFIED QUADRATIC PRIOR) . . , . 

G-l 

H. 

SUBPROGRAMS 

OF 

A91: 

BAYESIAN DOT ALLOCATION 
(ADAPTIVE PRIOR) 

H-l 

I. 

SUBPROGRAMS 

OF 

A88: 

BAYESIAN DOT ALLOCATION, 

MAJORITY RULE LABELING 

1-1 

J. 

SUBPROGRAMS 

OF 

A90: 

NON-BAYESIAN SEQUENTIAL 
DOT ALLOCATION, 

MAJORITY RULE LABELING 

J-l 


y 



Ik. 


1. PURPOSE 


The purpose of the software for cluster-based proportion estima- 
tion is to provide a standard tool for testing, evaluating, and 
predicting the performance of Procedure 1 replacements. 

In the Supporting Research Branch of the Earth Observations Divi- 
sion there are two projects underway to investigate new techniques 
to classify Landsat data for crop inventory purposes. These 
projects are based on different philosophies concerning the way 
in which an efficient inventory might be conducted. One approach, 
which supports the current corn and soybeans project, is to deter- 
mine if increasing the number of type 1 dots in Procedure 1 will 
decrease the variance of the classifier to the point that machine 
classification will provide a more accurate proportion estimate. 

The second approach, which supports a longer range effort to arrive 
at a replacement for Procedure 1, will study the performance of 
various clustering algorithms in combination with various strati- 
fied area estimations or cluster-labeling techniques. The attempt 
is to determine if unsupervised clustering followed by various pro- 
portion estimation strategies (such as cluster labeling or propor- 
tion allocation based on dots within clusters) will provide more 
efficiently proportion estimates having acceptable accuracy. The 
software described in this document is used to evaluate proposed 
replacements of Procedure 1 and thus supports the second approach.* 

In Procedure 1, dots are labeled by an analyst-interpreter. In 
the tests performed by this software system, the computer labels 
the dots using ground truth,, Thus, the variability resulting 
from the analysts is eliminated. In a normal Procedure 1 opera- 
tion, dots are chosen from a grid. In these simulations, dots 
are chosen pseudorandomly from clusters. The clusters are formed 

♦This software was developed in the Test and Evaluation Group, 

P. J. Aucoin, Jr., Group Leader. 
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using one of three algorithms. Each pixel of a cluster is equally 
likely to be selected as a dot. On the basis of cluster size and 
label (or dot label) , the proportion estimate of small grain for 
the area is then computed and compared with the ground truth 
proportion. 

The implemented system consists of a total of nine different dot 
allocation and labeling schemes. One purpose of implementing 
these schemes is to test them in combination with cluster maps 
produced by various clustering algorithms. Currently, three clus- 
ter algorithms are being tested: CLASSY, AMOEBA, and ISOCLS. The 

objective is to answer the following questions: 

a. Which dot allocation and labeling scheme is the most effective? 

b. Which of the various clustering algorithms is most effective? 

c. What are the adverse effects if fewer dots are used? 

A way to answer these questions is to study the statistical prop- 
erties of the proportion estimates. The system has a built-in 
function that allows a user to select a different pseudorandom 
sequence for dot allocation. By repeating a dot allocation and 
labeling scheme over a reasonably large number of different 
pseudorandom sequences, the following statistics may be obtained: 

a. The bias, mean squared error (MSE) , reduction in MSE, average, 
variance, and variance reduction of an estimate with respect 
to the ground truth proportion 

b. The above statistics as functions of the number of dots 
allocated 

c. The mean and variance of the number of dots allocated in 
cases where this number varies 


y* 




Another purpose of this software is to provide a basic library 
of routines to be used to construct new dot allocation and label- 
ing schemes for future testing. 

To interpret and analyze the proportion estimates is beyond the 
scope of this document. The specification, however, does estab- 
lish the baseline configuration of the software for cluster-based 
proportion estimation. 


2. SCOPE 


The system contains 11 independent programs (2 for conversions 
and 9 for dot allocation and labeling) and 1 utility subroutines 
package. These programs are written to be executed in an inter- 
active environment on the PDP 11/45 image or support processor. 
However, batch versions are available. (A description of the 
hardware used is outside the scope of this document.) Table 2-1 
lists the independent programs as they are named in their Fortran 
source files and in the executable task files. 

2.1 CONVERSION PROGRAMS 

Two conversion programs have been implemented: one for converting 

or stripping the available cluster map file, the other for the 
ground truth map file. The acceptable format of input files is 
Landsat Universal format, Fortran readable and sequentially 
accessible. The outputs are "stripped" files, Fortran readable 
and directly accessible. The stripped files are more efficient 
to access because they are only about half the size of their 
parent files and the direct access operation skips unnecessary 
disk inputs and outputs (I/O's). The name convention is as 
follows : 

a. For cluster maps, the file element is replaced by "STP" 
(stripped). For example, cluster map file, 100519101 .DT2 
is stripped to 100519101. STP. 

b. For ground truth map files, only the first four characters 
are used. The file element is replaced by "STP," For 
example, a ground truth file, 100577278 .GT0 is converted 
to 1005. STP. 




TABLE 2-1.- SOFTWARE FOR CLUSTER-BASED PROPORTION ESTIMATION 


Fortran file 

Task file 

Comment 

A10.FTN 

A81.TSK 

To "strip" down a ground truth map file 

All.FTN 

A86.TSK 

To "strip" down a cluster map file 

A01.FTN 

A82.TSK 

(A82B.TSK) 

Proportional dot allocation 

A08.FTN 

A83.TSK 

Proportional dot allocation, 
majority rule labeling 

A09.FTN 

A84.TSK 

Bayesian dot allocation (uniform prior) 

A12.FTN 

A85.TSK* 

Bayesian dot allocation (no prior) 

A13.FTN 

A87.TSK* 

Bayesian dot allocation 


(A87B.TSK) 

(quadratic prior) 

A16.FTN 

A89.TSK* 

Bayesian dot allocation 
(modified quadratic prior) 

A18.FTN 

A91.TSK* 

Bayesian dot allocation 
(adaptive prior) 

A15.FTN 

A88.TSK 

Bayesian dot allocation (uniform prior), 
majority rule labeling 

A17.FTN 

A90.TSK* 

Non-Bayesian sequential dot allocation, 
majority rule labeling 


* 

Not in original job order; added as requested. 




2 . 2 DOT ALLOCATION AND LABELING PROGRAMS 


There are nine independent programs representing nine different 
dot allocation and labeling schemes. These programs have a simi- 
lar logic flow, accept thw same input format, produce the same 
output format, and use the same utility subroutines package. A 
typical program consists of the following: 

1. Interactive input of job parameters 

2. Repetition capability to compute the statistics of a propor- 
tion estimate 

3. Intermediate summary or grand summary or both 

4. Resetting the pseudorandom number generator to any desired 
starting point 

5. Optional detailed dot files 

6. Optional status messages at the terminal 

The available dot allocation and labeling schemes are listed in 
the Comment column of table 2-1. 

2. 3 UTILITY SUBROUTINES PACKAGE 

A package consisting of utility subroutines for handling disk 
read, dot generation, and other basic operations is implemented. 
Table 2-2 lists the subroutines available. Reading the cluster 
map from its disk file is handled by subroutine CLMPCS (cluster 
map cluster sizes) and its other entry points CLMPXY (cluster 
map position — x,y) and CLMPLC (cluster map size — lines, 
columns). Reading the ground truth map from its disk file is 
handled by subroutine GTMPLB (ground truth map label). The selec- 
tion of a dot and associated bookkeeping is done by subroutine 
GETDOT. The subroutine MR (majority rule) finds a majority labei 
on a f irst-come-f irst-assigned basis. The subroutine LBLITP 
(label interpretation) interprets a given label and increments 
a counter if the label is a small grain. 


TABLE 2-2.- CONTENT OF THE UTILITY SUBROUTINES PACKAGE 


Subroutine 

Entry 

Usage 


CLMPCS 


Reads the duster map, counts the total 
number of pixels, counts the number of 
clusters, counts how many pixels each 
cluster has, stores types of clusters. 

Must be 
called 
before 
calling 



Internally, sets up a tablet number of 
pixels ae a function of line number and 
cluster number. 

CLKPXY 


CLMPXY 

Returns the line and column number of a 
given dot in any cluster using the 
internal table set up previously. 



CLMPLC 

Sets the size of the cluster map. 

default to 
117 lines 
by 19 5 
columns 

GETDOT 

RANST 

(random 

•tart) 

Chooses pseudorandomly with equal proba- 
bility a pixel from a cluster and assigns 
it as a dot. A dot array containing pre- 
viously generated dots is consulted in 
order to avoid double selection. Updates 
the dot array. 

Sets the starting point of the pseudoran- 
dom sequence. 


GTMPLB 


Returns the ground truth label of a pixel 
given its line and column position. 


MR 


Returns a majority label for a group of 
labels. 


LBLITP 


Interprets a given label as small grain or 
other. Increments a counter if the label 
is small grain. 



y* 
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2.4 ORGANIZATION OF DOCUMENTATION 


Each of the following sections (3-14) discusses one of the soft- 
ware components. The discussion begins with a description of the 
program, including its purpose and structure. This is followed 
by the program's linkage, interfaces, inputs, and outputs. A 
brief form of the algorithm is included, as is a flow chart. The 
section concludes with the listing of the program. Subprograms 
of the 10 programs having them are discussed in appendixes A-J. 

2. 5 APPLICABLE DOCUMENTS 

An applicable document is the requirements, which were informally 
transmitted under TIRF 79-0009, March 1979. 


3. A81: TO STRIP A GROUND TRUTH MAP 


3.1 DESCRIPTION 

This program converts the file format of a ground truth map to a 
form that permits faster data access. Several properties of the 
ground truth map are assumed: single channel, size 351 lines 

and 392 columns, written in Landsat Universal format, Fortran 
readable by sequential access. 

The output is a stripped ground truth map file that is Fortran 
readable by direct access. It contains no header record, and 
only 117 data records. Each record contains only 196 bytes. 

The conversion results in a reduction to one-third the number of 
lines and one-half the number of columns. The 6-to-l pixel reduc 
tion is based on a first-come-f irst-assigned majority rule. The 
lack of a header record and trailing zero-filled bytes make the 
output file much smaller, about 8-percent the size of the input 
ground truth map file. 

3.2 LINKAGE 

A subroutine named MR6 is used to compute the majority label for 
every six pixels. 

3.3 INTERFACE 
N/A 


f 


3.4 INPUTS 


The following should be entered at the terminal by the user 
(interactive) or by card image file (batch) : 

Card 

or Default Default 

line Parameter Format input value 

1 Ground truth map file name A13 None 

(However, if "STOP" 
is typed, program 
exits. ) 

The following disk file is needed: the ground truth map file. 

3.5 OUTPUTS 

This program echoes the input specifications and outputs status 
messages. 

A disk file (named with the first four characters of the input 
file name followed by .STP) is generated to contain the stripped 
ground truth map. 

3.6 BRIEF ALGORITHM: TO STRIP A GROUND TRUTH MAP 

1. Read ground truth map file name; if "STOP" is read, program 
stops. 

2. Open the ground truth map file and the output file. 

3. Repeat a through d for LINE = 1 to 117: 

a. Read three lines from ground truth map into a buffer. 
Sequential access. 

b. For K = 1 to 196, do (1) and (2) : 

(1) Extract the pixels (i, j); i = 1, 2, 3; 
j = [ (K - 1) *2 + 1] to 2 * K. 
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(2) Call MR6 to compute the majority label for the 
six pixels. 

c. Write on output file the 196 majority labels. Direct 
access. 

d. If line count reaches a certain number, output a status 
message. 

4. Write a job-ending message. 

5. Close both input and output files. 

6. Go back to 1. 


XT 

// 


3.7 LISTING 


C PROGRAM NAME* A10.FTN (A91.TSK). 

C CONCERTS GROUND TRUTH IMAGE (IN DISK FILE) OF 

C SIZE 3S1 LINES BY 392 COLUMNS TO LAC IE SEGMENT 

C SIZE 11? LINES BY 196 COLUMNS. MAJORITY RULE IS 

C USED FOR EACH GROUP OF 6 SUBPIXELS. 

C WRITTEN AND EDITTED BY NIM-YAU CHU. 

C 

BYTE Jl(3060),Li(2998),J2(S40).L2(468), 

* J3(S40) ,L3(468) , J4(S40) ,L4(468) 

BYTE IG(6).NST0P 
BYTE NAME (IS) >NGT( 10) 

EQUIUALENCE ( J1 (73) .LI ) , < J2C73) .L2) . ( J3C73) , L3) . ( J4C73) ,L4) 
EQUIUALENCE (J4(71),LINEX) 

DATA NGT(S)/’ .’/,NGT(6)/’SV,NGT(7)/’T’/,NGT(8)/’P’/ 

DATA NSTOP/’S’/ 

C 


100 CONTINUE 

WRITE (8. 201 ) 

201 FORMAT ( ’ PROGRAM *A10(A81.TSK). TO REDUCE GROUND TRUTH V 
X * RESOLUTION. OUTPUT IN STRIPPED FORMAT’/ 

X ’ TYPE GROUND TRUTH FIIE NAME’/ 

X ’ AAAAAAAAAAAAA’ ) 

READ ( 7 . 202 ) ( NAME ( K ) . K » 1 , 1 3 ) 

202 FORMAT (13A1) 

IFCNAME(l) .EQ.NSTOP) STOP 
URITE(B.203)(NAME(K).K«l.i3) 

203 FORMAT (/’ GROUND TRUTH FILE FROM ’,13A1,/ 

X ’ LABEL BY MAJORITY RULE. .. EXECUTION BEGINS. . .WAIT’ ) 

C 

OPEN (UNIT* 2. NAME* NAME. TYPE* ’OLD’ . 

X FORM- ’ UNFORMATTED ’ . READONLY . ERR -901 ) 

DO 241 K-i,4 
241 NGT ( K) -NAME ( K ) 

OPEN (UNIT-3, NAME -NGT , TYPE- ’ NEW * , FORM- ' UNFORMATTED ’ , 

X ACCESS-’ DIRECT’ ,RECURDSIZE-49,MAXREC-11?) 

C 

C 

READ(2)J1 

NLINE-117 

NCOL-196 

K3-2XNC0L+72 

KS-NLINE/6+1 

C 

DO 231 L-l.NLINE 
READ(2)(J1CK),K«1,K2) 

READ(2)(J2CK?.K-1.K3) 

READC2) (J3(K),K=i»K3) 

IX) 221 K-l.NCOL 
K1-(K-1)X2+1 
K2-KT+1 
IG(l)-LKKi) 

IG(2)»L1(K2) 

IG(3)-L2(K1) 

IG(4j-L2(K2) 

IG(5) -L3(K1 ) 

IGC6) -L3CK2) 

C WRITE(8,99101)L,K, IG 
99101 FORMAT ( ’ L.K. IG-’ , 1514) 

CALL MRS ( I G i L4 < K ) ) 

221 CONTINUE 

WRtTE(3’L) (L4(IO ,K-l«NCOL.) 

C UR I TE ( 8 , 99201 ) C J4 (K ) . K*6S, 1O0 ) 

99201 FORMAT (’ OUPUT- ' . 1614) 

I F ( MOD ( L , K5 ) . EQ . 0 ) WRITE(8.22S)L. (NGT(K) .K-1,8) 

82S FORMAT (IS,’ LINES WRITTEN ON ’ . 8A1 ) 

231 CONTINUE 


& 


ORIGINAL PAG“ is 
OF POOR QUALITY 






fiEPRODUCl/Ji. .* 

original page 


Y OP mB 

is poor 


c 

WRITEC8. 351 INLINE. ICOL. (NGTOO ,K«1 , 1©) 

351 FORMAT <’ STRIPPED MAP OF SI2E> , 14, * LINES BY*. M. ’ COLS’ . 

* ' PRODUCED IN * . 10A1 . // > 

* ’ NEXT GROUND TRUTH mf PLEASE. TO STOP. .JUST TVPE STOP’ ) 

CLOSE (UN IT-2) 

CLOSE C UNIT-3) 

GOTO 100 
C 

901 WRITECS. 911 ) 

911 FORMAT C ’ ERROR OCCURRED WHILE READING INPUT GROUND TRUTH FILE.’) 

991 STOP 
END 

C 

c 

SUBROUTINE MR6( IG. LABEL ) 

C JG f DIM 6) CONTAINS 6 G.T. SUBPIXELS 

C ON RETURN. LABEL WILL CONTAIN THE MAJORITY RUl.E LABEL. 

C 

BYTE IGC 1 ) . JG(6) . L.GC6) . LABEL 

C CLEAR BUFFER 

DO 201 1-1.6 

201 LG(I)-0 

Jl-1 

.10(1 ) * 1 Of .11 ) 

LG ( J 1 ) - 1 

C - COUNT SUBPIXELS HAUEING SAME LABELS 

DO 251 1=2.6 

DO 221 J-1.J1 

C WRITEC8. 99108) I » J. IGC I ) , JGC J) . J1 

99108 FORMAT ( ’ I, J. IG. JG. J1-' , 10M) 

IFCIGc I).EQ.JG(J)) GOTO 231 
221 CONTINUE 
J1=J1+1 
JG(J1)=IG( I ) 

LG(J1)-1 
GOTO 251 

231 LGCJ)»LG(J)+1 
£51 CONTINUE 
C 

C CHOOSE THE MAJORITY LABEL 

MAX=0 

DO 281 1=1.6 

JFiLGU ) .LE.MAX) GOTO 281 
MAX-LGC I ) 

J1 = I 

281 CONTINUE 

LABEL-JG(Jl) 

RETURN 

END 






4. A86 * TO STRIP A CLUSTER MAP 

4.1 DESCRIPTION 

This program converts the file format of a cluster map to a form 
that permits faster data access. Several properties of the cluster 
map are assumed: single channel, size 117 lines and 196 columns, 

written in Landsat Universal format, Fortran readable by sequential 
access. 

The output is a stripped cluster map file that is Fortran readable 
by direct access. It contains no header record, and only 117 data 
records. Each record contains only 196 bytes, corresponding to 
the 196 pixels of a line. 

The lack of a header record and trailing zero-filled bytes make 
the output file much smaller, about half the size of the input 
cluster map file. 

4.2 LINKAGE 
N/A 

4.3 INTERFACE 
N/A 

4.4 INPUTS 

Fortran-formatted input of the following parameter is needed: 


Card 
or 
i 


Parameters 


Format 


Default 

input 


Default 

value 


Cluster map file name 

A13 

None 




(However, 

if "STOP 



is typed, 

program 



exits. ) 




The following disk file is needed: the cluster map file 


4 . 5 OUTPUTS 

This program echoes the input specifications and outputs status 
messages. 

A disk file (named with the first four characters of the input 
file name followed by .STP) is generated to contain the stripped 
cluster map. 

4.6 BRIEF ALGORITHM: TO STRIP A CLUSTER MAP 

1. Read cluster map file name; if "STOP" is read, program stops. 

2. Open the cluster map file and the output file. 

3. Repeat a through c for LINE - 1 to 117: 

a. Read from cluster map the appropriate number of bytes 
into a buffer. Sequential access. 

b. Write on output file the just read bytes. Direct access. 

c. If line count reaches a certain number, output a status 
message. 

4. Write a job-ending message. 

5. Close both input and output files. 

6. Go back to 1. 




4.7 LISTING 


C PROGRAM All.FTN (A86.TSK) 

C TO STRIP A CLUSTER MftP IN UNIUERSAL FORMAT OF SIZE 

C 117 LINES BY 196 COLS. 

C OUTPUT IS ft FORTRAN DIRECT- ACCESSIBLE FILE WITH NO HEADER 

C RECORD NAME AS *.STP U.E. SAME FILENAME BUT MUST HAUE 

C FILE ELEMENT AS .STP ( STANDS FOR STRIPPED ) . 

BYTE JXfNSTQPs NAMEC 1?) . NGTtl?) . JJ(400) ,LL(400) 

EQUIUAL.EICE (JJ(73),LL) 

DATA NGT(10)x* . V, HiSTC 11 ) S’ ✓ . MGT< 12) T’ 0, N6TQ3)/’ P 
DATA NSTOP/’SV 
C 

100 CONTINUE 

WRITE(8. 121 ) 

121 FORMAT C* PROGRAM A1 1 ( ASC. TSK‘> . TO STRIP A CLUSTER MAP FILE.' 

* /’ INPUT CLUSTER NAP FILE MANE. AAAAAAAAAftftAA ’ ) 
READ( 7. 122>(NAME(*0,K-1, 13) 

122 FORMAT C 13A1 ) 

C 

I F ( NAME (1 ) . EQ . NSTOP ) STOP 
DO 131 K-1,9 
131 NGT(K) «NAME<K) 

C 

OPEN (UNIT-2, NAME-NAME , READONLY , ACCESS- 1 SEQUENT I AL ’ , 

* TYPE- ’ OLD* > FORM- * UNFORMATTED* ) 

NLINE-117 

NCOL-196 

NWORD- (NCC'.-1)04 +1 

OPENCUNIT-3. NAME* NGT, ACCESS-’ DIRECT’ , 

* RECORDS 1 2E-NW0RD, MAXREC-NL I NE. 

* TYPE- ’ NEW* , FORM* ’ UNFORMATTED* ) 

READ ( 2 ) J>< 

K1-NC0L+72 
K2- HL I NE^to +1 
DO 301 L.-l iNLINE 
READC2) < JJOO ,K-J K1 ) 

WRITEO’L) CLL(K) .K-l.NCOL) 

I F ( 1 10D ( L . K2 ) . EQ . 0 ) WR I TE ( 8 . 251 ) L , NGT 
251 FORMAT US,’ L. I NES WE I TTEN ON ’,17Ai) 

301 CONTINUE 
iL 

WP I TE ( 3 , 35 1 ) NL I NE . NCOL , NGT 

FORMAT C’ STRIPPED MAP OF SIZE’, 1 4.' LINES BY*,I4,* COLS', 

* * PRODUCED IN ’ , 17A1 

X ’ NEXT CLUSTER NAP PLEASE. TO STOP, JUST TYPE STOP* \ 
CLOSE i UNIT-2) 

CLOSE f UN IT- 3) 

GOTO 100 
END 


jkT 


5. A82: PROPORTIONAL DOT ALLOCATION 


5.1 DESCRIPTION 

Implemented in this piece of software is a scheme in which dot 
allocation is proportional and labeling is direct from ground 
truth. 


Dots are selected from each cluster in a pseudorandom fashion with 
equal probability. The number of dots selected from cluster i 
is given by 


where n = total number of dots to be allocated to the entire 
scene 

N^ « number of pixels in cluster i 
N = total number of pixels in the entire scene 


The equation for the proportion estimate of small grain is 

m x. N. 


sg 


= y jl jl 
ffi n i N 


where m * total number of clusters 

x^ * number of dots labeled as small grain 
(notice x^ < n^) 


The algorithm for proportional dot allocation includes the fol- 
lowing strategies: 

a. If a cluster is too small to receive at least one dot, then 
this cluster is grouped with other small clusters at the end 
of the selection process. The number of dots for this mix- 
ture; cluster is computed according to equation (1) . 


/7 


b. Since the right-hand aide of equation (1) may be a fraction 
but the number of dote must be an integer, a rounding -off 
operation ia applied to n^. 

c. Because all n^ are rounded off to the nearest integer, En^ 
may not be equal to n, the total number of dots to be allo- 
cated to the entire acene. To correct for this, dots will 
be added to or subtracted from the clusters, starting with 
the clusters that have the most dots* 

This scheme, like other allocation and labeling schemes, may be 
repeated with different pseudorandom sequences. The main program 
resets the random number generator to a different starting point 
for each repetition. 

All computation of statistics such as mean and variance is per- 
formed in the main program. An intermediate summary is printed. 

Furthermore, the program permits entry of several values of n, 
the total number of dots to be allocated, at the job initializa- 
tion stage. A grand summary is produced. 

5.2 LINKAGES 

The PDP 11/45 Fortran library (for functions and subroutines such 
as SQRT) , the PDP 11/45 system library (for routines such as the 
pseudorandom number generator) , and the utility package (described 
in section 14) are required. 

Main program 
or subprogram 

MAIN 
SAE 

(stratified area estimation) 

GETDOT 


Su broutine or function required 

SAE, GETDOT, CLMPCS, GTMPLB, LBLITP 
GETDOT, CLMPCS, GTMPLB, LBLITP 

RAN 



5.3 INTERFACE 


Interface with other routines is through the common block PRTFLG 
(print flag — 2 bytes ) , which is used to control the optional 
printing of dot files and other information. 

5.4 INPUTS 

Fortran-formatted input of the following parameters is needed: 


Card 

or 

line 


Parameters 

Format 

Default 

input 

Default 

value 

1 

Cluster map file name 

A13 

None 


2 

a. 

Number of repetitions 

13 

0 or blank 

1 


b. 

Starting point of first 
pseudorandom sequence 

15 

0 or blank 

10 


c. 

Number of repetition 
printings 

13 

0 or blank 

5 

3 

a. 

Number of " total number 
of dots" 

13 

None 



b. 

1st total number of dots 

13 

None 



c. 

2nd total number of dots 
• 

13 

None 



r . 

• 

17th total number of dots 

13 

None 


4 

Number of status messages on 
terminal 

13 

0 or blank 

NO MSG 


Two disk files are needed, one containing the stripped cluster 
map, the other the stripped ground truth map. These files should 
have been created by conversion programs A81 and A86. 



5.5 OUTPUTS 


In addition to the normal echo of Input specifications and status 
messages on the terminal, the following will be output on a line 
printer. 


Report 


Contents 


Individual repetition 1. 

2 . 

3. 


Processor header 

Ground truth small-grain proportion 

A table stowing number of clusters, 
cluster name codes, cluster sizes, and 
number of dots assigned 


4. A dot file stowing dots chosen, their 
x-y position, and their ground truth 

labels 


Intermediate summary 


1. Processor header 

2. A table stowing the estimate and bias 
of each repetition 


3. Bias, MSE, reduction in MSE, average, 
variance, variance reduction 


Grand summary 


1. Processor header 

2. A table showing bias, MSE, reduction 
in MSE, average, variance, and vari- 
ance reduction as functions of total 
number of dots 


5.6 BRIEF ALGORITHM; FOR PROPORTIONAL DOT ALLOCATION AND 
DIRECT LABELING 

1. Type messages at terminal and read from it run specifications. 

2. Compute ground truth small-grain proportion. 


y* 

9 * 


For each total number of dotr specified, repeat a and bt 

a. Rapaat (1) through (5) until all rapatitlona ara finished: 

(1) Set starting point of the pseudorandom number 
generator. 

(2) Determine each cluster's dot allocation (proportional). 

(3) Select dote from the clusters. 

(4) Pick up small-grain labels. 

(5) Compute the proportion estimate. 

b. Compute bias, MSB, reduction in MSB, average, variance, 
and variance reduction, and print an intermediate summary. 

Print a grand summary. 



. 0 / 











SK 



Figure 5-1.— Concluded 




5.7 LISTING 



PROGRAM A01 (A82.TSK) * PROPORTION ESTIMATE OF SMALL SRAIN 

USING PROPORTIONAL DOT ALLOCATION 

LABELLING IS DIRECT BY GROUND TRUTH. 


C 

C 

BYTE NAME(1S>.NGTC13) 

INTEGER NDT ( 17 ) 

REAL PSG(200) , PM( 17) , PB( 17) » RR( 1?) 

C COMMON BLOCK IS A FLAG FOR PRINTING INDIVIDUAL RUNS 

OJttWi /prtflg/jflag 


c 

c 

c READING INPUT FROM TERMINAL 

C 

WRITEC8.121) 

121 FORMAT! ’ PROGRAM i A01 CA82.TSK).’ 

* /’ PROPORTIONAL ESTIMATION OF SMALL GRAIN’ 

* /’ DOT ALLOCATION IS PROPORTIONAL. ’ 

* /’ LABELLING IS DIRECT BY GROUND TRUTH.’ 

* /’ INPUT CLUSTLR MAP FILENAME’/’ AAAAAAAAAAAAA ’ ) 
READ C? , 122 ) C NAME (K ). K- 1 , 13 ) 

122 FORMAT (13A1) 

DO 125 k-1,4 

125 N3T ( K ) “NAME ( K ) 

NGT(5) » ’ . ’ 

NGT!6)«’S’ 

NGT(7)=’T’ 

NGT(8)»’P’ 


C 

WRITE(8» 131 ) 

131 FORMAT C ’ AN ASSIGNMENT OF PSEUDORANDOM DOTS IS A REPETITION RUN 

* /’ HOW MANY REPETITION RUNS? SEED FOR FIRST PSEUDO SEQ?’ 

* /’ HOW MANY PRINTINGS OF INDIVIDUAL RUSS?’ 

X /’ III 1 1 1 1 1 III’) 

READ!?, 1 32 )JU, JSEED. JPAGE 

132 FORMAT! I3« IX. IS. IX. 13) 

IFCJU.LE.0) JV-1 

IF< JSEED. LE.0) JSEED- 10 
I F ( JPAGE . LE . 0 ) JPAGE -S 
JSKIP- < JV-1 ) / JpAGE+ 1 
C 


WRIT£(8, 141) 

141 FORMAT! ’ YOU MAY SPECIFY MORE THAN ONE TOTAL DOT NO. ( TOTDOT ) ’ 

X FOR THE RUNS’/’ HOW MANY TOTDOT S? INPUT THE TOTDOT S. 

* /. 18C1X. ’Ill’) ) 

READC7, 142)NT . (NDT(K) , K-l >NT) 

142 FORMAT C 18(13. IX) ) 

C 

WRITECS. 151 ) 

151 FORMAT C ’ STATUS MESSAGES ON TERMINAL, HOW MANY?’/’ Ill’) 

READ (7,152) NMES 

152 FORMAT (13) 

C 

WRITEtS, 1C1 ) (NAME(K) , K»1 , 13) . JV, JSEED, NT , (NDT(K) , K-l.NT) 

181 FORMAT!//’ CLUSTER FILE GIVEN - ’.13A1. 

* /’ NO. OF REPETITION RUSS FOR EACH TOTDOT -’.13, 

* /’ THE FIRST REPETITION RUN STARTS WITH SEED - ’,16, 

* /’ NO. OF TOTDOT S SPECIFIED - ’,12, 

X THEY ARE'* ,/.2X, 1714 ) 


C 

C COMPUTE PROPORTIONAL ESTIMATE FOR GROUND TRUTH 

NLINE-117 

NCOL-196 

OPEN ( UN I T -2, NAME -NGT. TYPE -’OLD’ .READONLY, FORM- ’UNFORMATTED’ , 
X ACCESS- ’ DIRECT’ ) 

DO 231 L-l.NLINE 
DO 221 K-l.NCOL 
CALL GTMPLB C L , K , LABEL ) 

CALL LBL I TP ( LABEL , LB1 , IP) 



221 CONTINUE 

231 CONTINUE 

P-FLGATC IP)/NL1NE/NC0I. 

V C 

C OPEN CLUSTER FILE 

OPEN C UN I T-l, NAME-NAME. TYPE® ‘OLD’ , 

* READONLY , ACCESS 3 ’ D I RECT ' . FORM® ’ UNFORMATTED * ) 

CALL CLMPLCCNLINE.NCOL) 

C. 

C START ESTIMATION FOR EACH TOTDQT AND REPETITIONS 

JNES-0 

DC 361 K-l.NT 
JS-JSEED 
DO 331 J-l.JU 
JFLAG®0 

IF(HOD(J-1,JSKIP).EQ.0) JPLAG-1 
IFCJFLAG.ES. 1) UPITEC6 32DNAME.P 

32 1 FORMAT C ’ 1 ’ . 1 OX . ’ ACCURACY ACCES5EMENT SOFTWARE < 4-25-79 ) ’ , / , 2X , 

* ’PROGRAM A01 (A82.TSK) ‘ ESTIMATION OF SMALL GRAINS PROPORTION’ 

* />15X. ’DOT ALLOCATION sPROFORT I ONAL, ’ 

* /, 15X. 'DOT LABELL I NG > D I RECT BY GROUND TRUTH.’ 

* //,10X, ' II'fPUT CLUSTER MflP IS FILE M5A1. 

X /.10X,’ TROFORT I ONAL ESTIMATE FOR GROUND TRUTH IS’,F8.5) 

IFCJFLAG.ES. 1) HR I TE ( 6 > 322 ) J . JSEED 

322 FORMATC/, 2X. ’ REPETITION RUN*’ , 13. 

* ’ RANDOM DOT SEQUENCE STARTS WITH’ , 18. ’ ’ ) 

CALL RANSTC JSEED) 

CALL SAECNDTOO.PSGCJ)) 

JMES=JMES+1 

I F t JMES . LE . NMEf ) WRITEC6. 333 ) NDT ( K> , J > JSEED , PSG < J ) 

323 FORMAT C * TOTDOT® ’ . 13. ' REPETITION RUN-’. 13, ’ SEED®’, 16, 

* ' ESTIMATE®’ ,F8.5) 

IFCJFLAG.EQ. 1) WRITEC6, 324)NDT(K ) . PSGC J) 

324 FORMAT C / ’ FOR TOTAL. DOT NO.®’ > 14; ’ PROPORTION ESTIMATE® ’ , F8.5) 
331 JSEED=JSEED+150 

C 

WRITE (6, 321 )NAME. P i 

WR I TE C 6 > 334 ) NDT C K ) ) 

334 FORMATC/’ SUMMARY OF REPETITION RUNS FOR TOTDOT NO. ®’,I4. 

* ’ ’ . • . 10X» 

* ’REPETITION PSEDUO SEQ SMALL GRAIN BIAS’, /,10X, 

X ’ RUN SEED ESTIMATE WRTG.T.’) 

r COMPUTE BIAS AND M.S.F.. 

PBCK)®0. 

PMCfO-0. 

DO 341 J-l.JU 
TEMP = PSG ( J ) -P 
PBCK) ®PBCK)+TEMP 
PM C IO = PM C K ) +TEMP**2 
WRITECS,33b)J, JS, PSGC J), TEMP 
336 FORMAT (13X, 13, 6X, I6,5X,F8.5,3X,F9.6) 

341 JS-JS+1S0 

PBao=PB(iO/'Ju 
PM(K)=FM\K)/JU 
AUERG = PB ( K ) -i-P 

I F C JU . EQ . 1 ) UAR® PM C K ) -PB C K ) *.*2 
IFCJU.GT.l) UAR® CPMCkO-PBCK)**2)*JU/(JU-l) 

C COMPUTE UARIANCE REDUCTION 

RR C K ) = PM C K i / C P* < 1 . -F ) /NDT < K ) ) 

REDUAR=UAR/C P* C 1 . -F ) /ND’I CIO ) 

WR I TE C 6 , 357 ) PB C IO , PM C K ) , RR C K ) . AUERG , UAP. , REDUAR 
3S7 FORMATC/’ BIAS® ’ ,F10. 6, ’ M. S. E. = ’ , F10. £, ’ REDUCTION® ’. F10. £• 

* /’ AUERAGE®’ ,F10.6, ’ UARIANCE®’ ,F10. 6. ’ REDUCTION®’ .F10. 6) 
361 CONTINUE 

C 

C— PRINT GRAND SUMMARY FOP THIS JOB 

WR1 TECS, 321 )NAME.P 
WRITE) 6, 371) NT, JU 


OR!GI NAl 
0F POOR 



> 






371 FORMAT (.', 14X» ' GRAND NUMMARY QF THIS JOB - 

* /,1SX.' NO, OF DIFFERENT TOTDOT S «',13, 

* A14X, ' NO. OF REPETITION RUNS PER TOTDOT -‘.13. 

* a/, IX. ’TOTDOT BIAS M.S.E. REDUCTION’ « 

* * AUERASE VARIANCE MAR REDUCTION' ) 

DO 331 K-l.NT 

AUERG-PBOO+P 

IF(JU.EQ.l) UAR- PM<K)-PB(K)**2 

IF< JU.GT. 1) UAR- <PM(K>-PB(K)X*2)*JU/(JU-1> 

REDUAR-UAR/( P*C l.-P)XNDT(K) ) 

331 t«ITE<6,382>NDT<K),PB<K),PM(K),RR<K),AUERG.UAR.REDUAR 

332 FORMAT* 2X, 13, 2X.F3.5, 4(F10.6). 3X.F10.6) 

C 

C 

WRIT£<6,401> 

401 FORMAT ( ’ 1 END OF THIS JOB * ) 

C 

STOP 

END 

r 

c 

SUBROUTINE SAE(NDOT.PSG) 


C STRATIFIED AREAL ESTIMATION USING PROPORTIONAL DOT ALLOCATION. 

C GIUEN ANY TOTOAL NO. OF DOTS IN NDGT, PROGRAM WILL RETURN 

C THE PROPORTIONAL ESTIMATE IN PSG. NOTICE THE FOLLOWING FEATURES « 

C (1) ALTHOUGH DOTS ARE ALLOCATED PROPORTIONAL TO CLUSTER SIZES, 

C ACTUAL NO. OF DOTS ARE ROUND OFF. THOSE CLUSTERS GET LESS 

C THAN ONE COT ARE GROUPED TOGETHER AND COTS ARE REALLOCATED 

C TO THIS EXTRA CLUSTER. AFTER THIS, DEDUCT FROM OR ADD DOTS 

C TO LARGEST CLUSTERS TO MARE FINAL NO. OF DOTS ALLOCATED 

C EQUAL TO NO. OF HOTS DESIRED (NDGT) . 

C (2) THE RANDOM NO. GENERATOR CAN BE SET TO ANY STARTING POINT 

C BY CALL FANSTC I START ) . DOT GENERATION WILL BE USING 

C THE I START TH RANDOM NUMBER. 

C IF NO CALL TO RANST IS WIDE, RANDOM GENERATOR ALWAYS STARTS 

C FROM Jl-0 AND J2-0. 

C 

C WRITEN AND EDITED BY N.Y. CHU ON 4-2-79. 


C 

C 

INTEGER MM(51),NN(51).NDARY<200) 

INTEGER LL(Sl), JJC10) 

BYTE LABEL 

C COMMON BLOCK IS A FLAG FOR PRINTING INDIUIDUAL PUNS 

COMMON /PRTFLG/JFLAG 
DATA HSTART'0/ 

I F ( NDOT . GT . 200 ) GOTO 911 

C 

C GET PIXEL COUNT FOR EACH CLUSTER 

I F ( NSTART . NE . 0 ) GOTO 121 
NSTART-1 

CALL CLMPCSCNPIXEL.M.m.LL) 

IFCM.GE.51) GOTO 901 
121 CONTINUE 


C DETERMINE NO. OF COTS FOR EACH CUSTER 

C 

c BEGIN PROPORTIONAL DOT ALLOCATION 

C 


Kl-0 

DO 211 K=i-n 

NN t K ) -FLOAT C NDOT ) #MM ( K ) .-'HP I XEL+0 . 5 
K1-K1+NNCK) 

C WRITE(b> 99101 )K> NPIXEL, MN(K) > NNiK) 
99101 FORMAT ( ’ K, NPIXEL-' .214, ’ MN,NN-’,2I4) 
;?tl CONTINUE 



C uKPu 1 h jUC'ICILITY OF THE 

C- CHECK AND GROUP THOSE CLUSTERS GETTING NO DOT ORIGINAL PAGE IS POOR 

K2-0 

K3-0 

DO 221 K-l.M 
1F<HNCK).GT.0) GOTO 221 
K2-K2+MMCK) 

K3-K3+1 

221 CONTINUE 
K3-M-K3 
MX*M 

I" CREATE AN EXTRA CLUSTER AND AS THE 01+1) TH CLUSTER 

IFCK2.ECJ.0) GOTO 231 

K =NDOT*FLOAT < K2 ) /NP I XEL+0 . S 

IFCK.EQ.0) GOTO 231 

MX-MX+1 

MMCMX)«K2 

NN(MX)»K 

Kl-Kl+K 

231 CONTINUE 
C 

C CHECK IF NO. OF ft: SIGNED DOTS EQUAL TO THE SPECIFIED 

I F (ND0T-K1 ) 241 . 281 . 24S 
241 INC— 1 

GOTO 249 
245 INC-1 

249 CONTINUE 
C 

C ADD TO OR DEDUCT FROM LARGE CLUSTERS TO 

C MAKE NO. OF DOTS ASSIGNED EQUAL TO THE SPECIFIED 

K1.-IABSCND0T-K1) 

L-0 

250 IFCK1.EQ.0) GOTO 281 
L-L+l 

Kl-Kl-1 
I TEMP-0 
DO £61 K-l.M 
IFfL.-EQ.l) GOTO 252 
K2-L-1 

DO 251 J-1.K2 
IFCK.EQ.JJ(J)) GOTO 281 

251 CONTINUE 

252 IF<NN(K).LE. ITEMP) GOTO 261 
ITEMP-NN(K) 

JJCL) *K 

261 CONTINUE 

NNlJJ(L) )»NNCJJ(L))+INC 
C UR I TE l 6 , 99108 ) L , NDOT . J J CL ) . NN < JJ ( L ) ) 

99106 FORMAT ( ’ 99106 L.NDOT. JJtL) .NN-’ ,514) 

IFtL.EQ.K3) L-0 
GOTO 250 

281 CONTINUE 
C 

C WRITE CLUSTER INFORMATION 

IFCJFLAG.EQ. 1) WRITEC6. 282) NDOT, NPIXEL, MX 

282 FORMAT (//.EX, ’ TOTAL NO. OF DOTS AUAILABLE -’,16, 

X x.SX.’ TOTAL NO. OF PIXEL IN THE CLUSTER MAP -’,16, 

* /',5X, ’ TOTAL NO. OF CLUSTERS -M6, 

* /'.10X. ’ CLUSTER CLUSTER NO. OF NO. OF DOTS’ 

* 10X, ’ NO. CODE PIXELS ASSIGNED’) 

DO 203 K-l.MX 

IFCJFLAG.EQ. 1) WRITEC6. 284 )K, LLCK) ,MM(K) ,NN<K) 

283 CONTINUE 

204 FORMAT! 13X. 12, 5X. 14, 3X, 16, 4X, 13) 

1FCMX.EQ.M) GOTO 288 
KX-0 

DO 285 K-l.M 

IFCNNCK) .NE.0) GOTO 285 

KX-KX+1 

NDARY ( KX ) -K 


47 


205 CONTINUE 

IFCJFLAG.EQ. 1) WRITE(6»286>MX,<NDARYa0.K-l,K>O 

286 FORMAT</*L THOSE SMALLER CLUSTERS GROUPED (AS CLUSTER*. 13. 

* ’ ) ARE i * . 3(/, 16(13. ',*)) ) 

268 CONTINUE 

C 

C END OF PROPORTIONAL DOT ALLOC I AT I ON 

C 

C WR1T£(6. 99005) (NN(K), K-l.MX) 

99605 FGRMATC NN-’.ISM) 

C 

C FOR EACH CLUSTER. ASSIGN DOTS 

ND-1 

DO 291 K-l.MX 
Kl-NN(K) 

IF(Kl.EQ.O) GOTO 291 
DO 289 L-1.K1 

289 CALL GETDOT(MMCK) .L.NDARY(ND) . K3) 

99191 KK1-ND+K1-1 

C WRITEC6, 99192) (NDARY(L) .L-ND.KK1) 

99192 FORMAT! * NDARY-* ,1514) 

291 ND-ND+K1 
C 

C FOR EACH CLUSTER. FIND LINE AND COL NO. OF EACH ASSIGNED DOTS 

IFtJFLAG.EQ. 1) WRITEC6.29S) 

295 FORMAT (//’ DOT FILE * 

* CLUSTER DOT NUMBER POSITION GROUND TRUTH LABEL 

* /’ NO. UIRT ITS GLISTER LINE. COL RAW, CODE*) 
PSG-0. 

ND-1 

DO 331 K-l.MX 
Kl-NNCK) 

IFCK1.EQ.0) GOTO 331 
ISG-0 

DO 321 L-1.K1 
K3-NDARVCND) 

c to TEST FOR THE EXTRA CLUSTER 

IFCK.LE.M) GOTO 311 

C ALG. TO GET CORRECT LABEL FOR THE EXTRA CLUSTER 

DO 305 J-l.M 

IFCNNC J) .NE.0) GOTO 305 

IF(K3.LE.MMCJ>) GOTO 307 

K3-K3-MMCJ) 

rnNT t ni if 

307 CALL CLMPXYC J.K3.NDL.NDC.LL) 

CALL GTMPLBCNDL.NDC, LABEL) 

CALL LBLITPCLABEL.LB1, ISG) 

IFCJFLAG.EQ. 1) WRITEC6. 300) J. K3. NDL, NDC. LABEL. LB1 

308 FORMAT (3X. 12. 8X. 16. 7X, 13, 14. 5X, 14, IS) 

GOTO 321 

311 CALL CLMPXY ( K , K3 > NDL • NDC . LL ) 

CALL GTMPLB ( NDL . NDC . LABEL ) 

CALL LBL I TP C LABEL . LB 1 , ISG) 

IFCJFLAG.EQ. 1) WRITEC6, 308)K, K3. NDL. NDC. LABEL, LB1 
321 ND-ND+1 

c COMPUTE S.G. ESTIMATE 

PSG-PSG+ FLOAT ( ISG )/Ki*MMCK)YNPIXEL 
331 CONTINUE 
GOTO 990 
C 

901 WRITECS. 902)M 

902 FORMAT ( ’ YOU HAUE’ » 13. ’ TOO MANY CLUSTERS (MAX -50) * ) 

GOTO 990 

911 WRITECS. 912) NDOT 

914 FORMAT ( ’ YOU SPECIFIED’ , 12, ’ TOO MANY DOTS (MAX' -200) ’ ) 

C 

990 RETURN 
END 


^ i if 


6. A83 ; PROPORTIONAL DOT ALLOCATION, MAJORITY RULE LABELING 


6,1 DESCRIPTION 

Implemented in this piece of software is a scheme in which dot 
allocation is proportional and labeling is by ground truth using 
majority rule. 

Dots are selected from each cluster in a pseudorandom fashion 
with equal probability. The number of dots selected from a 
cluster i is given by 

N i 

n. = n (6-1) 

where n = total number of dots to be allocated to the entire 
scene 

N^ = number of pixels in cluster i 

N * total number of pixels in the entire scene 

For each cluster, all the dot labels are compared, then the 
cluster is labeled using majority rule. 


The equation for the prooortion estimate of small grain is 



cluster i 
labeled as 
small grain 


(6-2) 


The algorithm for proportional allocation includes the following 
strategies: 

a. If a cluster is too small to receive at least one dot, then 
this cluster is grouped with other small clusters at the end 
of the selection process. The number of dots for this mix- 
ture cluster is computed according to equation (1) . 




b. Since the right-hand aide of equation (1) may be a fraction 
but the number of dots must be an integer, a rounding-off 
operation is applied to n^. 

c. Because all n^ are rounded off to the nearest integer, En^ 
may not be equal to n, the total number of dots to be allo- 
cated to the entire scene. To correct for this, dots will 
be added to or subtracted from the clusters, starting with 
the clusters that have the most dots. 

This scheme, like other allocation and labeling schemes, may be 
repeated with different pseudorandom sequences. The main program 
resets the random number generator to a different starting point 
for each repetition. 

All computation of statistics such as mean and variance is per- 
formed in the main program. An intermediate summary is printed. 


Furthermore, the program permits entry of several values of n, 
the total number of dots to be allocated, at the job initializa- 
tion stage. A grand summary is produced. 


6.2 LINKAGES 

The PDP 11/45 Fortran library (for functions and subroutines such 
as SQRT) , the PDP 11/45 system library (for routines such as the 
pseudorandom number generator) , and the utility package (described 
in section 14) are required. 


Main program 

or subprogram Subroutine or function required 


MAIN 

SAEMR 

GETDOT 


SAEMR, GETDOT, CLMPCS, GTMPLB, LBLITP 
GETDOT, CLMPCS, GTMPLB, LBLITP, MR 
RAN 



6.3 INTERFACE 


Interface with other routines is through the common block PRTFLG 
(2 bytes), which is used to control the optional printing of dot 
files and other information. 

6.4 INPUTS 

Fortran- formatted input of the following parameters is needed: 
Card 

or Default Default 


line 


Parameters 

Format 

input 

value 

1 

Cluster map file name 

A13 

None 


2 

a. 

Number of repetitions 

13 

0 or blank 

1 


b. 

Starting point of first 
pseudorandom sequence 

15 

0 or blank 

10 


c. 

Number of repetition 
printings 

13 

0 or blank 

5 

3 

a. 

Number of "total number 
of dots" 

13 

None 



b. 

1st total number of dots 

13 

None 



c. 

2nd total number of dots 
• 

13 

None 



r. 

• 

• 

17th total number of dots 

13 

None 


4 

Number of status messages on 

13 

0 or blank 

NO MSG 


terminal 

Two disk files are needed, one containing the stripped cluster 
map, the other the stripped ground truth map. These files should 
have been created by conversion programs A81 and A86. 



6.5 OUTPUTS 


In addition to the normal echo of input specifications and status 
messages on the terminal, the following will be output on a line 
printer. 


Report 


Contents 

Individual repetition 

1. 

Processor header 


2. 

Ground truth small-grain proportion 


3. 

A table showing number of clusters, 
cluster name codes, cluster sizes, 
and number of dots assigned 


4. 

A dot file showing dots chosen, their 
x-y position, and their ground truth 
labels 


5. 

A table showing the majority labels 
for the clusters 

Intermediate summary 

1. 

Processor header 


2. 

A table showing the estimate and bias 
of each repetition 


3. 

Bias, MSE, reduction in MSE, average, 
variance, variance reduction 

Grand summary 

1. 

Processor header 


2. 

A table showing bias, MSE, reduction 
in MSE, average, variance, and vari- 
ance reduction as functions of total 
number of dots 



6.6 BRIEF ALGORITHMS FOR PROPORTIONAL DOT ALLOCATION AND 
MAJORITY RULE LABELING 


1. Type messages at terminal and read from it run specifications. 

2. Compute ground truth small-grain proportion. 

3. For each total number of dots specified, repeat a and bt 

a. Repeat (1) through (6) until all repetitions are finished: 

(1) Set starting point of the pseudorandom number 
generator. 

(2) Determine each cluster's dot allocation (proportional). 

(3) Select dots from the clusters. 

(4) Pick up small-grain labels. 

(5) Find the majority label for each cluster. 

(6) Compute the proportion estimate according to (5). 

b. Compute bias, MSB, reduction in MSE, average, variance, 
and variance reduction, and print an intermediate summary. 

4. Print a grand summary. 
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Figure 6-1. — Concluded. 



6.7 LISTIHG 


C 

C PROGRAM A0Q?A63.TSK ) » _ PROP ORTION ESTIMATE OF SMALL GRAIN 

C — DOT ALLOCATION IS PROPORTIONAL 

t— DOT LABELLING IS MAJORITY RULE BY GROUND TRUTH 

C 

c 

c 

BYTE NAME(1S).NGT(13) 

INTEGER NDTU7) 

REAL PSG(20O)*Ft1<17)»PB(l?)»RR(l7' 

C — COMMON BLOCK IS A FLAG FOR PRINT ENABLED I SABLE 

COMMON ’'PRTFLG/'JFLAG 

DATA NOT<S>/\MN6T?6)/*S’/,N6T<7)/*TV.NGT(8)/*PV 

L 

C READING INPUT FROM TERMINAL 

C 

URIT£?8, 121 • 

121 FORMAT? • PROGRAM* P«98(rte3.TSK> . * 

* x* PROPORTIONAL ESTIMATION OF SMALL GRAIN* 

* x* DOT ALLOCATION IS PROPORTIONAL. ’ 

* /' LABELLING IS MAJORITY RULE FROM GROUND TRUTH. * 

* /» INPUT CLUSTER MAP FILENAME’/* AAAAAAAAAAAAA* ) 

READ < 7 , 1 22 H NAME ( K ) , K ■ 1 . 1 3 ) 

122 FORMAT? 13Ai ) 

DO 125 K-1.4 

12S NGT(K)»NAME(K) 

C 

UR1TE?8. 131 ) 

131 FORMAT? * AN ASSIGNMENT OF PSEUDORANDOM OUTS IS A REPETITION PUN 

* s' HOH MANY PEPFT I T I Ot RUNS? SEED FOR FIRST PSEUDO SEQ?* 

* s' HOW MANY PRINTINGS OF IND1UIDUAL RUNS' 7 ’ 

*• s' in mu in*) 

READ ( 7 . 1 3? ) JU , JSEED . JPAGE 

132 FORMAT? 13. IX. 15. IX. 13) 

IF(JU.LE.C) JU-1 
IF(JSEED.LE.G) JSEED-10 
IF? JPAGE.LE.0) J PAGE'S 
JSK I P- ? J?> 1 ) ' JFA6E+ 1 

C 

WRITECS. 141) 

141 FORMAT? * VOU MAY SPECIFY MORE THAN ONE TOTAL DOT NO. ? TQTDOT ) * 

* . * FOR THE RUNS’. '* HOW MANY TOTDOT S' 7 INPUT THE TOTDOT S. 

* /. IS? IX.* 1 1 1*) ) 

READ?7.142)NT. INDT0O.K-1.NT) 

142 FOPMAT? IB? 13. IX) ) 

C 

WRITE?S.l£i) 

151 FORMAT?* STATUS MESSAGES ON TERMINAL. HOW MANY Vs' HI’ ) 

READ? V. 1S2)NMES 

152 F0W1AT? 1 3) 

C 

WRITE? 8. 181 ) (NAME(K) • K- 1 • 13) • JU. JSEED.NT. <NDT?K) .K-l .iVi ) 

181 FORMAT?//' CLUSTER FILE GIUFN • * . 13A1, 

* /* NO. OF REPETITION PUNS FOR EACH TOTDOT -M3. 

* /’ THE FIRST REPETITION RUN STARTS WITH SEED - * . 16. 

* x* NO. OF TOTDOT S SPECIFIED ■ M2. 

* *. THEY ARE* ' ./.2X, 1714 ) 

C 

C COMPUTE PROPORTIONAL ESTIMATE FOR GROUND TRUTH 

NLINE-11? 

NCOL-196 

0PEN?UN:T»2.Nf*1E-MGT.TYPE-*0LD* . READONLY. FOPM-* UNFORMATTED* . 
t ACCESS- 'DIRELT' ) 

DO 231 L-l.NLINE 
DO 221 K-l.NDJL 
CALL GTMPLB'L.K. LABEL) 

CALL LBLITPCLABEL.LB1. IP) 


& 


reproducibility op 
original Pag I*' t ° ^ 


the 


221 CONTINUE 

231 CONTINUE 

P - FLOAT UP) UHL I NE-'NCQL 

c 

OPEN C UN I T= 1 . WIE-NfttlE > TYPE* ’ OLD’ , FORM* ’ UNFORMATTED’ . 

* READON_Y,ACCESS«’ DIRECT’ ) 

CALL CLMPLCCNLINE.NCOL) 

C 

r, START ESTIMATION FOR EACH TOTDGT PNC REPETITIONS 

JMES-0 

DO 361 K-l.NT 
JS-JSEED 
DO 331 J-l.JU 
JFLAfi-0 

IFCL .3 ( J , JSK I P ) . EQ . 1 ) JFLAG-1 
I Ft .Fi. AG. in. 1) URITEC6.321)NAME,P 

321 FOR i . C ’ 1 ’ . 10X. ’ ACCURACY ACCESSEMENT SOFTWARE C 4-28-79 ) ’ 

* . , ’PROGRAM A08: ESTIMATION OF SMALL GRAINS PROPORTION’ 

Y '■ , 1SX , ’ DOT ALLOCAT I ON : PROPORT I CM. , ’ 

* / % 15X, ’ DOT LABELLING: MAJORITY RULE BY GROUND TROTH.’ 

* v,l 0X,’ INPUT CLUSTER MAP IS FILE ’ ,1SAI, 

* /', 10X, ’ PROPORTIONAL ESTIMATE FOR GROM) TRUTH IS* ,F8.5) 
IFCJFLAG.EQ.l) URITEC6, 322) J, JSEED 

322 FORMAT C U 2X > ’ REPETITION RUN- ’, 13, 

* ’ RANDOM DOT SEQUENCE STARTS WITH’, 18, ’ ’) 

CALL RANSTCJSEED) 

CALL SAEMRCNDTCK) , FSG( J) ) 

JMES-JMES+1 

IFCJMES.LE.NMES) URITEC8, 323)NDTCK) . J, JSEED.PSGC J) 

323 FORMAT C ’ TGTDOT- ’ ■ 13, ’ REPETITION RUN- ’, 13, ’ SEED-’, IS, 

* ’ ESTIMATE-’ ,F8. 5) 

IFCJFLAG.EQ.l) MRITEC6, 324)NDT(K) > PSGC J) 

324 FORMAT ( / ’ FDR TOTAL DOT NO.-’, 14,’ PROPORTION ESTIMATE- ’, F8. 5) 

331 JSEED- JSEED+ 150 

C 

UR I TE C G , 321 ) NAME > P 
UR I TE C fo , 334 ) NDT C K ) 

334 FORMAT C/’ SUMMARY OF REPETITION RUNS FOR TOTDOT NO.- ’,14, 

Y ’ REPET 1 T I ON ‘ ’ PSEDUO SEQ SMALL GRAIN BIAS’ , /, 10X , 

Y ’ RUN SEED ESTIMATE WRTG.T. ’) 

C COMPUTE BIAS AND M.S.E. 

PB(K)=0. 

FT1CK) =0. 

DO 341 J-l.JU 
TEMP = PSG C J ) -P 
PB CIO -PB C K ) +TEMP 
PMCK)=PMCK)+TEMP*Y2 
UR I TE C 6 , 336 ) J , JS , PSG C J ) , TEMP 
336 FORMAT f. 13X, 13, SX, IS, SX. F8. S, 3X, F9. 6) 

341 JS-JS+150 

RBCK)=PBCK)/JU 
PM ( K ) = PM C K ) / JU 
AUERG-PB C K ) +P 

IFCJU.EQ.l) UAR- FMCIO-PBCK)YY2 
IFCJU.GT.l) UAR- C PM C K ) -PB C K ) Y Y2 ) * JU/ ( Ji J- 1 ) 

C COMPUTE UARIANCE REDUCTION 

RRCK)=PMCK)/C PY: C 1 . -P ) /NOT C K ) ) 

REDUAR-UAR/ C PY C 1 . -P ) /'NDT C K ) ) 

UR I TE C 6 , 35? ) PB C K ) , PM C K ) , RR C K ) . AUERG , UAR , REDUAR 
357 FORMAT C ’ BIAS- ’ ,F10.6, ’ M.S.E. -’ ,F10. 6, ’ REDUCTION- ’, F10. 6, 

Y /’ AUERAGE-’ ,F10.6, ’ UARIANCE- ’, F10. 6, ’ REDUCTION- ’, F10. 6) 
361 CONTINUE 

C 


- PRINT GRAND SUMMARY FOR THIS JOB 
UR I TE C S > 32 1 ) NAME , P 
UR ITE C 6 , 371 ) NT , JU 


37 


371 FORMAT ( / . 1 4X , ’ GRAND SUMMARY OF THIS JOB 

* /.1SX, ’ NO. OF DIFFERENT TOTDOT S -M3. 

* /. 14X, ’ NO. OF REPETITION RUNS PER TOTDOT -M3. 

* //.IX. * TOTDOT BIAS M.S.E. REDUCTION’. 

* ’ AUERAGE UARIANCE UAR REDUCTION’ ) 

DO 381 K-l.NT 

AUERG«PB(K)+P 

IF(JU.EQ.l) UAR- PM(K)-PB(K)*X2 
IF(JU.GT.l) UAR- (PM(K)-PB(K)*#2)*JU/(JU-1) 
REDUAR-UAR/C P*C1.-P)/NDT(K) i 

381 WRITE(6,382)NDT(K) , PB(K) «PM(K) »RR(K) .AUERG.UAR.REDUAR 

382 F0RMAT(2X. 13. 2X.F9.5. 4(F10.6), 3X.F10.6) 

C 

c 

WRITECS. 401 ) 

401 FORMATS 1 END OF THIS JOB ’) 

C 

STOP 

END 

C 

C 

SUBROUTINE SAEMR C NDOT . PSG ) 


i; STRATIFIED AREAL ESTIMATION USING PROPORTIONAL DOT ALLOCATION, 

C AND MAJORITY LABELLING BY GROUND TRUTH. 

C GIUEN ANY TOTOAL NO. OF DOTS IN NDOT, PROGRAM WILL RETURN 

C THE PROPORTIONAL ESTIMATE IN PSG. NOTICE THE FOLLOWING FEATURES < 

C (1) ALTHOUGH COTS ARE ALLOCATED PROPORTIONAL TO CLUSTER SIZES. 

C ACTUAL NO. OF DOTS ARE ROUND OFF. THOSE CLUSTERS GET LESS 

C THAN ONE DOT ARE GROUPED TOGETHER AND DOTS ARE REALLOCATED 

C TO THIS EXTRA CLUSTER. AFTER THIS, DEDUCT FROM OR ADD DOTS 

C TO LARGEST CLUSTERS TO MAKE FINAL NO. OF DOTS ALLOCATED 

C EQUAL TO NO. OF NOTS DESIRED (NDOT). 

C (2) THE RANDOM Id. GENERATOR CAN BE SET TO ANY STARTING POINT 

C BY CPU. RANSTC I START) . DOT GENERATION WILL BE USING 

C THE I START TH RANDOM NUMBER. 

C IF NO CALL TO RANST IS MADE, RANDOM GENERATOR ALWAYS STARTS 

C FROM Jl-0 AND J2-0. 

C 

C WRITFN AND EDITED BY N.Y. CHU ON 4-26-79. 

C 


C 

I INTEGER MM(S1).NN(51), NDARY ( 200 ) , LL ( 5 1 ) 

INTEGER JJ(30).LG(60) 

BYTE LABEL . LBL ( 200 ) > JR<60) 

C COMMON BLOCK IS A FLAG FOR PRINT ENABLE/DISABLE 

COMMON /PRTFLG/JFLAG 
DATA NSTART/0/ 

C 

I F ( NDOT . GT . 200 ) GOTO 911 
C 

C GET PIXEL COUNT FOR EACH CLUSTER 

I F ( NSTART . NE . 0 ) GOTO 121 
NSTART-1 

CALL CLMPCS ( NP I XEL . M , MM , LL ) 

IFCM.GE.51) GOTO 901 
121 CONTINUE 
C 

C DETERMINE NO. OF DOTS FOR EACH CLUSTER 

C 

C BEGIN ROPORT I ONAL DOT ALLOCATION 

C 

Kl-0 

DC 211 K-l.M 

NN < K ) « FLOAT ( NDOT ) *MM ( K ) /NP I XEL+0 . 5 
Kl-Kl+NNCK) 

211 CONTINUE 
C 

C CHECK AND GROUP THOSE CLUSTERS GETTING NO DOT 

K2-0 

K3-0 



DO 221 K-ltN 
IF(NN(K).GT.0) GOTO 221 
K2-K2+MHCK) 

K3-K3+1 

221 CONTINUE 
K3-M-K3 
MX-M 

C — CREATE AN EXTRA CLUSTER AND AS THE (M+l) TH CLUSTER 

IF(K2.EG!.0) GOTO 231 
K -NDOT*FLQAT (K2)/NPI XEL+0 . 5 

IFCK. EQ.0) GOTO 231 
MX-MX+1 
MM(MX)-K2 

ISN(MX) *K 
Kl-Ki+K 

231 CONTINUE 
C 

C CHECK IF NO. OF ASSIGNED DOTS EQUAL TO THE SPECIFIED 

IF(NDOT-Kl ) 241. 281 >245 
241 INC— 1 

GOTO 249 
24S INC-1 

249 CONTINUE 
C 

C ADD TO OR DELICT FROM LARGE CLUSTERS TO 

C MAKE NO. OF DOTS ASSIGNED EQUAL TO THE SPECIFIED 

Kl-IABS(NDOT-Ki) 

L-0 

250 IF(K1.EQ.0) GOTO 281 
L-L+l 

Kl-Hl-1 
I TEMP-0 
DO 261 K-l.M 
IF(L.EQ.l) GOTO 252 
K2-L-1 

DO 251 J-1.K2 
I F ( K . EQ . JJ ( J ) ) GOTO 261 

251 CONTINUE 

252 IF(NNCK) .LE. ITEMP) GOTO 261 
ITEMP-NN(K) 

JJ(L)=K 

261 COST I NUE 

NNC JJ(L ) ) -NN ( J J ( L ) ) + I NC 

IFCL. EQ.K3) L-0 
GOTO 250 

281 CONTINUE 
C 

C WRITE CLUSTER INFORMATION 

IF(JFLAG.EQ.l) WRITEC6, 282)ND0T, NPIXEL. MX 

282 FORMAT < X/-, 5X. ’ TOTAL NO. OF DOTS AUAILABLE -’.16, 

X / , 5X , ' TOTAL NO. OF PIXEL IN THE CLUSTER MAP -’.16, 

* x.SX, ’ TOTAL NO. 0^ CLUSTERS -’,16, 

% /, 10X, ’ CLUSTER CLUSTER NO. OF NO. OF DOTS’ 

* /, 10X, ’ ISO. CODE PIXELS ASSIGNED’) 

DO 283 K-l.MX 

IF(JFLAG.EQ.l) URITE(6.284)K.LL(K) ,MM(K) .NISCK) 

283 CONTINUE 

284 FORMAT ( 13X, 12. 5X. 14, 3X. 16, 4X, 13) 

IF(MX.EQ.M) GOTO 288 

KX-0 

DO 285 K-l.M 

IFCNN(K) .NE.0) GOTO 285 

KX-KX+1 

NDARY(KX) -K 

285 CONTINUE 

I F ( JFLAG . EQ . 1 ) WRITE(6. 286)MX, (NDARYCK) , K-l . KX) 

286 FORMAT!/’ L THOSE SMALLER CLUSTERS GROLIPED (AS CLUSTER’ , 13, 

* ’ ) ARE:’ , 3(x, 16(13, ’ , ’ J ) ) 

288 CONTINUE 

C 

C 
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C END OF PROPORTIONAL DOT ALLOC l AT I ON 

C 

C 

C FOR EACH CLUSTER. ASSIGN DOTS 

ND-1 

DO 291 K“1 1 MX 
Kl-NNCK) 

IFCK1.EQ.0) GOTO 291 
DO 289 L-1.K1 

289 CALL 6ETD0T(MM(K).L.NDARYCND),K3) 

291 ND-ND+Ki 
C 

c FOR EACH CLUSTER. FIND LINE AND COL I'd. OF EACH ASSIGNED DOTS 

IFCJFLAG.EQ. 1) URITEC6.295) 

295 FORMAT (//’ DOT FILE ’ 

* /’ CLUSTER DOT NUMBER POSITION GROUND TRUTH LABEL’ 

X /’ NO. URT ITS CLUSTER LINE. COL RAW, CODE’) 

PSG-0. 

ND-1 

DO 331 K-l.MX 
Kl-NNCK) 

1FCK1 . EQ.0) GOTO 331 
ISG-0 

K3-NDARYCND) 


c T0 TEST FOR THE EXTRA CLUSTER 

IF(K.LE.M) GOTO 311 

c AL.b. TO GET CORRECT LBLCND) FOR THE EXTRA CLUSTER 


DO 305 J-l.M 
IF(NN(J).NE.0) GOTO 305 
IFCK3.LE.MMC J) ) GOTO 30? 

K3-K3-MMCJ) 

305 CONTINUE 

.307 CALL CLMPXYCJ.K3.NDL.NDC.LL) 

CALL GTMPLB ( NDL . NDC .LBLCND) ) 

CALL LBLITPCLBLCND) , LB1 . ISG) 

IFCJFLAG.EQ. 1) WRITEC6.308 j J. K3.NDL.NDC.LBL.CND) .LF1 
308 FORMAT (3X, 12, 8X, 16, 7X, 13, 14, 5X, 14, IS) 

GOTO 321 

311. CALL CLMPXY C K , K3 , NDL , NDC , LL. ) 

CALL GTMPLB C NDL . NDC , LBL C ND ) ) 

CALL LBL I TP ( LBL ( ND 1 ,LB1, ISG) 

IFCJFLAG.EQ. 1) URITE(6,308)K.K3. NDL, NDC, LBLCND) ,LB1 
321 ND-ND+1 
331 CONTINUE 
C 

C COMPUTE MAJORITY LABEL FOR EACH CLUSTER 

IFCJFLAG.EQ. 1) URITEC6.421) 

421 FORMAT C///’ MAJORITY LABELLING RULE BY GROUND TRUTHS’/ 
X ’ CLUSTER NO. OF LABELS CHOSEN LABEL’ ) 

ND-1 
P5G-0. 

DO 471 K-l.MX 

I F C NN C K ) . EQ . 0 ) GOTO 471 

CALL MR C LBL C ND ) . NN C K ) . LABEL , J 1 . JG , LG ) 

ISG-0 

CALL LBL I TP C LABEL. , LB1 > ISG > 

IFC ISG. EQ. 1 ) PSG= PSG+FLOAT C MM C K ) )/NPIXEL 
I F C JFLAG . EQ . 1 ) HR I TE C 6 . 451 ) K . J1 . LABEL , L. B1 
451 FORMAT C GX , 13, 9X, 13, 3X, 14, IX, 14) 

ND-ND+NNCK) 

471 CONTINUE 
GOTO 990 
C 

901 UR I TE C 6 . 902 ) M 

302 FORMAT C’ YOU HAUE’ , 13, ’ TOO MANY CLUSTERS (MAX -50) ’ ) 
GOTO 990 

911 UR I TE ( 6 » 91 2 ) NDOT 

912 FORMAT C’ YOU SPECIFIED’ , 14, ’ TOO MANY DOTS C MAX -200) ’ ) 

C 

990 RETURN 
END 




7. A84 ; BAYESIAN DOT ALLOCATION (.UNIFORM PRIOR) 


7.1 DESCRIPTION 

Implemented in this piece of software is scheme in which dot 
allocation is sequentially Bayesian and labeling is direct from 
ground truth. 


In the Bayesian dot allocation scheme , dots are allocated sequen- 
tially to those clusters that may reduce the variance of the seg- 
ment proportion estimate the most. In order to do this, we must 

2 

first compute the expected change in variance. La , for each 
cluster. The expected change in variance for cluster i is 
defined as 


Ao.' 


■w 


'( ni + 5) 2 (n. 2 + 7n L + 8) 


( ni - l)n i (n i + 2) (n*. + 3)' 


+ 1) (n ± 


x ± + 1) 


(7-1) 


where N^ = number of pixels in cluster i 

N = total number of pixels in the entire scene 

n^ = number of dots previously allocated to cluster i 

x. = number of dots previously allocated to cluster i 
which are labeled as small grain 

(notice x^ < n^) 

2 . 

Next, a dot is allocated to the cluster whose La ^ is the largest. 
Then, for this chosen cluster, n i is updated to n. + 1, That 
dot's label is read from the ground truth file. If the label is 

small grain, x. is updated to x. +1. After Miis, the segment 

2 1 1 
variance, a , is computed as 





4t. 


(7-2) 


where m » total number of clusters 


P. 

1 


+ 1 
n i + 2 


(7-3) 


2 

This a is printed and the above allocation scheme is repeated 
until the number of dots allocated reaches a value specified by 
the user. 


The proportion estimate is computed as 






(7-4) 


This scheme, like other allocation and labeling schemes, may be 
repeated with different pseudorandom sequences. The main program 
resets the random number generator to a different starting point 
for each repetition. 


All computation of statistics such as mean and variance is per- 
formed in the main program. A grand summary is printed. 


7.2 LINKAGES 

The POP 11/45 Fortran library (for functions and subroutines such 
as SQRT) , the PDP 11/45 system library (for routines such as the 
pseudorandom number generator) , and the utility package (described 
in section 14) are required. 


Main program 

or subprogram Subroutine or function required 


MAIN 

SAEB1 

GETDOT 


SAEB1, GETDOT, CLMPCS, GTMPLB, LBLITP 

GETDOT, CLMPCS, GTMPLB, LBLITP 

RAN 



7.3 INTERFACE 


Interface with other routines is through the common block PRTFLG 
(2 bytes) , which is used to control the optional printing of dot 
files and other information. 

7.4 INPUTS 

Fortran- formatted input of the following parameters is needed: 


Card 

or 

line 


Parameters 

Format 

Default 

input 

Default 

value 

1 

Cluster map file name 

A13 

None 


2 

a. 

Number of repetitions 

13 

0 or blank 

1 


b. 

Starting point of first 
pseudorandom sequence 

15 

0 or blank 

10 


c. 

Number of repetition 
printings 

13 

0 or blank 

5 

3 

a. 

Maximum number of dots 

13 

None 




that can be allocated 





b. 

Number of dots initially 
assigned to each cluster 

13 

None 


4 

Number of status messages on 
terminal 

13 

0 or blank 

NO MSG 


Two disk files are needed, one containing the stripped cluster 
map, the other the stripped ground truth map. These files should 
have been created by conversion programs A81 and A86. 
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7.5 OUTPUTS 


In addition to the normal echo of input specifications and status 
messages on the terminal, the following will be output on a line 
printer. 


Report 


Contents 


Individual repetition 1. Processor header 

2. Ground truth small-grain proportion 

3. A dot file showing dots chosen, their 
x-y position, their ground truth labels, 
segment variance, and proportion 
estimate 


Grand summary 


4. A table showing number of clusters, 
cluster name codes, cluster sizes, 
and number of dots assigned 

1. Processor header 

2. A table showing bias, MSE, reduction 
in MSE, average, variance, and vari- 
ance reduction as functions of number 
of dots 


7.6 BRIEF ALGORITHM: FOR BAYESIAN DOT ALLOCATION 

(UNIFORM PRIOR) 


1. Type messages at terminal and read from it run specifications. 


2. Compute ground truth small-grain proportion. 

3. Repeat a through g until all repetitions are finished: 


a. Set starting point of the pseudorandom number generator. 

b. Clear a dot counter. 

c. Increment the dot counter. 


d. If the dot counter indicatea that the current dot should 

be allocated aa an initial dot, then go to the next step, 

2 

Otherwise compute Ao< for every cluster and choose the 

x 2 

cluster with maximum Ao^ . 


e. Allocate a dot -o the chosen cluster. 

f. Pick up the small-grain labels, compute segment variance 
and proportion estimate, and print a line of the dot file 


g. If the dot count does not exceed specified value, go to c 


Compute bias, MSE, reduction in MSE, average, variance, and 
variance reduction, and print a grand summary. 






Figure 7-1.— Flow chart for Bayesian dot allocation 


(uniform prior) . 


2* 




7.7 LISTING 


C 

c 


C — PROGRAM A09> PROPORTION ESTIMATE OF SMALL GRAIN 

C- USING BAYES (UNIFORM PRIOR) DOT ALLOCATION 

C LABELLING IS DIRECT BY GROUND TRUTH. 


C 

BYTE NAME( 15) »NGT(13) 

REAL PSG(20O)«PM(200).PB(200) 

C COMMON BLOCK IS A FLAG FOR PRINT ENABLE/D I SABLE 

COMMON ^PRTFLG/JFLAG 

DATA NGT(5)/’ . ’/.NGT(6)/’S’/,N(aT(?)/’T’/.Nl3Tce)/’P’/ 

C 

C 

C READING INPUT FROM TERMINAL 

C 

URITE(8> 121) 

121 FORMAT ( ’ PROGRAM* AGS.’ 

t /’ PROPORTIONAL ESTIMATION OF SMALL GRAIN' 

* /’ DOT ALLOCATION IS BAYES (UNIFORM PRIOR) .* 

* /’ LABELl.K'tG IS DIRECT BY GROUND TRUTH.' 

* INPUT CLUSTER MAP FILENAME'/' AAAAAAAAAAAAA’ ) 

READC7. 122) CNAME(K) . K-l . 13) 

122 FORMAT (13A1) 

,r 

UR I TECS. 131 ) 

131 FORMAT ( ' AN ASSIGNMENT OF PSEUDORANDOM DOTS IS A REPETITION RUN’ 

* /• HOU MANY REPETITION RUNS? SEED FOR FIRST PSEUDO SEQ?’ 

* /’HOI MANY PRINTINGS OF INDIUIDUAL RUNS?' 

* /’ III IIIII III’) 

READ<7. 132) JU. JSEED. JPAGE 

132 FORMAT C 13. IX. IS. IX. 13) 

IFCJU.LE.0) JU-1 

IFC JSEED. LE.0) JSEED- 1© 

IF(JPAGE.lE.O) JFAGE-S 
JSKIP- ( JU-1 )/JPA6E+l 
C 

UPITEC8. 141 ) 

141 FORMAT C ' SPECIFY MAXIMUM NO. OF DOTS IN THE REPETITIONS RUNS. ’/ 

* ' AND M3. OF DOTS I NT I ALLY ASSIGNED TO EACH CLUSTER. */ 

* ’mill’) 

READ C 7. 142)MAXD0T.NNIX 

142 FORMAT ( 13. IX. 13) 

CALL INITl(HNIX) 

C 

WRITEC8.1S1) 

151 FORMAT ( 1 STATUS MESSAGES ON TERMINAL. HOU MANY?’/’ Ill’) 

READC7, 152)M1ES 

152 FORMAT (13) 

C 

URITE(8. 181) CNAME(K) .K-l. 13). JU. JSEED. MAXDOT 
181 FORMAT (//’ CLUSTER FILE GIUEN - *.13A1. 


* /’ NO. OF REPETITION PUNS FOR EACH TOTDOT •’ . 13. 

* /’ THE FIRST REPETITION RUN STARTS UITH SEED • ’.16. 

* /’ MAXIMUM NO. OF DOTS IN THE REPETITION RUNS-’, 16) 
C 

C COMPUTE BAYES (UNIFORM PRIOR) ESTIMATE FOR GROUND TRUTH 

NLINE-117 


NCOL-196 
DO 21S K-l. 4 
215 NGT(K)-NAMECK) 

OPEN (UNIT*??. NAME-NGT . TYPE- ’ OLD ’ . READONLY . FORM- ’ UNFORMATTED ’ . 
* ACCESS-’DIRECT’ ) 

DO 231 L-l.NLINE 
DO 221 K-l.NCa 
CALL GTMPLBO-.K. LABEL) 

CALI LBLITP(LABEL,LB1. IP) 

221 CONTINUE 

231 CONTINUE 

P- FLOAT C I P ) /M_ I NE/NCOL 


c 


OPEN ( UN 1 T- 1 . NA ME - NAME. TYPE -’ OLD* , FORM- - UNFORMATTED ’ , 
* READONLY. ACCESS- ’DIRECT’) 

CALI. CLMPLC ( NL 1 NE » NCOL ) 


C 

c 

311 

321 


322 


383 
C 

341 

'361 

C 

C 

371 


301 

382 

C 

c 

401 

C 


- START ESTIMATION FOR EACH TOTDOT AND REPETITIONS 
M 311 K-l. MAXDOT 

PB(K)-©. 

PMCK) -O. 

JMES-0 
JS-* JSEED 
DO 361 J-l.JU 
JFL.AG-0 

I F ( MOD < J . JSK 1 P ) . EQ . 1 ) JFLAG-1 
IFC.JFLAG.EQ.l) WRITEC6. 321 )NAME, P 

FORMAT i*l’* 10X, ’ ACCURACY ACCESSEMENT SOFTWARE ( S-9-7S ) ’ 

* x,2X. ’PROGRAM A09> ESTIMATION OF SMALL 588 l NS. PRO PORT I ON’ 

* /. 1SX. ’DOT ALLOCATION • BAYES (UNIFORM PRIOR) ,’ 

* /.ISX.’DOT LABELLING 'DIRECT BY GROUND TRUTH.' 

* //. 10X, ’ INPUT CLUSTER MAP IS FILE ’»15A1» 

* .10X, ’ PROPORTIONAL ESTIMATE FOR GROUND TRUTH IS'.FS.S) 
IF( JFLAG.EQ. 1) WRITE<6. 322) J. JSEED 

FORMAT C/.2X. * REPETITION RUN-'. 13, 

* ’ RANDOM DOT SEQUENCE STARTS WITH’ , 18. ’ — -* ) 

CALL PANST (JSEED) 

CALL SAEB1 ( MAXDGT . PSG ) 

JMES-JMES+1 

IF(JMES.LE.NMES) UR I TE ( 8 , 323 ) MAXDOT . J . JSEED 

FORMAT ( ’ MAXDOT-*, J3.’ REPETITION RUN-*, 13.’ SEED-*. IS) 

JSEED- JSEED+150 

— COMPUTE BIAS AND M.S.E. 

DO 341 K-l, MAXDOT 
TEr.P-PSG(K)~P 

PBiK ) -PB(K) i TEMP 
PM ( K ) - PM ( K ) +TEMP#*2 
com I NUE 


- PRINT GRAND SUMMARY FOR THIS JOB 
WRITE(6.321)NAME,P 
WR I TEC6. 371 ) MAXDOT, JU, JS 

FORMAT'./. 14X. * GRAND SUMMARY OF THIS JOB ’, 

* '.SX,’ MAX. NO. OF DOTS IN EACH REPETITION RUN-’, 13, 

* /.14X.’ NO. OF REPETITION RUNS 13, 

* 14X, ’ RANDOM DOTS START WITH SEED-’.I6, 

* '/.IX.’ DOT BIAS M.S.E. REDUCTION’. 

*. ’ AUERAGE UARIANCE UAR REDUCTION’ ) 

DO 331 K-l. MAXDOT 
PB(K)-PB<K)XJU 
PMCK) -PMC JO/JU 
RR-PHC K >/( P* ( 1 . -P ) y'K ) 

AUERG»PBtK)+P 

IFc JU.EQ. 1) UAR- PM(K)-PB(K)**2 
IF(JU.GT.l) MAR- (PM(K)-re(K)JM^)'4;JUX(JU-l) 

REDUAR-UAR/( P*(1.~P)/K ) 

WR I T E ( 6 . 382 ) K . PB ( K ) . PM ( K ) . RP . AUERG . UAR . REDUAR 
FORMAT ( 2X » 13, 2X.F9.S, 4CF10.6), 3X.F10.6) 


WRITE (6. 401) 

FORMAT ( ’ 1 END OF THIS JOB ’ ) 

STOP 

END 


C 

C 

C 

SUBROUTINE SAEBlt MAXDOT. PSG) 

C STRATIFIED AREAL ESTIMATION USING BAYES DOT ALLOCATION. 

C FORMULA HERE ARE FOR BAYES WITH UNIFORM PRIOR. 

C GIUEN ANY MAX NO. OF DOTS IN MAXDOT. PROGRAM WILL RETURN 

C THE ESTIMATES ARE IN ARRAY PSG ON RETURN. 




c 

C UR 1 TEN AND EDITED BY N.Y. CHU ON 5-9-79. 

C 

INTEGER NDAPY(20©>tMM(Sl)*NN(51)«NX(Sl)tLL<Sl) 

REAL UU.PSG(l) 

BYTE LABEL 

C COMMON BLOCK IS A FLAG FOR PRINT ENABLED I SABLE 

COMMON /PRTFL6/JFLAG 
DATA NSTART/0/ , MINSZ/S/ 

DATA INI DOT/3/ 

C 

1F(MAXDOT.GT.SQ0) GOTO 911 
C 

C GET PIXEL COUNT FOR EACH CLUSTER 

IFCNSTART.NE.C) GOTO 121 
NSTART-1 VfrO fa 

CALL CLMPCSCNPIXEL.MiMM.LL) 

IFCM.GE.51 ) GOTO 901 ** 

121 CONTINUE 
C 

C BEGIN RAVES DOT AU OCATIlTN 

IFC JFLwG. EO. 1 ) URITEC6.141) 

141 KOK1AT ( // « 25X. * DOT FILE ’ 

* /’DOT CLUSTER INFORMATION POSITION GROUND TRUTH* 

* . ! SEGMENT PROPORTIONAL* 

* /* NUMBER NO. PIXEL DOTS S.G. LINE. COL PAL. CODE* 

■ * .* UARIANCE ESTIMATE*) 

C 

C MAKE SURE EACH CLUSTER HAS A MINIMUM NO. OF PIXELS 

Kl-0 

NDS»0 

DO 1S1 K-l.M 

NXCK)-0 

HNCKJ-0 

IF(MMCK) . GE.MINSZ) NDS-NDS+INIDOT 
IF(MMCK) .LT.MINS2) K1-K1+MMCK) 

IFCI71CIO . LT.MJNS2) NNOO — 1 
151 CONTINUE 
MX*M 

IFCK1.LT. MINS2) GOTO 161 
MX-M+1 
MMCMX) »K1 
NDS«NDS+INIDOT 
161 CONTINUE 
C 

C FROM HERE TO STATEMENT 261* LOOP TO ASSIGN 1 DOT AT A TIME 

C 

C IN ORDER TO ACHJEUE ASSIGNING 3 DOTS INITIALLY OR LESS 

C TO EACH CLUSTER. 

c THE DELTA UARIANCE COMPUTATION IS BYPASSED FOR THE FIRST 

C NDS-INIDOTYMX DOTS 

C 

KG1-0 

NDOT-0 

191 NDOT«NDOT+1 

IFCNDOT.GT.NDS) GOTO 200 
195 KG1-KG1+1 

IF(KGl.GT.MX) KG1*1 
K1-KG1 

IFCNNCKGl ) .GT.-i ) GOTO 214 
GOTO 195 
C 

c 

C FOP EACH CLUSTEP. COMPUTE DEl.TA UARIANCE AND CHOOSE THE LARGEST 

200 AMAX>0. 

DO 211 K-l.MX 
J»NN(K) 

IFCJ.LE.l) GOTO 211 

UAF- C FLdATCMM(10)/NPIXEL*(.J+5) >**2 






* *(FL0ATCJ)«j!7.fJ+8. )*(NXcK)+l. )4t(J-f«(K)+l. ) 

* /( ( J-l. )*J*(J+2. >#*2 ) 

C WRITE<6. 99121 )K.MMCK). J.NXOO.UAR 

99121 FORMAT ( ’ K.MMOO . J.fWlO.UAR-S 14. 18.21S.E16.6) 

IF(UAR.LT.AMAX) GOTO 211 
K1*K 

AMAX-UAR 

2ii corn 1 nue 

r, 

L - ASSIGN A DOT TO THE CLUSTER WITH MAX DELTA UAR1ANCE 

IFiNNLKlJ.K.MMtKin WRITE(6.213)ND0T.K1.MM(K1> 
IF(NN<Kl).GE.MMlKl> ) GOTO 2-41 

213 FORMAT ( * ATTEMPT TO ASSIGN DOT- SI 3. ’ WHILE NO PIXELS 

* ’ IS AVAILABLE. CLUSTER-SI3.’ TOT PIXEL- SI6) 

C POSITION THE POINTER TO NDARY 

214 ND»1 

DO 215 K-1.K1 
IFiNN(K).LE.C) GOTO 215 
ND-ND+NNOO 

215 CONTINUE 

C- MAKE ROOM FOR THE NEW DOT IN NDARY 

Nl— 1 

nn ?i*? ♦'^wrjnr.NP N) 

cit 

I-T'-ND-NN(K1 ) 
tf4<Kl»-NN(Kl)+l 
C ASSIGN A DOT 

CALL GETDGT(MM(K1).NN(K1>.NDARY<ND>.ID0T> 

«• WRIT£(6.99141)K1.HD.NN(K1).NX(K1). (NDARY<lO.K-1.44) 

99141 FORMAT C K1.ND-NN.NX-S4I4. 4(Y,2X. 1514) ) 

IFCH.LE.M) GOTO 227 
DO 225 Ki-l.M 
IF'-NNCKD.GT.C) GOTO 225 
IFi iror.LT.IlKKi.' ■ GOTO 227 
IDOT-IPOT-MMOC1) 

225 CONTINUE 

227 C hLL CLMPXY iKi.IDOT.NL.NC.LL> 

CALL GTMPLB' NL.NC. LABEL > 

CAl L LBL I TP t LABEL . LB1 . NX ( K1 ) ) 

PS&iNDOTJ-C*. 

DO 229 JX-i.MX 
lF(NN(JVO.LE.0j GOTO 229 

PSG < NDOT i * PSG ( NDOT > + FLOAT U « C JX) >YNN(.JX) *MM ( JX > /NP I XEL 
229 CONTINUE 

C 

c 

C- CONFUTE SEGMENT UAPIANCE 

241 UU» 0. 

DO 251 K-l.MX 
IF'.NNUO.LE. 1) GOTO 251 
P* <NX(K) + 1. J/'tNNtF )+2. ) 

UU-UiJ+fFLOAUMMiK) i/TIPIXEL )**2 *P*(1.-P-Y( NN<K)-1. > 

251 CONTINUE 

IF( JFl-AG. EQ. 1 1 WR I TE 1 6 • 255 ) r (DOT • K.i « IDOT.NNiKl ) .NXCK1 ) 
t , NL. NC. LABEL. LB l.L*',PS«5t NDOT) 

255 FORMAT 1 2X. 13. 3: 13. IX. IS- IX. IJ.2X. 13. 2X.2I4. 2X.214. 2F11.6) 

IFiNDOT.LT.MAXDOT) GOTO ;91 

C 

DO 261 1*1. 11>. 

! F ( HH > F i . LE*: -i ) HN « K ) -0 
261 CONTIfiUF 


C- - END OF BAYES TOT ALLOCATION 



J 


°f the 

pr) n n 


c 

c 

C WRITE CLUSTER INFORMATION 

IF(JFLAG.EQ.l) WRITEC6, 282)MAXD0T,NPIXEL,MX 
£82 FORMAT C/ 1 / 1 5X> * TOTAL NO. OD DOTS AUAILABLE •’ . 16. 

* x.SX.’ TOTAL NO. OF PIXEL IN THE CLUSTER MAP «’ , 16. 

* /.SX,’ TOTAL NO. OF CLUSTERS -’.IS, 

* /. 5X> ’ CLUSTER CLUSTER NO. UF NO. OF DOTS LABELLED* 

* /,SX. ’ NO. CODE PIXELS ASSIGNED SMALL GRAIN’) 

DO 293 K-l.MX 

IF(JFLAG.EQ.l) URITE(6,284)K.LL(IO ,MM(K) ,W(K) ,NX(K) 

283 CONTINUE 

284 F0RMATC8X, I2.SX, I4.3X, I6.4X, 13, 10X, 13) 

IF(MX.EQ.M) GOTO 288 

KX-0 

DO 285 K-l.M 

IFCNN(K)'. NE. 0) GOTO 285 

KX-KX+i 

NDARV(KX)-K 

285 CONTINUE 

IFCJFLAG.EQ.l) WRITEC6, 286)MX, (NDARY(K) ,K-i,KX) 

286 FORMAT (/’ THOSE SMALLER CLUSTERS GROUPED (AS CLUSTER’ . 13, 

X ’ ) ARE:’ , 3C/.16< 13, ’ , ’ ) ) ) 

288 CONTINUE 

C 

C 

C 

C 

GOTO 990 
C 

901 WRITE (G, 902 )H 

902 FORMA 1 ( ’ YOU HAUE’ , 13, ’ TOO MANY CLUSTERS (MAX=50) ’ ) 

GOTO 990 

911 WRITE(6, 912)ND0T 

912 FORMAT (’ YOU SPECIFIED’ , 14, ’ TOO MAN-' DOTS (MAX =200) ’ ) 

C 

C ENTRY TO DEFINE I NT I AL NO. OF DOTS TO BE ASSIGNED 

ENTRY IHITl(NNIT) 

INIDOT-NNIT 
990 RETURN 
END 


8. A85: BAYESIAN DOT ALLOCATION (NO PRIOR) 

8.1 DESCRIPTION 

Implemented In this piece of software is a scheme in which dot 

allocation is sequentially Bayesian (no prior) and labeling is 

direct from ground truth. This scheme uses the same algorithm 

as does A84 (Bayesian dot allocation, uniform prior) but has 

2 

slightly different formulas for Acr^ and segment variance, as 
will be noted below. 


In the Bayesian dot allocation scheme, dots are allocated sequen- 
tially to those clusters that may reduce the variance of the seg- 
ment proportion estimate the most. In order to do this, we must 
first compute the expected change in variance, A o , for each 
cluster. The expected change in variance for cluster i is 
defined as 


? / N -V 
Ao i = (fr) 


V n i - x i ) 

Tn ± - l) J n. 2 (n i + l) 2 


(8-1) 


where N^ = number of pixels in cluster i 

N = total number of pixels in the entire scene 

n^ = number of dots previously allocated to cluster i 

x i = number of dots previously allocated to cluster i 
which are labeled as small grain 


(notice x^ < n^) 


2 

Next, a dot is allocated to the cluster whose Aa i is the largest. 
Then, for this chosen cluster, n^ is updated to n^ + 1. That 
dot's label is read from the ground truth file. If the label is 



a small grain , x. is updated to x 4 + 1. After this, the segment 
2 i 1 

variance, o , is computed as 


o 


2 




(8-2) 


where m = total number of clusters 


P. 

x 



(8-3) 


2 

This a is printed and the above allocation scheme is repeated 
until the number of dots allocated reaches a value specified by 
the user. 


The proportion estimate is computed as 


P 

sg 





(8-4 ) 


This scheme, like other allocation and labeling schemes, may be 
repeated with different pseudorandom sequences. The main program 
resets the random number generator to a different starting point 
for each repetition. 

All computation of statistics such as mean and variance is per- 
formed in the main program. A grand summary is printed. 

8.2 LINKAGES 

The PDP 11/45 Fortran library (for functions and subroutines such 
as SQRT) , the PDP 11/45 system library (for routines such as the 



pseudorandom number generator), and the utility package (described 
in section 14) are required. 


Main program 

or subprogram Subroutine or function required 


MAIN SAEB2, GETDOT, CLMPCS, GTMPLB, LBLITP 

SAEB2 GETDOT, CLMPCS, GTMPLB, LBLITP 

GETDOT RAN 


8.3 INTERFACE 

Interface with other routines is through the common block PRTFLG 
(2 bytes) , which is used to control the optional printing of dot 
files and other information. 


8.4 INPUTS 

Fortran-formatted input of the following parameters is needed: 
Card 

or Default Default 


line 


Parameters 

Format 

input 

value 

1 

Cluster map file name 

A13 

None 


2 

a. 

Number of repetitions 

13 

0 or blank 

1 


b. 

Starting point of first 
pseudorandom sequence 

15 

0 or blank 

10 


c. 

Number of repetition 
printings 

13 

0 or blank 

5 

3 

a. 

Maximum number of dots 
that can be allocated 

13 

None 



b. 

Number of dots initially 
assigned to each cluster 

13 

None 



4. 

Number of status messages 
on terminal 

13 

0 or blank 

NO MSG 




Two disk files are needed/ one containing the stripped cluster 
map, the other the stripped ground truth map. These files should 
have been created by conversion programs A81 and A86. 


8.5 OUTPUTS 

In addition to the normal echo of input specifications and status 
messages on the terminal, the following will be output on a line 
printer. 


Report 


Contents 


Individual repetition 1. Processor header 

2. Ground truth small-grain proportion 

3. A dot file showing dots chosen, their 
x-y position, their ground truth labels, 
sagment variance, and proportion 
estimate 


Grand summary 


4. A table showing number of clusters, 
cluster name codes, cluster sizes, 
and number of dots assigned 

1. Processor header 

2. A table showing bias, MSE, reduction 
in MSE, average, variance, and vari- 
ance reduction as functions of number 
of dots 


8.6 BRIEF ALGORITHM: FOR BAYESIAN DOT ALLOCATION (NO PRIOR) 

1. Type messages at terminal and read from it run specifications. 

2. Compute ground truth small-grain proportion. 

3. Repeat a through g until all repetitions are finished: 

a. Set starting point of the pseudorandom number generator. 

b. Clear a dot counter . 



c 


Increment the dot counter 


d. If the dot counter indicates that the current dot should 

be allocated as an initial dot, then go to the next step. 

2 

Otherwise compute Ao. for every cluster and choose the 

1 2 

cluster with maximum Ao^ . 

e. Allocate a dot to the chosen cluster. 

f. Pick up the small-grain labels, compute segment variance 
and proportion estimate, and print a line of the dot file 

g. If the dot count does not exceed specified value, go to c 

Compute bias, MSE, reduction in MSE, average, variance, and 
variance reduction, and print a grand summary. 


& 



Figure 8-1.- Flow chart for Bayesian dot allocation (no prior) 
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8.7 LISTING 


r, 

c • 

C PROGRAM A12(A85.TSK) « PROPORTION ESTIMATE OF SMALL GRAIN 


C USING BAYES (NO PRIOR) DOT ALLOCATION 

C LABELLING IS DIRECT BY GROUND TRUTH. 

C 

BYTE NAT1EC 15) , NGTC 13) 

REAL PSG(200).PM(200),PB(200> 

C COMMON BLOCK IS A FL.AG FOR PRINT ENABLE/DISABLE 


COMMON /PRTFLG/JFLAG 

DATA NGT(S)/’ . * /.M3T(6)/’S’/»N6T(7)/*T’/,M3T(8)/’P’/ 

C 

C 

C READING INPUT FROM TERMINAL 

C 

WRITE(8,121) 

121 FORMAT ( ’ PROGRAM* A12CA8S. TSK) . ’ 

* /’ PROPORTIONAL ESTIMATION OF SMALL GRAIN* 

* /’ DOT ALLOCATION IS BAYES (NO PRIOR) . ’ 

* /’ LABELLING IS DIRECT BY GROUND TRUTH.’ 

X /’ INPUT CLUSTER MAP FILENAME’/’ AAAAAAAAAAAAA’ ) 

READ(7,122)(NAME(K).K-1.13) 

122 FORMAT ( 13A1 ) 

C 

WRITE(8, 131) 

131 FORMAT ( ’ AN ASSIGNMENT OF PSEUDORANDOM DOTS IS A REPETITION RUN’ 

* /’ HOW MANY REPETITION RUSS? SEED FOR FIRST PSEUDO GEO'?’ 

X /’ HOW MANY PRINTINGS OF INDIUIDUAL RUNS'?’ 

x /’ in iiiii Tiro 

READ ( 7 > 1 32 ) JU , JSEED . JPAGE 

132 FORMAT (13. IX, IS. IX. 13) 

IFCJU.EO.0) JU«1 

I F ( JSEED . EQ . 0 ) JSEED- 10 
I F ( JPAGE . EQ . 0 ) JPAGE -5 
JSK I P« ( JU-1 ) /JPAGE+1 
C 

WR1TE(8,141) 

141 FORMAT (’ SPECIFY MAXIMUM NO. OF DOTS IN THE REPETITIONS PUNS.’/ 

X ’ AND NO. OF DOTS I NT I ALLY ASSIGNED TO EACH CLUSTER.’.' 

X * III III’ ) 

READ ( 7 , 142)MAXD0T,NN1X 

142 FORMAT ( 13, IX. 13) 

CALL INITS(NNIX) 

C 

WR I TE ( 8 , 151) 

151 FORMAT (’ STATUS MESSAGES ON TERMINAL, HOW MANY?’/’ Ill’) 

READ(7, 1S2)NMES 

152 FORMAT (13) 

C 

WR I TE ( 8 , 1 81 ) ( NAME (K) ,K*1 . 13) , JU, JSEED > MAXDCT 
181 FORMAT (//’ CLUSTER FILE GIUEN = ’.13A1, 

X /’ NO. OF REPETITION RUNS FOR EACH TfJTDOT -’.13, 

X /’ THE FIRST REPETITION RUN STARTS WITH SEED - ’,16, 

X /’ MAXIMUM HO. OF DOTS IN THE REPETITION RUNS-’, 16) 

C 

C COMPUTE BAYES (NO PRIOR) ESTIMATE FOR GROUND TRUTH 

NLINE-117 
NCOL-196 
DO 2 i 5 K-1,4 
215 NGT(K) -NAME(K) 

OPEN ( UN I T - 2 . NAME - NGT . TYPE - ’ OL D ’ , READONLY , FOR1 1 - ' UNFORMATTED ’ , 

X ACCESS- ’ DIRECT’ ) 

UU iJi L-l.NLlNE 
DO 221 K-l.NCOL 
CALL GTMPLBCL.K. LABEL) 

■ CALL LBLITP(LABEL, LB1 , IP) 

221 CONTINUE 

231 CONTINUE 


F-FLQhT< lF^NLINE.'NCOI. 


OPEN C UN I T ■ 1 > NOME ■ NAME » TYPE ■’ OLD ‘ » FORM ■ ’ UNFORMATTED ’ , 

* READONLY . ACCESS- ’ D I RECT ’ ) 

CALL CLMPLCCNLINE.NCOL) 

C START ESTIMATION FOR EACH TOT DOT AND REPETITIONS 

DO 311 K-l, MAXDOT 
PB(K) *0. 

311 pmoo-s. 

■ JME3-0 
JS-JSEED 
DO 361 J-l.JU 
JFLAG-0 

IF(MOD( J, JSKIP) .EQ. 1 ) JFl.AG-1 
I F i JFLAG . EQ . 1 ) WRITE (6. 32 i ) NAME . P 

321 FORMAT ( ’ 1 ’ . 10X, ’ ACCURACY ACCESSEMENT SOFTWARE (5-1 0-79 ) ’ 

* /. 2::. ’PROGRAM Al2i ESTIMATION OF SMALL GRAINS PROPORTION’ 

* /.15X, ’DOT ALLOCATION i BAYES (NO PRIOR) , ’ 

* / . 1SX , ’ DOT LABELL I NG « D I RECT BY GROUND TRUTH . ’ 

* //. 10X, ’ INPUT CLUSTER MAP IS FILE ’.1SA1, 

* 10X, ’ PROPORTIONAL ESTIMATE FOR GROUND TRUTH IS’,F8.5) 
IFcJFLAG.EQ.i) WRITEC6. 322) J, JSEED 

322 FORMATS, 2X. ’ REPETITION RUN- ’, 13, 

* ’ RANDOM DOT SEQUENCE STARTS WITH’ , 10. ’ ’ ) 

CALL RANSTC JSEED) 

CALL 5A£B2 < MAXDOT , PSG ) 

JMES-JMES+1 

I F < JMES . LE . NMES 1 WR I TE ( 8 , 323 ) MAXDOT . J , JSEED 
223 FORMAT k’ MAXDOT* *, 13, ’ REPETITION RUN-’ , 13, ' SEED-’, IS) 

JSEED -JSEED+150 

C COMPUTE BIAS AND M.S.E. 

DO 341 K-l, MAXDOT 
TEMP-PSG(K)-P 

pboo-pbckj+temp 

341 PM CIO - PI-1 ( K ) +TEMPYY2 

361 CONTINUE 
C 

C PRINT GRAND SUMMARY FOR THIS JOB 

WR I TE ( 6 > 321 ) NAME , P 

WR I TE ( 6 , 371. > MAXDOT . JU , JS 

371 FORMAT ( o 14X, ’ GRAND SUMMARY OF THIS JOB 

* ' « 5X . ’ MAX. NO. OF DOTS IN EACH REPETITION RUN- ‘.13, 

t s . 14X, ’ NO. OF REPETITION RUNS -’.13, 

* -'r 14X, ’ RANDOM DOTS START WITH SEED- ’,16, 

% / -• , IX , ’ DOT BIAS M.S.E. REDUCT I ON ’ , 

* ’ AVERAGE VARIANCE VAR REDUCTION’ ) 

DO 331 K-l. MAXDOT 

PBCK)«PB(K>''JU 
FM'.K)*PM<:K).- JU 
RP*PM( K ' • C P*(1.-R)'K > 

AUERG-PBiK)+P 

IF(JU.EQ.l) VAR- FM(K)-PBCK)**2 
I F ( JU . GT . 1 ) VAR- ( PM < K) -PB C K) Yl:£ ) *JU/ C JU-1 ) 

REDUAR-VAR/C P*<l.-P)/K ) 

381 WP I TE < 6 • 382 ) K . PB (K ) , PM ( IO , PR , AVERG , UAR . REDUAR 

332 FORMAT C2X. 13. 2X.F9.5, 4CF10.6), 3X.F10.6) 
r; 

C 

WPITEC6, 401 ) 

401 FOFMhT C ’ 1 END OF THIS JOB — - ’ ) 

C 

STOP 

END 

C 

C 

r 

SL ©POUT I HE SAEB2 1 MAXDOT , PSG ) 


V STRATIFIED AREAL ESTIMATION USING BAYES DOT ALLOCATION. 

C FORMULA HERE ARE FOR BAYES WITH NO PRIOR. 

C GIVEN ANY MAX NO. OF DOTS IN MAXDOT. PROGRAM WILL RETURN 



ORIGINAL rt.Cl X 
OF POOR QUALiTf 


RM«5* 




C — — THE ESTIMATE IS IN ARRAY PSlo ON RETURN. 

C 

C UR I TEN AND EDITED BY N.Y. CHU ON 5-10-79. 

C 

C 

INTEGER NDARY(200).MM(51)»NN<51)»NX(Sl).LL(51) 

REAL UU.PSG<1) 

BYTE LABEL 

C COMMON BLOCK IS A FLAG FOR PRINT ENARLE/D ' SABLE 

COMMON /PRTFLG/JFLAG 
DATA NSTART/0/. MINSZ/5/ 

DATA IN I POT/3/ 

C 

IF'MAXDOT.GT.200) GOTO 911 
C 

C GET PIXEL COUNT FOR EACH CLUSTER 

IF(NSTART.NE.O) GOTO 121 
NSTART-1 

CALL. CLMPCS(NPIXEL.M.MM.LL) 

IF(M.GE.Sl) GOTO 901 
121 CONTINUE 
C 

C : BEGIN BAYES DOT ALLOCATION 

C 

IF(JFLAG.EQ.l) WRITE(6. 141) 

141 FORMAT (/v ■ 2SX, ' DOT FILE 1 

* /'DOT CLUSTER INFORMATION POSITION GROUND TRUTH 

* . ' SEGMENT PROPORTIONAL’ 

* /’ NUMBER NO. PIXEL DOTS S.G. LINE. COL RAW. CODE' 

* UARIANCE ESTIMATE’) 

C 

C MAKE SURE EACH CLUSTER HAS A MINIMUM NO. OF PIXELS 

K1*0 

NDS-0 

DO 1E1 K*1.M 

NX(K)«0 

NNOO-0 

IFCMMCK) .GE.MINSZ) NDS-NDS+INIDOT 
IFCMMCK) .LT.MINSZ) Kl-Kl+rtlCK) 

IFCMMCK) .LT.MINSZ) NN(K)— 1 
151 CONTINUE 
MX^M 

IFCKl. LT.MINSZ) GOTO 161 
MX-M+1 
MM(MX) «K1 
NDS-NDS+INIDOT 
161 CONTINUE 
C 

c FROM HERE TO STATEMENT 261 » LOOP TO ASSIGN 1 DOT AT A TIME 

C 

C IN ORDER TO ACHIEUE ASSIGNING 3 DOTS INITIALLY OR LESS 

C TO EACH CLUSTER. 

C THE DELTA UARIANCE COMPUTATION IS BYPASSED FOR THE FIRST 

C NDS«INIDOT*MX DOTS 

C 

KG1-0 

NDOT-0 

191 NDOT *N'- >1 

IFCHDGT GT.NDS) GOTO 2O0 
195 KG1-KG1+1 

IF(KGl.GT.MX) KG1-1 
K1-KG1 

IFCNNCKG1).GT.-1) GOTO 214 
GOTO 195 
C 
C 



C FOR EACH CLUSTER, COMPUTE DELTA 'VARIANCE AND CHOOSE THE LARGEST 

see amax-0. 

DO 211 K-l.MX 
J-NNCK) 

IFCJ.LE.l) GOTO 211 

MAR- < FLOAT(MM(K) IXNPIXELXC J-l. )AJ/CJ+1. ) )**2 

* *NXCK)#C J-NXCK) ) 

C HR I TEC 6, 99121 )K»MMCK) , J.NXCK) »UAR 

99121 FORMATC’ K.MMCK), J.NXCK), UPR-’ . 14, I8.2I5.E16.6) 

IFCUAR.LT. AMAX) GOTO 211 
Kl-K 

AMAX-UAR 
211 CONTINUE 

C 

C ASSIGN A EOT TO THE CLUSTER WITH MAX DELTA UARIANCE 

IFCM4CK1 ) .GE.MMCK1) ) WRITECS. 213>ND0T,Kl.MMCKi) 

IFchtSCKl ),GE.MMCKl) ) GOTO 241 

213 FORMAT ( ’ ATTEMPT TO ASSIGN DOT-’, 13,’ WHILE NO PIXEL’, 

* ’ IS AUAILABLE. CLUSTER" ' , 13, ' TOT PIXEL-’ , 16) 

C POSITION THE POINTER TO NDARY 

214 ND-1 

DU 216 K-l.Kl 

IFCNNCK) .LE.0) GOTO 21S 

ND-ND+NNCK) 

215 CONTINUE 

C MAKE ROOM FOR THE NEW DOT IN NDARY 

Nl— 1 

DO 217 K"NDOT,ND,Nl 
217 NDARYCK+1)-NDARYCK) 

ND-ND-NNCK1) 

NNCK1)-NNCK1)+1 
C ASSIGN A DOT 

CALL GETDOTCMMCK1 ) .NNCK1 ) .NDARYCND) , IDOT) 

C WRITECS. 99141 )K1, ND. NNCK1 ), NXCK1 ), (NDARY(K).K-l, 44) 

99141 FORMAT ( ’ Kl.ND.NN.NX"’ .414. 4C/.2X.15I4) ) 

ifcki.ljlm) goto 227 

DO 225 Kl-l.M 
IFCNNCK1).GT.0) GOTO 225 
IFCIDOT.LT.MMCK1)) GOTO 227 
I DOT • I DOT-MM ( K1 ) 

225 CONTINUE 

227 CALL CLMPXY CK1 » I DOT , > NC » LL ) 

CAL L GTMPLB < NL . NC . LABEL ) 

CALL LBLITPCLABEL.LBI.NXCKD) 

PSG(NDOT)«0. 

DO 229 JX-l.MX 
1FCNNCJX) .LE.0) GOTO 229 

PSG ( NDOT ) • PSG ( NDQT ) + FLOATCNXC JX) )/mC JX) *MM( JX5/NPIXEL 
229 CONTINUE 
C 
C 

C CCffIPUTE AND STORE SEGMENT UARIANCE 

241 UU-0. 

DO 251 K-l.MX 
IFCNNCK). LE.l) GOTO 251 
P- FLOATCNXCK) )/MS(K) 

UU-UU+CFL0ATCI-T1CK))/'NPIXEL )**2 *F*C1.-P)/C NN(K)-1. ) 

251 CONTINUE 

IFCJFLAG.EQ.i) WRITECS, 2SS)ND0T.Kt. IDOT, NNCK1) ,NXCK1) 

* .NL.NC. LABEL. LBl.UU.PSGCNDOT) 

255 F0RI1ATC2X. 13, 3X. 13. IX. 16. IX. 10, 2X. 13, 2X.2I4, 2X.2I4. 2F11.6) 
1FCNDOT.LT. MAXDOT) GOTO 191 
C 

DO 261 K-l.MX 
IFCNNCK). LE.-l) NN(K)-0 
261 CONTINUE 



END OF BAYES DOT ALLOCATION 


C 

C 

c 


C- WRITE CLUSTER INFORMATION 

IF(JFLAG.EO.l) WP I TE ( 6 . 282 ; MAXDOT . NP I XEL . MX 
282 FORMATS. 5X, ' TOTAL NO. OD DOTS AUAILABLE -M6. 

* X.SX, * TOTAL NO. OF PIXEL IN THE CLUSTER MAP -M6. 

* /.SX. * TOTAL NO. OF CLUSTERS -MS. 

* -'.5X, * CLUSTER CLUSTER NO. OF NO. OF DOTS LABELLED* 

* S'SX , ’ NO. CODE PIXELS ASSIGNED SMALL GRAIN' ) 
DO 283 K-l.MX 

IF( JFLAG.EQ. 1 > WRITE(6.284)K.LL(K) ,MM(K> .NN(K) .NX(JO 
203 CONTINUE 

284 F0RMAT(8X, I2.5X, I4.3X, I6.4X. 13. 10X. 13) 

IF(MX.EQ.M) GOrO 288 

Kx-e 

DO 285 K-l.M 

IF(NN(K) .NE.0) GOTO 28S 

KX-KX+l 

NDARY(KX)-K 

285 CONTINUE 

IF(JFLAG.K). 1) I JR I TECS. 286 )MX. (NDARYOO.K-l.KX) 

28S FORMATS* THOSE SMALLER CLUSTERS GROUPED (AS CLUSTER* • 13# 

* *) ARE • * . 3(/M6(I3. ’.*)) ) 

288 CONTINUE 

C 

C 

C 

c 

GOTO 890 
C 
C 

901 WRITE(6.902)M 

90S FORMAT ( * YOU HAUE* . 13. * TOO MANY CLUSTERS (MAX* SC ) * ) 

GOTO 99C 

911 WRITE(6. 912JND0T 

912 FORMAT ( ’ YOU SPECIFIED* . M. ’ TOO MANY DOTS (MAX* 200 j ’ ) 

C 

C 

C ENTRY TO DEFINE INTI AL NO. OF DOTS TO BE ASSIGNED 

ENTRY INIT2(NNIT) 

INIDOT-NNIT 
990 RETURN 
BSD 
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9. A87t BAYESIAN DOT ALLOCATION (QUADRATIC PRIOR) 


9.1 DESCRIPTION 

Implemented In thle piece of software is a scheme in which dot 
allocation is sequentially Bayesian (quadratic prior) and label- 
ing is direct from ground truth. This scheme uses an algorithm 

similar to that used by A84 and A85 but has different formulas 
2 

for do and segment variance. These two parameters are computed 
with the help of function THl. 


In the Bayesian dot allocation scheme, dots are allocated sequen- 
tially to those clusters that may reduce the variance of the seg- 
ment proportion estimate the most. In order to do this, we must 

2 

first compute the expected change in variance, do , for each 
cluster. The expected change in variance for cluster i is 
defined as 


2 /** i\^ 

- (r) 


- ) ®n i ,x 1 *n i +l,x i +l <1 " ^n^+l.Xj+l* 


n i' x i 


- 1 ) 




(l - e_ „ )$_ „ a - e_ „ ) 

n^,x^ n^+i.x^ n^+i,x^ 


where N^ 
N 


- number of pixels in cluster i 
■ total number of pixels in the entire scene 


:m nia : ms s « : 1 1 s inn sift : li h it : me t M-t ii 


(9-1) 


(9-2) 


a - 6.31924, b - -8.19799, c - 2.99258 

n^ ■ number of dote previously allocated to cluster i 

X| - number of dots previously allocated to cluster i 

which are labeled as small grain 

(notice x A < n^) 


It should be noted that 0’ „ 

n# x 

THl. 


la computed via a function called 


2 

Next, a dot is allocated to the cluster whose Ao^ is the largest. 
Then, for this chosen cluster, n^ is updated to ^ + 1. That 
dot's label is read from the ground truth file. If the label is 
a small grain, x^ is updated to x^ + 1. After this, the segment 
variance, o*, is computed as 



where m ■ total number of clusters 
2 

This o is printed and the above allocation scheme is repeated 
until the number of dots allocated reaches a value specified by 
the user. 


The proportion estimate is computed as 

m N a a 

P ■ ” -L e 

89 frl N n i' x i 


(9-4) 


This scheme, like other allocation and labeling schemes, may be 
repeated with different pseudorandom sequences. The main program 
resets the random number generator to a different starting point 
for each repetition. 

All computation of statistics such as mean and variance is per- 
formed in the main program. A grand summary is printed. 



9.2 LINKAGES 


The PDP 11/45 Fortran library (for functions and subroutines such 
as SQRT) , the PDP 11/45 system library (for routines such as the 
pseudorandom number generator) , and the utility package (described 
in the section 14) are required. 


Main program 

or subprogram Subroutine or function required 


MAIN 

SAEB3 

GETDOT 

THl 


SAEB3, GETDO. - CLMPCS, GTMPLB , LBLITP 
GETDOT, CLMPCS, GTMPLB, LBLITP, THl 
RAN 
None 


9.3 INTERFACE 

Interface with other routines is through the common block PRTFLG 
(2 bytes), which is used to control the optional printing of dot 
files and other information. 


9.4 INPUTS 

Fortran- formatted input of the following parameters is needed: 


Card 

or 

line 


Parameters 

Format 


Default 

input 

Default 

value 

1 

Cluster map file name 

A13 

None 


2 

a. 

Number of repetitions 

13 

0 

or blank 

1 


b. 

Starting point of first 
pseudorandom sequence 

15 

0 

or blank 

10 


c. 

Number of repetition 
printings 

13 

0 

or blank 

5 




3 a. Maximum number of dots 13 None 

that can be allocated 

b. Number of dots initially 13 None 

assigned to each cluster 

4 Number of status* messages on 13 0 or blank NO MSG 

terminal 

Two disk files are needed, one containing the stripped cluster 
map, the other the stripped ground truth map. These files should 
have been created by conversion programs A81 and A86. 


9.5 OUTPUTS 

In addition to the normal echo of input specifications and status 
messages on the terminal, the following will be output on a line 
printer. 


Report 


Contents 


Individual repetition 1. Processor header 


2. Ground truth small-grain proportion 

3. A dot file showing dots chosen, their 
x-y position, their ground truth labels, 
segment variance, and proportion 
estimate 

4. A table showing number of clusters, 
cluster name codes, cluster sizes, and 
number of dots assigned 


Grand summary 1. Processor header 

2. A table showing bias, MSE, reduction 
in MSE, average, variance, and vari- 
ance reduction as functions of number 
of dots 



9.6 BRIEF ALGORITHM: FOR BAYESIAN DOT ALLOCATION 

(QUADRATIC PRIOR) 


1. Type messages at terminal and read from it run specifications. 

2. Compute ground truth small-grain proportion. 

3. Repeat a through g until all repetitions are finished: 

a. Set starting point of the pseudorandom number generator. 

b. Clear a dot counter. 

c. Increment the dot counter. 

d. If the dot counter indicates that the current dot should 

be allocated as an initial dot, then go to the next step. 

2 

Otherwise compute A o. for every cluster and choose the 

1 2 

cluster with maximum Aa^ . 

e. Allocate a dot to the chosen cluster. 

f. Pick up the small-grain labels, compute segment variance 
and proportion estimate, and print a line of the dot file. 

g. If the dot count does not exceed specified value, go to c. 

4. Compute bias, MSE, reduction in MSE, average, variance, and 
variance reduction, and print a grand summary. 
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1 


Call CLNPCS to initial 
cluster mp road 


ID 



Compute Ao' 
Select the 
with max Ac 


Select the fluster 


( 


Assign and label a dot 
Compute segment variance 
and proportion estimate 


\ 


Print dot file 




> 



Figure 9-1.— Flow chart for Bayesian dot allocation 

(quadratic prior). 
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Figure 9-1.— Concluded. 




9.7 LISTING 


c 

c 

C PROGRAM A13 (A87.TSK) < PROPORTION ESTIMATE OF SMALL GRAIN 

c US I M3 BAYES (QUADRATIC PRIOR) DOT ALLOCATION 

C LABELLING IS DIRECT BY GROUND TRUTH. 

C 

BYTE NAME(15).NGTC13) 

REAL PSG(200).PM(200),PB(200) 

C COMMON BLOCK IS A FLAG FOR PRINT ENABLE/DISABLE 


COMMON /PRTFLG/JFLAG 

DATA NGT(5)/’ . ’/,NGT(6)/’S*/,NGT(?)/’T’/,NGT(8)/’P’/ 

C 

C 

C READING INPUT FROM TERMINAL 

C 

WRITEO, 121) 

121 FORMATC ’ PROGRAM A13 (A87.TSK).’ 

* /’ PROPORTIONAL ESTIMATION OF SMALL GRAIN' 

* /’ DOT ALLOCATION IS BAYES (QUADRATIC PRIOR).’ 

* /’ LABELLING IS DIRECT BY GROUND TRUTH.’ 

* /’ INPUT CLUSTER MAP FILENAME’/’ AAAAAAAAAAAAA’ ) 

READ ( 7 » 122) (NAME(K) ,K=1, 13) 

122 FORMAT (13A1) 

C 

WRITEO, 131 ) 

131 FORMAT ( ’ AN ASSIGWENT OF PSEUDORANDOM DOTS IS A REPETITION RUN’ 

* /’ HOW MANY REPETITION RUNS? SEED FOR FIRST PSEUDO SEQ?’ 

* /’ HOW MANY PRINTINGS OF INDIUIDIIAI PUNA?’ 

* /’ III 1IIII III’) 

READ ( 7 , 132 ) JU . JSEED . JPAGE 

132 F0RMATU3, IX. 15, IX, 13) 

IF( JU, LE. 0) JU-i 

IF( JSEED.LE.0) JSEED” 10 
IF( JPAGE. LE.0) JPAGE-5 
JSK I P- ( JU-1 ) /JPAGE+1 
C 

WRITEO, 141 ) 

141 FORMATC SPECIFY riAXIMUM NO, OF DOTS IN THE REPETITIONS RUNS,’/ 

* ’ AND NO. OF DOTS I NT I ALLY ASSIGNED TO EACH CLUSTER.’/ 

X ’ III IIP) 

READC7, 142)MAXD0T.NNIX 

142 FORMAT ( 13, IX, 13) 

CALL INIT3(NNIX) 

C 

WRITEO, 151 ) 

151 FORMATC STATUS MESSAGES ON TERMINAL. HOW MANY?’/’ Ill’) 

READ (7,152) NMES 

152 FORMAT (1 3) 

C 

WR I TE ( 8 , 18 1 ) ( NAME (K) , K-l, 13) , JU, JSEED , MAXBOT 
181 FORMAT (//’ CLUSTER FILE GIUEN - ’,13A1. 

* /’ NO. OF REPETITION RUNS FOR EACH TOTDOT ■’,13, 

* /’ THE FIRST REPETITION RUN STARTS WITH SEED » ’,16, 

* /’ MAXIMUM NO. OF DOTS IN THE REPETITION RUNS”’, 16) 

C 

C COMPUTE BAYES (QUADRATIC PRIOR) ESTIMATE FOR GROUND TRUTH 

NLINE”117 
NCOL-196 
DO 215 K-1,4 
215 NGT(K)”NAME(K) 

OPEN (UNIT-2, NAME-NGT . TYPE- ’ OLD ’ . READONLY , FORM- ’ UNFORMATTED ’ , 

* ACCESS- ’ DIRECT’ ) 

DO 231 L-l.NLINE 

DO 221 K-l.NCOL 
CALL GTMPLB (L, K, LABEL ) 

CALL LBL I TP ( LABEL >LB1> IP) 

221 CONTINUE 

231 CONTINUE 


7 ° 


reproducibility of the 

ORIGINAL PAGE IS POOR 


P«FLOAT< I P > /N . I NE/NCOL 
C 

OPEN(UNIT"l, NAME-NAME, TYPE-’ OLD’ , FORM- * UNFORMATTED ’ , 

* READONLY , ACCESS- * D I RECT* ) 

CALL CLMPLC ( I NE > NCQ_ ) 

C 

C START ESTIMATION FOR AN TGTDOT AND REPET IT I CMS 

DO 3X1 K-l. MAXDOT 
PB(K)-0. 

311 Pfl(K) =0. 

JMES-0 
JS-JSEED 
DO 361 J-l.JU 
JFLAG-0 

IF<MOD< J-l, JSK1P) .EQ.0) JFLAG-l 
IFCJFLAG.EQ.l) WRITE(6,321)NAME,F 

321 FORMAT C ’ 1’ , 10X, ’ACCURACY ACCESSEMENT SOFTWARE (S-l 1-79) ’ ,^,2X» 

X ' PROGRAM A13 CA87.TSK)' ESTIMATION OF SMALL GRAINS PROPORTION’ 
X ' > 1SX . ’ DOT ALLOCAT I ON « BAYES ( QUADRAT I C PRIOR),’ 

X 1SX, ’DOT LABELLING 'DIRECT BY GROUND TRUTH.’ 

X / ', 10X, ’ INPUT CLUSTER MAP IS FILE ’ , 1SA1 . 

X ✓, 10X,’ PROPORTIONAL ESTIMATE FOR GROUND TRUTH IS’.FS.S) 

I F < JFLAG . EQ . 1 ) WRITERS, 322)J.JSEED 

322 FORMAT ( / , 2X , ’ REPETITION RUN-’. 13. 

* ’ RANDOM DOT SEQUENCE STARTS WITH’ , 18, ’ ---’ ) 

CALL RANSTC JSEED) 

CALL SAEB3 ( MAXDOT , PSG ) 

JMES-JMES+l 

I F ( JMES . LE . NMES ) WR I TE ( 8 . 323 ) MAXDOT , J , JSEED 

323 FORMAT ( ’ MAXDOT- ’, 13, ’ REPETITION RUN-’ , 13, ’ SEED-’, 16) 

JSEED* JSEED+ 150 

C COMPUTE BIAS AND M.S.E. 

DO 341 F-l, MAXDOT 
TEMP - PSG (K)-P 
PB ( K) =PB ( K ) +TEMP 
341 FM<K)-PMlK)+TEMP**2 
361 CONTINUE 
C 

C PRINT GRAND SUMMARY FOR THIS JOB 

WRITEC6, 321 )NAME, P 

WR I TE C 6 . 371 ) MAXDOT , JU , JS 

371 FORMAT (/, 14X, 5 GRAND SUMMARY OF THIS JOB ’, 

* /,SX,’ MAX. NO. OF DOTS IN EACH REPETITION PUN-’ , 13, 

X /, 14X, ’ NO. OF REPETITION RUNS 13, 

* 14X, ’ RANDOM DOTS START WITH SEED*’, 16, 

* vV,lX,’ DOT BIAS M.S.E. REDUCTION’, 

X ’ AUERAGE UARIANCE UAR REDUCTION’ J 

DO 381 E-i. MAXDOT 
PBCK)-PB(K)/JU 
PMOO»PM(K)/JU 
RR-PMCO/C PX<1.-P)/K ) 

AUERG-PBGO+P 

IF'JU.EQ.l) UAR* PM ( K ) -PB ( IO XX2 
IFCJU.GT.l) UAR- (PMiK)-PB(K)XX2)*JU/'(JU--l) 

PEDUAR=UAR/( P*(l.-P)/K ) 

381 WR I TE C 6 , 382 ) K , PB ( K ) , PM ( K ) , RR , AUERG , UAR . REDUAR 

382 FORMAT (2X. 13, 2X.F9.5, 4(F10.6), 3X.F10.6) 

C 

C 

WRITEC6, 401 ) 

401 FORMAT ( ’ 1 END OF THIS JOB ’) 

C 

STOF 

END 

C 

c 

c 
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SUBROUTINE SAEB3 ( MAXDOT . PSG ) 


C STRATIFIED AREAL ESTIMATION USING BAYES DOT ALLOCATION. 

C FORMULAS HERE ARE FOR BAYES WITH QUADRATIC PRIOR. 

C GIUEN ANY MAX NO. OF DOTS IN MAXDOT. PROGRAM WILL RETURN 

C THE PROPORTIONAL ESTIMATE IS IN ARRAY PSG. 


C 

C WRITEN AND EDITED BY N.Y. CHU ON 5-11-79. 

C 

c 

INTEGER NDARYC200) »MM(51) .NNC51) »NX(51) .LLC51) 
REAL UV.PSG(l) 

BYTE LABEL 

C COMMON BLOCK IS A FLAG FOR PRINT ENABLE/DISABLE 

COMMON /PRTFLG/JFLAG 
DATA NSTART/0/. MINSZ/5/ 

DATA INI DOT/3/ 

C 

IF r MAXDOT. GT. 200) GOTO 911 
C 

C GET PIXEL COUNT FOR EACH CLUSTER 

1FCNSTART.NE.0) GOTO 121 
NSTART-1 

CALL CLMPCS ( NF I XEL • M > MM • LL ) 

IF(M.GE. 51 ) GOTO 901 
121 CONTINUE 
C 


C BEGIN BAYES DOT ALLOCATION 

C 

IFCJFLAG.EQ. 1) WRITEC6, 141) 

141 FORMAT ( // • 25X. ’ DOT FILE ’ 

* /’DOT CLUSTER INFORMATION POSITION GROUND TRUTH’ 

* SEGMENT PROPORTIONAL’ 

* /’ NUMBER NO. PIXEL DOTS S.G. LINE. COL RAW. CODE’ 

* , ’ VARIANCE ESTIMATE’ ) 

C 

C MAKE SURE EACH CLUSTER HAS A MINIMUM NO. OF PIXELS 

Ki-0 


NDS-0 

DO 1S1 K-l.M 

NXCK1-0 

NN(K)-0 

IF(MMCK).GE.MINSZ) NDS-NDS+INIDOT 
IFCMMCK) . LT.MINSZ) K1-K1+MMOO 
IF(MM(K). LT.MINSZ) NN(K)— 1 
151 CONTINUE 
MX»M 

IFCKl. LT.MINSZ) GOTO 161 
MX-M+i 
MM(MX) «K1 
NDS-NDS+INIDOT 
161 CONTINUE 
C 

C FROM HERE TO STATEMENT 261 « LOOP TO ASSIGN 1 DOT AT A TIME 

C 

C IN ORDER TO ACHIEVE ASSIGNING 3 DOTS INITIALLY TO EACH CLUSTER 

C THE DELTA VARIANCE COMPUTATION IS BYPASSED FOR THE FIRST 

C NDS* IN1DQTKMX DOTS 

C 

KG1-0 

NDOT-0 

191 NDOT-NBOT+1 

IF(NDOT.GT.NDS) GOTO 200 
195 KG1-KG1+1 

IFCKG1.GT.MX) KG1-1 
K1-KG1 

IFCNN(KGl).GT.-l) GOTO 214 
GOTO 195 
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c 

c 

C FOE EACH CLUSTER. COMPUTE DELTA UARIANCE AND CHOOSE THE LARGEST 

200 AMAX-0. 

DO 211 K-l.MX 
J-NN(K) 

JX-NX(K) 

IF(J.LE.l) GOTO 211 
C 

UAR- ( FLQAT(MMCK) )/NPIXEL )**2 

* *< THii'J. JX)#(1.-TH1(J. JX))/(J-1. ) 

* -THKJ. JX)*TH1< J+l. JX+-1)*(1.-TH1(J+1. JX+1) )AJ 

* -( 1 . -TH1C J. JX) )XTH1( J+l. JX)*( 1, -TH1 ( J+l. JX) )/J ) 

C 

C URITE(6,99121)K,MM(K).J,NX(K).UAR 

99121 FORMAT C * K. MMCK) , J. NXCK) .UAR- ’ , 14, 18. 215, E16. 6) 

IFCUAR.LT.AMAX) GOTO 211 

ri*K 

AMAX-UAE 
211 CONTINUE 
C 

C ASSIGN A DCIT TO THE CLUSTER WITH MAX DELTA UARIANCE 

IF(NN(K1).GE.MM(K1)) WRITEC6. 213)N1)0T, K1 .HM(Kl) 

IFCNNCK1 ) . 6E. MM(K1 ) ) GOTO £41 

213 FORMAT f ’ ATTEMPT TO ASSIGN DOT*’ , 13. ’ WHILE NO PIXELS 

* ’ IS AUAILABLE. CLUSTER*’ , 13. ’ TOT PIXEL*’ . 16) 

C POSITION THE POINTER TO NDARY 

214 ND*1 

DO £15 K-l.Kl 
IFCMN(K) . LE. 0) GOTO 215 
ND-ND+NNOO 
216 CONTINUE 


C MAKE ROOM FOR THE NEW DOT IN NDARY 

N1-- -1 

DO 217 K-NDOT.ND.N1 
217 NDARYCK+1) -NDARY (K) 

ND-ND-NNiKl) 

NN(K1 ) *NN(K1 )+l 
C ASSIGN A DOT 

CALL GETDOTCMM(Kl) .NN(Kl) , NDARYIND) . IDQT) 

C HR I TE < 6 , 99 1 4 1 ) K 1 , ND . NN < K 1 ) . NX ( K 1 ) , t NDARY ( K ) , K - 1 , 44 ) 


99141 FORMAT < ’ K1.ND.NN.NX-’ .414. 4C/.2X.1514) ) 

IFCK1.LE.M) GOTO 227 
DO 225 Kl-l.M 
IF(NN(K1).GT.0) GOTO 225. 

IFCIDOT.LT.MM(Kl)) GOTO 227 
IDOT* IDOT-MM(Kl ) 

225 CONTINUE 

227 CALL CLMPXYCK1 . IDOT. NL . NC, LL) 

CALL GTMFLBCNL.NC, LABEL) 

CALL LBLITPCLABEL.LBl.NXOa) ) 

PSG(NDOT) *0. 

DO £29 JX- 1 . *4X 
I F ( NN ( JX ) . L. E . 0 ) GOTO 229 

PSG ( NDOT ) -PSG C NDOT ) + FLOAT f MI1 C JX ) ) /NP I XELXTH1 C NN < JX ) . NX ( JX ^ ) 
229 CONTINUE 
C 

c 

C COMPUTE SEGMENT UARIANCE 

£41 UU-0. 

DO 251 K-l.MX 
IF(NN(K).LE. 1) GOTO 251 

UU-UU+ ( FLOAT ( MM NO ) /NP I XEl. ) **2 XTH1(NN<K) .NX(K) ) 

* YC 1 . -TH1 <NN(K) . NXCK > ) ) / <!NN(K )-l . ) 

251 CONTINUE 

IFCJFLAG.EQ. 1) HRITEC6, 255) NDOT. K1 . IDOT.NN(Kl) .NXCK1 ) 

T , NL . NC . LABE! . . LB 1 . UU . PSG < NDOT ) 

255 FORMAT C2X, 13. 3X, 13, IX. 16. IX. 13. 2X. 13. 2X.2I4, 2X.2I4. 2F11.6) 
I F ( NDOT . LT . I4AXDOT ) GOTO 191 
C 

DO 261 K-l.MX 
IF(NNCK).LE.-l) NN(K)«0 





261 CONTINUE 

C 

C 

C END OF BAYES COT ALL OCATION 

C 

c 

C WRITE CLUSTER INFORMATION 

IF(JFLAG.EQ.l) WRITEC6,2aT)MAXD0T,NPIXFL.MX 
282 FORMAT ( x/', 5X» * TOTAL NO. OD DOTS AVAILABLE , 16, 

* x.SX,’ TOTAL NO. OF PIXEL IN THE CLUSTER MAP -’.16, 

X x.5X. ’ TOTAL NO. OF CLUSTERS , 16. 

* x,5X, ' CLUSTER CLUSTER NO. OF NO. OF DOTS LABELLED' 

* /.SX, ' NO. CODE PIXELS ASSIGNED SMALL GRAIN’ ) 


DO 283 K-l.MX 

IF(JFLAG.EQ.l) WRITEC6.c84)K,LL(K) ,MM(K) ,NN(K) ,NXCK) 

283 CONTINUE 

284 F0RW)T(8X, I2.5X, I4.3X, I6.4X, 13, 10X, 13) 

IF(MX.EQ.M) GOTO 288 

KX-0 

DO 28S K-l.M 

IF(NNCK) .NE.W) GOTO 28S 

KX-KX+1 

NDARY(KX)-K 

285 CONTINUE 

IF(JFLAG.EQ.l) WRITEC6,286)MX, (NDARYOO .K-l.KX) 

286 FORMAT ex' THOSE SMALLER CLUSTERS GROUPED (AS CLUSTER’ , 13, 

* ’) ARE 1 ’ , 2<.', 16(13, ’.’)) ) 

288 CONTINUE 

C 

C 

C 

C 

GOTO 990 

C 

c 

901 WRITEC6, 902)M 

’302 FORMAT ( ’ YOU HAUE’ , 13, ’ TOO MANY CLUSTERS (MAX -50) ’ ) 

GOTO 990 

911 WRITE(6.912)ND0T 

912 FORMAT C YOU SPECIFIED’ , 14, ’ TOO MANY DOTS (MAX -200) ’ ) 

C 

C 

C ENTRY TO DEFINE INTIAL NO. OF DOTS TO BE ASSIGNED 

ENTRY INIT3(NNIT) 

INIDOT-NNIT 
990 RETURN 
END 
C 
C 

FUNCTION TH1CN.NX) 


C FUNCTION TO EUALUATE CLUSTER PROPORTIONAL ESTIMATE 

C N-NO. OF DOTS ASSIGNED FOR A CLUSTER, 

C NX- NO. OF DOTS LABELLED AS SMALL GRAIN 

C BOTH N AND NX ARE SIMPLE INTEGER (2 BYTES IN PDP) 

C TH1 (REAL NO.) CONTAINS THE ESTIMATE CM RETURN. 


C 

INTEGER N, NX 
REAL TH1 
C 

DATA A/6.0Y, B/-7.8??/, C/2.9385X 
C 

X-NX 

TH1 - 

* ( A*(X+1. )#CX+2. )*CX+3. ) + B*(X+1. )#(X+2. )*(N+4. ) 
X + CXCX+l. )*(N+3. )*(N+4. ) ) ' 

* ( A*(X+1. )#CX+2. )*CN+4. ) + B*<X+1. )*(N+3. )*<N+4. ) 
X + C*(N+2. )*(N+3. )*(N+4. ) ) 

RETURN 

END 
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10. A89 1 BAYESIAN DOT ALLOCATION (MODIFIED QUADRATIC PRIOR) 


10.1 DESCRIPTION 

Implemented in this piece of software Is a scheme in which dot 

allocation is sequentially Bayesian (modified quadratic prior) 

and labeling is direct from ground truth. This scheme is similar 

to those in A84, A85, and A87, but it uses different formulas for 
2 

La and segment variance. The constant used in the function TH2 
is also modified during execution. 


In the Bayesian dot allocation scheme, dots are allocated sequen- 
tially to those clusters that may reduce the variance of the seg- 
ment proportion estimate the most. In order to do this, we must 

2 

first compute the expected change in variance, La , for each 
cluster. The expected change in variance for cluster i is 
defined as 


La 


2 

i 



V x i 


(i - e 


n< 


(n j - 1) 



A A A 

9 n i' x i 9 n i' 4 ' 1 ' x i +1 * 1 “ ®n i +l,x i +l* 

- ■ 


( 10 - 1 ) 

where N^ » number of pixels in cluster i 

N * total number of pixels in the entire scene 


n,x 


• (x ♦ 1) (x ♦ 2) (x ♦ 3) ♦ b (x ♦ 1) (x ♦ 2) (n ♦ 4) + c (x ♦ 1) (n ♦ 3) In ♦ 4) 
«ix I nix ♦ IRnT TF * b ' lx + ' 1) fn l ij ' h ♦ i c )r T V iHn * i{ hi "» if ' 


( 10 - 2 ) 


For initia. dots* a ■ 6, b - -7.877, c - 2.9345 

After initial dots* a - 6, b - 12(P #9 - 1) , c ■ 5 - 6P # ^ if 

0.211 < P„ s 0.789 
■9 

a - 6, b - 12(0.211 - 1), c - 5 - 6 x 0.211 


if 0.211 > P.„ 

■9 

a - 6, b - 12(0.789 - 1), c - 5 - 6 x 0.789 

if 0.789 < P.„ 

•9 

- number of dots previously allocated to cluster i 

■ number of dots previously allocated to cluster i 
which are labeled as small grain 

(notice x ^ s n^) 


It should be noted that __ is computed via a function called TH2 

n, x 


2 

Next, a dot is allocated to the cluster whose Ao^ is the largest. 
Then, for this chosen cluster, n^ is updated to n^ + 1. That 
dot's label is read from the ground truth file. If the label is 
a small grain, x^ is updated to x^ + 1. After this, the segment 
variance, a 2 , is computed as 



( 10 - 3 ) 


where m * total number of clusters 


2 

This a is printed and the above allocation scheme is repeated 
until the number of dots allocated reaches a value specified by 
the user. 


The proportion estimate is computed as 



This scheme , like other allocation and labeling schemes, may be 
repeated with different pseudorandom sequence*. The main program 
resets the random number generator to a different starting point 
for each repetition. 


All computation of statistics such as mean and variance is per- 
formed in the main program. A grand summary is printed. 


10.2 LINKAGES 

The PDP 11/45 Fortran library (for functions and subroutines such 
as SQRT) , the PDP 11/45 system library (for routines such as the 
pseudorandom number generator), and the utility package (described 
in section 14) are required. 


Main program 

or subprogram Subroutine or function required 


MAIN 

SAEB4 

GETDOT 

TH2 


SAEB4, GETDOT, CLMPCS, GTMPLB, LBLITP 
GETDOT, CLMPCS, GTMPLB, LBLITP, TH2 
RAN 
None 


10.3 INTERFACE 


Interface with other routines is through the common block PRTFLG 
(2 bytes), which is used to control the optional printing of dot 
files and other information. 


10.4 INPUTS 

Fortran-formatted input of the following parameters is needed: 


Card 

or 

line Parameters 

1 Cluster map file name 

2 a. Number of repetitions 


Default 

Default 

Format input 

value 

A13 None 

13 0 or blank 

1 





b. 

Starting point of first 
pseudorandom sequence 

15 

0 or 

blank 

10 


c. 

Number of repetition 
printings 

13 

0 or 

blank 

5 

3 

a. 

Maximum number of dots 
that can be allocated 

13 

None 




b. 

Number of dots initially 
assigned to each cluster 

13 

None 



4 

Number of status messages on 
terminal 

13 

0 or 

blank 

NO MSG 


Two disk files are needed, one containing the stripped cluster 
map, the other the stripped ground truth map. These files should 
have been created by conversion programs A81 and A86. 


10.5 OUTPUTS 

In addition to the normal echo of input specifications and status 
messages on the terminal, the following will be output on a line 
printer. 


Report Contents 

Individual repetition 1. Processor header 

2. Ground truth small-grain proport io.. 

3. A dot file showing dots chosen, their 
x-y position, their ground truth labels, 
segment variance, and proportion estimate 


4. A table showing number of clusters, 

cluster name codes, cluster sizes, and 
number of dots assigned 



Grand summary 1. Processor header 

2. A table showing bias, MSE, reduction 
in MSE, average, variance, and vari- 
ance reduction as functions of number 
of dots 

10.6 BRIEF ALGORITHM: FOR BAYESIAN DOT ALLOCATION 

(MODIFIED QUADRATIC PRIOR) 

1. Type messages at terminal and read from it run specifications 

2. Compute ground truth small-grain proportion. 

3. Repeat a through h until all repetitions are finished: 

a. Set starting point of the pseudorandom number generator. 

b. Clear a dot counter. 

c. Increment the dot counter. 

d. If the dot counter indicates that the current dot should 

be allocated as an initial dot, then go to the next step. 

2 

Otherwise compute A a . for every cluster and choose the 

1 2 

cluster with maximum A . 

e. Allocate a dot to the chosen cluster. 

f. Pick up the small-grain labels, compute segment variance 
and proportion estimate, and print a line of the dot file 

g. If this is the end of initial dot assignment, compute the 
modified parameters a, b, c based on the current small- 
grain estimate, P g g. Call TH2PAR to set the modified 
parameters . 

h. If the dot count does not exceed specified value, go to c 

. Compute bias, MSE, reduction in MSE, average, variance, and 
variance reduction, and print a grand summary. 




4 



















10.7 


LISTING 


C PROGRAM A1S (A89.TSK) > PROPORTION ESTIMATE OF SMALL GRAIN 

C USING BAYES (MODIFIED QUAD PRIOR) DOT ALLOCATION 

C LABELLING IS DIRECT BY GROUND TRUTH. 

C 


BYTE NAME<15),NGT<13) 

REAL P5G(200).PM(200),PB(2O0) 

C COMMON BLOCK IS A FLAG FOR PRINT ENABLE/DISABLE 

COMMON /PRTFLG/JFLAG 

DATA NGT(5>/’ . ’/.NGTC6)/’S’/.NGT(?)/ , T , ^.NGT(8)/’PV 
C 
C 

C READING INPUT FROM TERMINAL 

C 

URITEC8. 121 ) 

121 FORMAT ( ’ PROGRAM: A16 (AS9.TSK) . * 

* /’ PROPORTIONAL ESTIMATION OF SMALL. GRAIN’ 

X /’ DOT ALLOCATION IS BAYES (MODIFIED QUAD PRIOR).’ 

* /’ LABELLING IS DIRECT BY GROUND TRUTH.’ 

X /’ INPUT CLUSTER MAP FILENAME’/’ AAAAAAAHAAAAA ’ ) 

READ (7. 122) (NAHE(K) ,K-1 , 13) 

122 FORMAT ( 13A1 ) 

C 

WRITEC8, 131) 

131 FORMAT ( ’ AN ASSIGNMENT OF PSEUDORANDOM DOTS IS A REPETITION PUN’ 

* /’ HOW MANY REPETITION RUNS? SEED FOR FIRST PSEUDO SEQ?’ 

* /’ HOW MANY PRINTINGS OF INDIUIDUAL RUNS?’ 

* /’ III Mill III’) 

READ ( 7 , 132 ) JU , JSEED , JPAGE 

132 FORMAT ( 1 3 < IX, IS, IX, 13) 

IFCJU.LE.0) JU=i 

IF( JSEED. LE.0) JSEED* 10 
IFf JPAGE.LE.0) JPAGE =5 
JSK I F = C JU- 1 ) / JPAGE+ 1 
C 

WRITE<8> 141 ) 

141 FORMAT (’ SPECIFY MAXIMUM NO. OF DOTS IN THE REPETITIONS RUNS, V 

* ’ AND NO. OF DOTS I NT I ALLY ASSIGNED TO EACH CLUSTER.’/ 

X ’ II I I II ’ ) 

READ (7,142) MAXDOT , NN I X 

142 FORMAT ( 13, IX, 13) 

CALL INITSCNNIX) 

r 

WRITE(8, 151) 

151 FORMAT ( * STATUS MESSAGES ON TERMINAL. HOW MANY?’/’ Ill’) 

READ (7, 152) MIES 

152 FORMAT ( 13) 

UR I TE ( 8 , 1 8 1 ) ( NAME ( K ) , K = 1 , 1 3 ) • JU , JSEED , MAXDOT 
181 FORMATi//’ CLUSTER FILE GIUEN = ’,13A1, 

X /’ NO. OF REPETITION RUNS FOR EACH TOTDQT *’ , 13, 

X /’ THE FIRST REPETITION RUN STARTS WITH SEED - ’ , 16, 

X /’ MAXIMUM NO. OF DOTS IN THE REPETITION RUNS*’, 16) 

C 

C COMPUTE BAYES (MODIFIED QUAD PRIOR) ESTIMATE FOR GROUND TRUTH 

INLINE* 117 
NC0L=196 
DO 215 1=1.4 
215 HGT(K) ■•TIAMEOO 

OPEN ( UN IT-2, NAME*NGT , TYPE- ’OLD’, READONLY , FORM* ’ UNFORMATTED ’ , 

X ACCESS* ’ DIRECT’ ) 

DO 231 L=1 .INLINE 
DO 221 K" 1 , NCOL 
CALL GTMPLB ( L , Y . , LABEL ) 

CALL LBL I TP ( LABEL ,LB1, IP) 

221 CONTINUE 

231 CONTINUE 

P * FLOAT ( I P ) /NL 1 1 IE/ NCOL 




rn 


0PJ^(UNIT»1. NAME-NAME, TYPE-SOLD*. FORM- 
* READONLY. ACCESS-’ DIRECT’ ) 

CALL CLMPLC ( NL I NE . bOX. ) 


’UNFORMATTED’ 


C 


311 


321 


- START EST I MAT I ON FOR AN TOT DOT AND REPET I T I CHS 
DO 311 K-l, MAXDOT 
PBCK) -0. 

PMOO-0. 

JMES-0 
JS-JSEED 
DO 361 J-l.JU 
JFLAG-0 

I F ( MOD ( J~1 > JSK l P ) . EQ . 0 ) JKLAG-1 
IFCJFLAG.EQ.l) URITEC6, 321 )NAME. P 
FORMAT t ’ 1’ , 10X. ’ACCURACY ACCESSEMENT SOFTWARE ( 5-25-79 ) */,2X 




* ’PROGRAM A16 (A89.TSK)i ESTIMATION OF SMALL GRAINS PROPORTION’ 

* /> 15X. ’ DOT ALLOCATimiBAYES (MODIFIED QUAD PRIOR).’ 

* /. 15X, ’ DOT LABELLING "DIRECT BY GROUND TRUTH.’ 

X //, 10X. ’ INPUT CLUSTER MAP IS FILE ’.1SA1, 

X /, 10X, ’ PROPORTIONAL ESTIMATE FOR GROUND TRUTH IS’.FS.S) 

I F ( JFLAG . EQ . 1 ) WR I TE ( 6 , 322 ) J , JSEED 

322 FORMAT C / . 2X > ’ REPETITION RUN-’ . 13. 

* ’ RANDOM DOT SEQUENCE STARTS WITH’ , 13, ’ ’) 

CALL RANSTC JSEED) 

CALL SAEB4 ( MAXDOT . PSG ) 

JMES-JMES+1 

I F ( JMES . LE . NMES ) WR I TE ( 8 , 323 ) MAXDOT , J , JSEED 

323 FORMAT ( ’ MAXDOT- ’. 13. ’ REPETITION RUN-’ , 13. ’ SEED-’, 16) 

JSEED- JSEED+150 

C COMPUTE BIAS AND M.S.E. 

DO 341 K-l, MAXDOT 
TEMP - PSG ( K ) -P 
PBCK) =PB(K)+TEMP 
341 PM C K ) * PM ( K ) +TFMPYY2 

361 CONTINUE 


331 

382 

C 

c 

401 

C 


C 

r 


— PRINT GRAND SUMMARY FOR THIS JOB 
UR I TE ( 6 , 321 ) NAME , P 
NR I TE ( 6 , 37 1 ) MAXDOT . JU , JS 

FORMAT C/.14X, ’ GRAND SUMMARY OF THIS JOB 

X X.SX,’ MAX. NO. OF DOTS IN EACH REPETITION PUN- 

X /, 14X, ’ NO. OF REPETITION RUNS -*,I3, 

X / , 14X, ’ RANDOM DOTS START WITH SEED- ’ , 16, 

X fv.lX,’ DOT BIAS M.S.E. REDUCTION’, 

X ’ AUERAGE UAPIANCE UAR REDUCTION’ ) 

DO 331 K-i, MAXDOT 
PBCK)»RB(K)/JU 
PMC K) -PMC K' 'JU 
RR-PMCK)XC R+:(l. -PVK .> 

AUERG « PB C K ) +P ' 

IFCJU.EQ.l) UAR- PMCIO -PB(K)**2 
IFCJU.GT.l) UAR- CPMCK)- PBCK)**2)*JU/( JU-1) 
REDUAK-UAR/C PX(1.-P)/K ) 

NR I TE C 6 , 382 ) K , PB ( K ) , PM C K ) . RR , AUERG , UAR , PEDUAR 
FORMAT C2X, 13. 2X.F9.5, 4CF10.6), 3X.F10.6) 


13, 


WRITE (8, 401) 
FORMAT C ’ 1 -- 

STOP 

END 


- END OF THIS JOB 
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SUBROUTINE SAEB4 ( MAXDOT . PSG ) 


C STRATIFIED AREAL ESTIMATION USING BAYES DOT ALLOCATION. 

C FORMULAS HERE ARE FOR BAYES WITH MODIFIED QUAD PRIOR. 

C GIUEN ANY MAX NO. OF DOTS IN NAXDOT. PROGRAM WILL RETURN 

C THE PROPORTIONAL ESTIMATE IS IN ARRAY PSG. 

C 

C WRITEN AND EDITED BY N.Y. CHU ON S-2S-79. 

C 


C 

INTEGER NDARY<200),MM(5i),NNCSl),NX(51),LL<Sl) 
REAL UU.PSG(l) 

BYTE LABEL 

C COMMON BLOCK IS A FLAG FOR PRINT ENABLE/DISABLE 

COMMON xPRTFLG/JFLAG 
DATA NSTARTX0/, MINSZ/5X 
DATA INI DOT/3/ 

C 

I F ( MAXDOT . GT . 200 ) GOTO 911 
C 

C GET PIXEL COUNT FOR EACH CLUSTER 

IFCNSTART.NE.0) GOTO 121 
NSTART-1 

CALL CLMPCS < MP I XEL . M . MM . LL ) 

IF(M.GE.Sl) GOTO 901 
121 CONTINUE 
C 

C BEGIN BAYES DOT ALLOCATION 

C 

IF(JFLAG.EQ.l) WRITEC6.141) 

141 FORMAT (xx, 25X. ’ DOT FILE * 


* /’ DOT CLUSTER INFORMATION POSITION GROUND TRUTH’ 

* . ’ SEGMENT PROPORTIONAL’ 

* x’ NUMBER NO. PIXEL DOTS S.G. LINE, COL RAW, CODE* 

* . ’ UARIANCE ESTIMATE’ ) 

C 

C MAKE SURE EACH CLUSTER HAS A MINIMUM NO. OF PIXELS 


Kl-0 

NDS-0 

DO 151 K-l.M 

NXCK1-0 

NN(K)-0 

IF(MMCK).GE.MINSZ) NDS-NDS+INIDOT 
IFCMM(K) . LT.MINSZ) Kl-Kl+MMCK) 

IF(MMCK).LT.MINSZ) NNCK)— 1 
151 CONTINUE 
MX=M 

IFCK1.LT. MINSZ) GOTO 161 
MX-M+1 
MMCMXl-Kl 
NDS-NDS+INIDOT 
161 CONTINUE 
C 

C FROM HERE TO STATEMENT 261 : LOOP TO ASSIGN 1 DOT AT A TIME 

C 

C IN ORDER TO ACHIEUE ASSIGNING 3 DOTS INITIALLY TO EACH CLUSTER 

I- Tup ptri TO Mr»pTpKH~TT rniid |TOj [ Q,,| re 0ypQCC;£ri FHR THF FTOQT 

C NDS -INI DOT #MX DOTS 

C 

KG1-0 

NDOT-0 

191 NDOT-NDOT+1 

IFCNDOT.GT.NDS) GOTO 200 
195 KG1-KG1+1 

IF(KGl.GT.MX) KG1-1 
K1-KG1 

IFCNN(KGl).GT.-l) GOTO 214 
GOTO 195 

C 

C 



C FOR EACH CLUSTER, COMPUTE DELTA VARIANCE AND CHOOSE THE LARGEST 

200 AMAX-0. 

DU 211 K-l.MX 
J-NN(K) 

JX-NX(K) 

IFCJ.LE.l) GOTO 211 
C 

UAR- ( FLOAT (MM CIOl/'NPIXEL )*#2 

* *C TH2(J,JX)*(1.-TH2(J,JX))/'(J-1. ) 

* -TH2 C J , JX ) *TH2 ( J+ 1 , JX+ 1 ) * ( 1 . -TH2 ( J+l , JX+ 1 ) ) ✓ J 

* - (1.-TH2CJ, JX) >*TH2(J+1. JX)*(1.-TH2(J+1, JX) )/J ) 

C 

C WRITEC6, 99121 )K,MM(K),J.NXCK),UAR 

99121 FORMAT (' K.MM(K) . J.NXCK) ,UAR* ’ , 14, IS,2IS,E16.6) 

IFCUAR.LT. AhHX) GOTO 211 
Kl-K 

AMAX-UAR 
211 CONTINUE 
C 

C ASSIGN A DOT TO THE CLUSTER WITH MAX DELTA UARIANCE 

IFCNN(Kl).GE.MMCKl)) WRITEC6,213)ND0T,K1,MM(K1) 

IFCNN(Kl) .GE.ftKKl) ) GOTO 241 

213 FORMAT ( ’ ATTEMPT TO ASSIGN DOT*’ . 13, ’ WHILE NO P»XEL’ . 

* ’ IS AUAILABLE. CLUSTER-’ , 13, ’ TOT PIXEL* 16) 

C POSITION THE POINTER TO NDARY 

214 ND*1 

DO 215 K-l.Kl 

I F < NN < IO . LE . 0 ) GOTO 215 

ND-ND+NNCK) 

215 CONTINUE 

C MAKE ROOM FOR THE NEW DOT IN NDARY 

Nl»-1 

DO 217 K*NDOT.ND,Nl 
217 NDARY (MU) « NDARY ( K ) 

ND*ND-NNCK1) 

NNCK1 ) =NN(K1 )+l 
C ASSIGN A DOT 

CALL GETDOTCMMCK1 ) . NN(K1 ) • NDARY (ND) , IDOT) 

C WR ITEC6, 99141 JK1.ND. NN(Kl) .NXCKt) , (NDARYCK) ,K“t .44) 

99141 FORMATC ’ K1 , ND, NN.HX- ’ , 414. 4</,2X.15I4) ) 

IFCK1.LE.M) GOTO 227 

DO 225 K1=1.M 

IF(NN(K1 ) .GT.0) GOTO 225 

IF( IDOT.L.T.MM(Kl) ) GOTO 227 

IDOT* IDOT-MMCK1 ) 

225 CONTINUE 

227 CAL L CLNPXY C K 1 . I DOT , NL . NC . LL ) 

CAl .L GTMHLB ( NL , NC , LABEL ) 

C ALL . LBL I TP ( LABEL , LB 1 , NX (K1 > ) 

PCS ( NDOT ) *0 . 

DO 229 JX=1 , MX 
IF(NN(JX).LE.0) GOTO 229 

PSG C NDOT ) * PSG ( NDOT ) + FLOAT C MM ( JX ) ) /NP I XELYTH2 ( NN ( JX ) , NX ( JX ) ) 

229 CONTINUE 

C 

c 

c COMPUTE SEGMENT UARIANCE 

241 UU-0. 

DO 251 ; . * 1 , MX 
IFCNN(K).LE.l) GOTO 251 

UU=UU+ ( FLOAT ( MM ( K ) )XNPIXEL>*#2 *TH2(NN(IO .NX(K) ) 


t. Y (1 . -TH2 ( NN ( K ) , NX ( K ) ) ) / ( NN l K ) - 1 . > 

4S1 CONTINUE 

IF(JFLAG.EQ.l) WRITE(6.255)ND0T.K1, ID0T.NN(K1),NX(K1) 

:Y « NL , NC , L.AEiEL , LB 1 , UU , PSG C NDOT ) 

25S FORMAT'. 2X. 13, 3X, I3.1X. 16, IX, I3.2X, 13. 2X.214, 2X,2I4, 2F11.6) 
C 

C MODIFY QUADRATIC PRIOR PARAMETERS UPON 

C FINISHING ASSIGNMENT OF INITIAL DOTS 

C 


& 


IFCNDOT.NE.NDS) GOTO 260 
A-6. 

PX-PSGCNDOT) 

IF C0.211.GT.PX) PX-0.211 
IF CPX.GT. 0.789) PX-0.789 
B - 12.*<PX-1. ) 

C • S.-6.*PX 
CALL THPPAR(A.B.C) 

IFCJFLHG.EQ. 1) WRITE(6,257)ND0T«PX*A,B»C 
2S7 FORMAT ( ’ AFTER ASSIGNING’ , 14. ’ DOTS AND ESTIMATE •’.F10.5,/', 

* ’ RESETTING PRIOR PARAMETERS ‘ A«’,F9.5,’ B-’.F9.5.’ C«’.F9.5) 

260 COST I NUE 
C 

C END OF MODIFYING PRIOR PARAMETERS 

C 

IF(NDOT.LT.MAXDOT) GOTO 191 

C 

DO 261 K-l.MX 
IF(NNCK).LE.-l) NNCK)-0 

261 CONTINUE 


C 

C 

C EISD OF BAYES DOT ALLOCATION 

C 

C 

C WRITE CLUSTER INFORMATION 

iFCJKLAG.EQ.l) WRITEC6, 2B2)MAXD0T. NPIXELiMX 
282 FORMAT ( // . 5X > ’ TOTAL NO. OF DOTS AUAIL.ABLE -’.16, 

* /,5X. ’ TOTAL NO. OF PIXEL IN THE CLUSTER MAP -’,16, 

X /,5X, ’ TOTAL NO. OF CLUSTERS -’.16, 

* /,SX, ’ CLUSTER CLUSTER NO. OF NO. OF DOTS LABELLED* 

* /,5X, ’ NO. CODE PIXELS ASSIGNED SMALL GRAIN’ ) 


DO 293 K-l.MX 

I F ( JFLAG . EQ . 1 ) WRITE(6,284)K,LL(K) ,MM(K) ,NN(K) ,NXdO 

283 CONTINUE 

284 F0RMATC8X. I2.5X, I4.3X, I6.4X, 13. 10X, 13) 

IF(MX.EO.M) GOTO 288 

KX-0 

DO 28S K-l.M 

IFCNN(K) .NE.0) GOTO 285 

KX-KX+1 

NDARY(KX) - K 

285 CONT I NUE 

IFCJFLAG.EQ. 1) WRITE(fc,286)MX, (NDARY'K) . K-l, KX) 

286 FORMAT THOSE SMALLER CLUSTERS GROUPED (AS CLUSTER’ , 13, 

* ’ ) ARE:’ , 3(/> 16 C 13, ’ ■ ’ ) ) ) 

283 CONTINUE 
C 
C 
C 


C 

901 WRITE(6,902)M 

902 FORMAT ( ’ YOU HAUE’ , 13, ’ TOO MANY CLUSTERS (MAX -50) ’ ) 
GOTO 990 

911 UPITE(6.912)HD0T 

312 FORMAT C ’ YOU SPECIFIED’ , 14. ’ TOO MANY DOTS C MAX = 200) ’ ) 

C 

C ENTRY TO DEFINE INTIAL NO. OF DOTS TO BE ASSIGNED 

ENTRY INITS(NNIT) 

INIDOT-NNIT 
RETURN 
END 



990 
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FUNCTION TH2CN.NX) 

- FUNCTION TO EUALUATE CLUSTER PROPORTIONAL ESTIMATE 

- N-NO. OF DOTS ASSIGNED FOR A CLUSTER. 

- NX- NO. OF DOTS LABELLED AS SMALL GRAIN 

- BOTH N AND NX ARE SIMPLE INTEGER (2 BYTES IN PDP) 

- TH2 (REAL NO.) CONTAINS THE ESTIMATE ON RETURN. 


INTEGER N. NX 
REAL TH2 


DATA A/6. 31924/. B/-8. 19799/. C/2. 99258/ 

DATA A/6. /.B/-?. 87?/. C/2. 9345/ 

X-NX 
TH2 - 

( A*(X+i. >*(X+2. )#(X+3. > + B*CX+1. )*CX+2. )*(N+4. ) 

+ CXCX+i. )*(N+3. )*<N+4. ) ) / 

< AKCX+l. )*<X+2. >#(N+4. > + B*(X+1. )*(N43. )*(Nf4. ) 
+ CKCN+2. )*<N+3. )*(N+4. ) ) 

RETURN 

- ENTRY TO ALLOW CHANGING OF PARAMETERS 

ENTRY TH2PAR ( AX . BX . CX) 

A-AX 

B-BX 

C-CX 

RETURN 

END 




11. A91 : BAYESIAN DOT ALLOCATION (ADAPTIVE PRIOR) 


11.1 DESCRIPTION 

Implemented in this piece of software is a scheme in which dot 
allocation is sequentially Bayesian (adaptive prior) and labeling 
is direct from ground truth. This scheme is similar to those in 
A84, A85, A87, and A89, except that here estimates based on two 
priors are computed for the initial dots, then one prior is chosen 
and all estimates are based on the selected prior. Such a scheme 
improves the estimate if the small-grain proportion is less 
than 0.21. 


In the Bayesian dot allocation scheme, dots are allocated sequen- 
tially to those clusters that may reduce the variance of the seg- 
ment proportion estimate the most. In order to do this, we must 

2 

first compute the expected change in variance, Ac , for each 
cluster. The expected change in variance for cluster i is 
defined as 





) 


5n i' x i 0n i +1 ' x i +1 * 1 ” 6 n i +l,x i +l^ 
— 


(1 " 6 n v > 0 n +1 * (1 " °n +1 * > 


(li-D 


where N^ ■ number of pixels in cluster i 

N = total number of pixels in the entire scene 

n^ - number of dots previously allocated to cluster i 

x^ * number of dots previously allocated to 
which are labeled as small grain 

(notice x^^ < n^) 


& 


cluster i 


A 

For quadratic prior, 0 ft x ia defined aa 

• 1IK - M i IIS t i ) t ill : US i il l s i H t )l ii i iia m 

For initial dots* a ■ 6, b « -7.877, c ■ 2.9345 


(l.L-2) 


After initial dots* a ■ 6, b ■ 12 (£„ - 1), c ■ 5 - 6P „ if 

ag sg 

0.211 < £ < 0.789 

*9 

a - 6, b - 12(0.211 - 1) , c - 5 - 6 x 0.211 
if 0.211 > # | 
a - 6, b - 12(0.789 - 1) , c - 5 - 6 x 0.789 
if 0.789 < $ 

sg 

where P g g * estimate of small-grain proportion after aasigning 
initial dots 


sg 


For exponential prior, §_ v is defined as 

Ji f 


x + 1 - a x + 1 ,, a wi ,i a * ,, a x 

.»!•» — ' 7 (1 11 s' (l m 1 (1 sm* 


0 


n,x 


l - (1 ' rrr M1 - - x-fi ’ 11 - rrr> 


1 - 4$ 


where a 


S3 


1 - 2P 
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(11-4) 


sg 


P sg * sstirost® of small-grain proportion after assigning 
initial dots 


For the initial dots assignment, estimates of small-grain propor- 
tion are computed on the basis of both priors. After all initial 
dots are assigned, the estimate based on the exponential prior is 
used to determine which prior is going to be used. If it is 


larger than 0.21, then the quadratic prior is selected; other- 
wise the exponential prior is selected. 


2 

Next, a dot is allocated to the cluster whose Ao^ is the largest. 
Then, for this chosen cluster, n^ is updated to n i + 1. That 
dot's label is read from the ground truth file. If the label is 
a small grain, x^ is updated to x^ + 1. After this, the segment 
variance, o 2 , will be computed as 


o 


2 



(1 - 9 . 


, x i n < 


n i - 1 



(11-5) 


where m ■ total number of clusters 


2 

This 0 is printed and the above allocation scheme is repeated 
until the number of dots allocated reaches a value specified by 
the user. 


The proportion estimate based on the quadratic prior is 
computed as 


m N, 


> ■ r -ie 

sg n n L , 


(11-6) 


The proportional estimate based on the exponential prior is 


A 

P 

sg 


m 

E 

i-1 




This scheme, like other allocation and labeling schemes, may be 
repeated with different pseudorandom sequences. The main program 
resets the random number generator to a different starting point 
for each repetition. 

All computation of statistics such as mean and variance is per- 
formed in the main program. A grand summary is printed. 



11. « LINKAGES 


The PDP 11/45 Fortran library (for functions and subroutines 
such as SQRT) , the PDP 11/45 system library (for routines such 
as the pseudorandom number generator , and the utility package 
(described in section 14) are required. 


Main program 
or subprogram 

Subroutine or function required 

MAIN 

SAEB5, GETDOT, CLMPCS, GTMPLB, LBLITP 

SAEB5 

GETDOT, CLMPCS, GTMPLB, LBLITP, TH3 

getdot 

RAN 

TH3 

None 

INTERFACE 



Interface with other routines is through the common block PRTFLG 
(2 bytes) , which is used to control the optional printing of dot 
files and other information. 

11.4 INPUTS 

Fortran- formatted input of the following parameters is needed: 


Card 

or 

line 


Parameters 

Format 

Default 

input 

Default 

value 

1 

Cluster map file name 

A13 

None 


2 

a. 

Number of repetitions 

13 

0 or blank 

1 


b. 

Starting point of first 
pseudorandom sequence 

15 

0 or blank 

10 


c. 

Number of repetition 
printings 

13 

0 or blank 

5 

3 

a. 

Maximum number of dots 
that can be allocated 

13 

None 






b* Number of dots initially 13 
assigned to each cluster 

None 


4 

Number of status messages on 13 
terminal 

0 or blank 

NO MSG 


Two disk files are needed, one containing the stripped cluster 
map, the other the stripped ground truth map. These files should 
have been created by conversion programs A81 and A86. 


11.5 OUTPUTS 

In addition to the normal echo of input specifications and status 
messages on the terminal, the following will be output on a line 
printer. 


report 

Individual repetition 


Grand summary 


Contents 

1. Processor header 

2. Ground truth small-grain proportion 

3. A dot file showing dots chosen, their 
x-y position, their ground truth labels, 
segment variance, and proportion esti- 
mates based on both quadratic and 
exponential priors 

4. A table showing number of clusters, 
cluster name codes, cluster sizes, and 
number of dots assigned 

1. Processor header 

2. A table showing bias, MSE, reduction 
in MSE, average, variance, and vari- 
ance reduction of the estimate based 
on the selected prior as functions of 
number of dots 


11.6 BRIEF ALGORITHM? FOR BAYESIAN DOT ALLOCATION 
(ADAPTIVE PRIOR) 


1. Type messages at terminal and read from it run specifications. 

2. Compute ground truth small-grain proportion. 

3. Repeat a through h until all repetitions are finished: 

a. Set starting point of the pseudorandom number generator. 

b. Clear a dot counter. 

c. Increment the dot counter. 

d. If the dot counter indicates that the current dot should 

be allocated as an initial dot, compute estimates based 

on both priors, then go to the next step. Otherwise com- 
2 

pute Ac . for every cluster and choose the cluster with 

1 a 2 

maximum Ao^ . 

e. Allocate a dot to the chosen cluster. 

f. Pick up the small-grain labels, compute segment variance 
and proportion estimate, and print a line of the dot file. 

g. If this is the end of initial dot assignment, do the 
following: 

(1) Determine if current exponential prior estimate 
exceeds 0.21, then call SETPRI to select to compute 
exponential prior estimate and call SETEP to set the 
value of a. 

(2) Otherwise, call SETPRI to select to compute quadratic 
prior estimate, then call TH3PAR to modify the prior 
parameters . 

h. If the dot count does not exceed specified value, go to c. 

4. Compute bias, MSE, reduction in MSE, average, variance, vari- 
ance reduction, and print a grand summary. 
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Figure 11-1.— Concluded. 
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11.7 LISTING 


C PROGRAM A18 (A91.TSK) ' PROPORTION ESTIMATE OF SMALL GRAIN 

C USING EAVES (ADAPTIUE PRIOR) DOT ALLOCATION 

C LABELLING IS DIRECT BY GROUND TRUTH. 

C 

C ADAPTIUE PRIOR MEANS 

C WHEN ESTIMATE IS GREATER TH3N 0.21. USE QUADRATIC PRIOR 

C OTHERWISE USE EXPONENTIAL PRIOR 

C 

c 


BYTE NAM£(1S).NGT(13> 

REAL PSG( 150) • PM( ISO) .PB( 150) 

C COMMON BLOCK IB A FLAG FOR PRINT ENHBLE/D I SABLE 

COMMON /PRTFLG ' JFLAS 

DATA NGT(5)/’.’/,NGT(6)/’S’/,NGT(7)/’T’/,NGT(8)/’P’/ 

C 

C 

C READING INPUT FROM TERM I NAL 

WRITECS, 121 ) 

121 FORMAT ( ’ PROGRAMS A18 (A91.TSK) . ’ 

* /’ PROPORTION ESTIMATION OF SMALL GRAIN’ 

* /’ COT ALLOCATION IS BAYES (ADAPTIUE PRIOR).’ 

* /’ LABELuING IS DIRECT BY GROUND TRUTH.’ 

* /’ INPUT CLUSTER MAP FILENAME*/’ AAAAAAAAAA AAA’ ) 

READ C 7 . 122 ) (NAME ( K > . K « 1 , 1 3 ) 

122 FORMAT ( 13A1 ) 

C 

WRITE(8. 131) 

131 FORMAT ( ’ AN ASSIGNMENT OF PSEUDORANDOM DOTS IS A REPETITION RUN’ 

* /’ HOW MANY REPETITION RUNS? SEED FOR FIRST PSEUDO SEQ?’ 

* /’ HOW MANY PRINTINGS UF INDIUIDUAL RUNS?’ 

* /’in iiiii nr) 

READ (7. 132) JU. JSEED. JPAGE 

132 FORMAT < 13. IX, IS. IX, 13) 

IF( JU.LE.0) JU-1 
IFCJSEED.LE.0) JSEED- 10 

I F ( JPAGE . LE . 0 ) JPAGE- 5 
JSK I P- ( JU- 1 ) /JPAGE+ 1 
C 

UP I TE ( 8 * 141) 

141 FORM AT (’ SPECIFY MAXIMUM NO. OF DOTS IN THE REPETITIONS RUNS,’/ 

* ’ AND NO. OF DOTS I NT I ALLY ASSIGNED TO EACH CLUSTER.’/ 

* ’ 1 1 1 1 1 1 ’ ) 

READ (7,142) MAXDOT , NN I X 

i<?2 FORMAT ( 1 3, IX, 13) 

CALL INITS(MNIX) 
r . 

WRITECS, 151) 

151 FORMAT (’ STATUS MESSAGES ON TERMINAL, HOW MANY?’/’ Ill’) 

READ (7, 152)NMES 

152 FORMAT (13) 

C 

WRITECS, 101 ) ( NAME ( R ) . K-l, 13) . JU, JSEED, MAXDOT 
181 FORMAT ( /<-' ’ CLUSTER FILE GIUEN - ’.13A1, 


* /’ NO. OF REPETITION RUNS FOR EACH TOTDOT "’,13, 

* /’ THE FIRST REPETITION RUN STARTS WITH SEED -- ’ , IG, 

* /’ MAXIMUM NO. OF DOTS IN THE REPETITION RUNS-’, 16) 

C 

C COMPUTE BAYES (ADAPTIUE PRIOR) ESTIMATE FOR GROUND TRUTH 

I ILINE- 117 


NCOL-196 
IX) 215 K-l, 4 
215 NGT(K)-NAMEOO 

OPEN ( UN IT » £ , NAME - NIST , TYPE = ’ OLD’ , READONLY, FORM- ’ UNFORMATTED’ , 
* ACCESS- ’ DIRECT’ ) 

DO 231 L-l.NLINE 
DO 221 K-l.NCOL 
CALL GTMPLB C L , K , LABEL ) 




CALL LBLITPtLABEL.LBl. IP) 

221 CONTINUE 

231 CONTINUE 

P • FLOAT ( I P ) XNL I NE/NCOL 
C 

OPENCUNIT-l, NAME-NAME. TYPE- ’OLD’ .FORM-’ UNFORMATTED* . 

* READONLY. ACCESS- ’DIRECT* ) 

CALL CLMPLC(NLINE.NCOL) 

C 

C START EST I MAT I CSS FOR AN TOT DOT AND REPETITIONS 

DO 311 K-l.MAXDOT 
PBCJO-0. 

311 PMCK)-0. 

JMES-0 
JS-JSEED 
DO 361 J-l.JU 
JFLAG-0 

IFCMOD(J-1.JSKIF).EQ.0) JFLAG-i 
IF(JFLAG.EQ.l) WRITEC6.321)NAME.P 

321 FORMAT* *1* . 10X, 'ACCURACv ACCESSEMENT SOFTWARE (6- 19-79) ’ /. 2X. 

* ’ PROGRAM A18 CA91.TSK) > ESTIMATION OF SMALL GRAINS PROPORTION’ 

* /.1SX. 'DOT ALLOCATION* BAYES (ADAPTIUE PRIOR) . ’ 

* /.1SX, ’DOT LABELLING 'DIRECT BY GROUND TRUTH. ’ 

* X/.10X,’ INPUT CLUSTER MAP IS FILE ’.15A1, 

* /.10X.’ PROPORTION ESTIMATE FOR GROUND TRUTH IS',F8.S) 
IFCJFLAG.EQ.l) WRITEC6. 322) J, JSEED 

322 FORMAT C/.2X. ’ REPETITION RUN-’. 13, 

* ’ RANDOM DOT SEQUENCE STARTS WITH’ , 18. * ’ ) 

CALL RANSTC JSEED) 

CALL SAEBS(MAXDOT.PSG) 

JMES-JMES+1 

IFCJMES.LE.NMES) WRITE(8.323)MAXD0T. J, JSEED. PSG(MAXDOT) 

323 FORMAT ( ’ MAXDOT* ’ , 13, ’ REPETITION RUN- ’, 13, ’ SEED-’. 16, 

* ’ LAST ESTIMATE- ’.F10. 6) 

JSEED- JSEED+150 

C COMPUTE BIAS AND M.S.E. 

DU 341 K-l.MAXDOT 
TEMF-PSG(K)-P 
PBCK) -PBlK)+TEMP 
341 PMCK)-PM(K1+TEMP**2 
361 CONTINUE 
C 

C PRINT GRAND SUMMARY FOR THIS JOB 

WRITEC6.32DNAME.P 

WR I TE ( C , 37 1 ) MAXDOT . JU . JS 

371 FORMAT (/, 14X. ’ - - GRAND SUMMARY OF THIS JOB ’. 

* /,SX,* MAX. NO. OF DOTS IN EACH REPETITION RUN-’. 13. 

* 14X, ’ NO. OF REPETITION RUNS -M3, 

* /, 14X, ’ RANDOM DOTS START WITH SEED-’, 16. 

% //,1X,’ DOT BIAS M.S.E. REDUCTION’, 

% ’ AUERAGE UARIANCE UAR REDUCTION’ ) 

DO 381 K-l. MAXDOT 
PB(K)»PE(K)/JU 
PMCK)-PM(K)/JU 
RR»PM(K)/( P*(l.-P)/K ) 

AUERG ■ PB C K ) +P 

IFCJU.EQ.l) UAR- PM(K)-PB(K)**2 

IFCJU.GT. 1) UAR- (PM(K)~PB(K)**2)*JU/(JU-i) 

REDUAR-UAR/C P*(l.-P)/K ) 

.381 WRITEC6,382)K,PB(K),PM(K).RR. AUERG. UAR.REDUAR 

382 FORMAT (2X. 13. 2X.F9.5, 4CF10.6), 3X.F10.6) 

C 

c 

WRITEC6.401) 

401 FORMAT ( ’ 1 END OF THIS JOB ’) 

C 

STOP 

END 


fi 
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C 

C 

c 

SUBROUTINE SAEB5 ( MAXDOT , FSG ) 

C STRATIFIED AREAL ESTIMATION USING BAYES DOT ALLOCATION. 

C FORMULAS HERE ARE FOR BAYES WITH ADAPT I UE PRIOR. 

C GIUEN ANY MAX NO. OF DOTS IN MAXDOT. PROGRAM WILL RETURN 

C THE PROPORTION ESTIMATE IS IN ARRAY PSG. 

C 

C WRITEN AND EDITED BY N.Y. CHU ON 6-19-79. 

C 

c 

I NTE6ER NDARY ( 1S0 ) , MM ( 51 ) . MN < SI ) . NX < 51 ) , LL( SI ) 

REAL UU.PSG'l) 

BYTE LABEL 
INTEGERS NMPR 

C COMMON BLOCK IS A FLAG FOR PRINT ENABLE/DISABLE 

COMMON -PRTFLG/JFLAS 
DATA NSTART/0/, M1HSZ/S/ 

DATA INI DOT/3/ 

DATA JFP/1/, JINIT/0/ 
r 

I F MfiXDOT . GT . 150 ) GOTO 911 

C 

C GET PIXEL COUNT FOR EACH CLUSTER 

I F ( NSTAPT . NE . 0 ) GOTO 121 
NSTART-1 

CALL CLMPCSCNPIXEL.M.MM.LL) 

IFOt.GE.Sl) GOTO 901 
121 CONTINUE 
C 

C BEGIN BAYES DOT ALLOCATION 

C 

IFCJFLAG.EQ. 1) WRITEC6.141) 

141 FORMAT C // > 25X, * DOT FILE ’ 

* /’DOT CLUSTER INFORMATION POSITION GROUND TRUTH’ 

* ,’ SEGMENT PROPORTION ESTIMATE’ 

* /’ NUMBER NO. PIXEL DOTS S.G. LINE, COL RAW, CODE’ 

* , ’ UARIAHCE QUAD EXP’ ) 

C 

r MAKE Si IRE FATM n WATER HAS A MINIMUM NO. OF PIXELS 

' ' ........ 

NDS-0 

DO 151 K-l.M 

NXCK)-0 

NNCK)=0 

I F ( MM t K ) . GE . M I NS2 ) NDS-NDS+ 1 N I DOT 
IF(MMCIO.LT. MINS2) Kl-Kl+MMCK) 

IFfMMCK) .LT.MINSZ) NNCK) ~-l 
1S1 CONTINUE 
MX-M 

IFCK1. LT.MINSZ) GOTO 161 
MX-M+1 
MM (MX) =1(1 
NDS-NDS 4-INI DOT 
161 CONTINUE 
C 

IFC..TFP.F.Q. 1) NMPR- ’QUAD’ 

IFCJFP.EQ.2) NMPR- ’EXP ' 

IFCJFLAG.EQ. 1) WRITEC6, 17UNDS, NMPR 
171 FORMAT C ’ I NI T I AL ASS I GNMENT ' ’ , 14, ’ DOTS . ’ , 

* ’ SEC! WAR IS BASED ON ’,A4,’ PRIOR.’) 

C 

C FROM HERE TO STATEMENT 261 « LOOP TO ASSIGN 1 DOT AT A TIME 

C 

C IN ORDER TO ALHIEUE ASSIGNING 3 DOTS OR LESS INITIALLY TO EACH CLUSTER 

C THE DELIA UARIHNCE COMPUTATION IS BYPASSED FOR THE FIRST 

C NDS=INIDOT*MX DOTS 

C 


J. i-lT 

n 


KG1-0 

NDOT-0 

191 NDOT-NDOT+1 

IFCNDOT.GT.NDS) SOTO 200 
195 KG1-KG1+1 

IFCKG1.GT.MX) KG1-1 
K1-KG1 

IFCNNCKGi) .GT.-l) GOTO 214 
GOTO 195 
C 

C FOR EACH CLUSTER. COMPUTE DELTA VARIANCE AND CHOOSE THE LARGEST 

200 CONTINUE 

C SET THE PRIOR TO BE USED 

CALL SETPRHJFP) 

AMAX-0. 

DO 211 K-l.MX 
J-NN(K) 

JX-NXCJO 

IF(J.LF.l) GOTO 211 
C 

UAR- ( FLOAT (MM CIO )/NPIXEL )**2 

* *C TH3CJ, JX)*C1.-TH3(J. JX)VCJ~1. ) 

* -TH3 C J . JX ) *TH3 C J+ 1 . JX+ 1 ) * ( 1 . -TH3 ( J+ 1 . -1X+ 1) V J 

* -C1.-TH3CJ, JXmTH3CJ+l, JX)*Cl.-TH3CJ+i,JX»/J ) 

C 

C URITEC6. 99121 )K.MMCK) . J-NXCIO .UAR 

99121 FORMAT ( ' K. MMCK) . J.NXCK) .UAR-' . 14. 10,215. E16. 6) 

IFCUAR.LT. AMfiX) GOTO 211 
Ki-K 

AMAX-UAR 
211 CONTINUE 
C 

C ASSIGN A DOT TO THE CLUSTER WITH MAX DELTA VARIANCE 

IFCNNCKD.GE.MiCKll) WRITEC6.213)ND0T.K1.MMCK1) 
IF(NNCKl).GE.MMCKl)) GOTO 241 

213 FORMAT C ’ ATTEMPT TO ASSIGN DOT"’ . 13. ’ WHILE NO PIXEL’. 

X ’ IS AVAILABLE. CLUSTER- ’. 13. ’ TOT PIXEL- ’. 16) 

C POSITION THE POINTER TO NDARY 

214 ND-1 

DO 215 K-1.K1 
IFCNNOO.LE.0) GOTO 215 
ND-ND+NNCK) 
rr,^T t if 

C MAKE MjON Fuk 1H£ haw Dui IN NDhkY 

Ml— 1 

DO 217 K-NDOT.ND.N1 
217 NDARY C K+l ) -NDARY C K ) 

ND-ND-NNCK1) 

NNCFJ )«NN(K1>+1 
C ASSIGN A DOT 

CALL GETD0TCMMCK1).NNCK1),NDARYCND). IDOT) 

C WRITECb.99141)Kl.ND.NNCKl).NXCKl). CNDARYCK) . K*1 . 44J 
99141 FORMATC’ K1.ND.NN.NX-’ .414. 4C/.2X.15I4) ) 

IFCKi.LE.M) GOTO 227 
DO 225 Kl-l.M 
IFCNNCK1) .GT.0) GOTO 225 
IFC IDOT.LT.MMCK1 ) ) GOTO 227 
IDOT* IDOT-MMCK1 ) 

225 CONTINUE 

227 CALL CLMPXY C K1 . I DOT . NL , NC . LD 
CALL GTMPLB C NL . NL- . LABEL ) 

CALL LBLITPCLABEL.LR1.NXCK1)) 




C COMPUTE PROP EST DUE TO QUADRATIC PRIOR 

u COMPUTE PROP EST HUE TO EXPONENTIAL PRIOR 

Pl-0. 

P2-0. 

DO 231 JX-l.MX 
IF(NN(JX).LE.0) GOTO 231 

C AMAX.Fi.Nl ARE TEMPORY VARIABLES 

AMAX- FL OAT ( MM ( JX ) )/NP I XEL 
Ml -NN(JX) 

J-NX(JX) 

C 

C FIRST. ACCUMULATE REST WITH QUAD PRIOR 

1FCNDOT.LE.NDS) CALL SETPRI(l) 

IF(NDOT.GT.NDS) CALL SETPRI ( JFP) 

PI- PI + AMAX*TH3(N1. J) 

C 

C NEXT. ACCUMULATE EST WITH EXP PRIOR 

IFCNDOT.GT.NDS) GOTO 231 
P2- P2 + AMAXXJ/N1 
231 CONTINUE 
C 

C STORE THE APPROPIATE PROP EST IN THEIR CORRECT PLACE 

IF(NDQT.GT.NDS) GOTO 233 
PSG( (NDOT-1 )*2+l )-Pl 
PSG< ND0TX2 ) -P2 
GOTO 235 

233 PSG(NDOT) -Pi 

IF(JFP.EQ.l) P2-0. 

I F ( JFP . EQ . 2 ) P2-P1 
IF( JFP. EQ. 2) Pl-0. 

235 CONTINUE 
C 

c 

C COMPUTE SEGMENT VARIANCE 

241 UU-0. 

CALL SETPRI (JFP) 

DO 251 K-l.MX 
lF(NN(IO.LE.l) GOTO 251 

UV-UV+ ( FLOAT ( MM ( K ) ) /NP I XEL ) XX3 *TH3 ( NT K K ) . NX ( K) ) 

* * ( 1 . -TH3 ( NN ( IO < NX ( K ) ) ) Y ( NN ( K) -1 . ) 

251 CONTINUE 

IF(JFLHG.EQ.l) WRITEC6, 25S1ND0T. K1 . I DOT . NN( K1 ) » NX (Kl ) 

* . NL , NC . LABEL . LB 1 . UU , 1 ’1 , P2 

255 FORMAT (2X. 13. 3X. 13. IX. 16. IX. 13. 2X. 13. 2X.2I4. 2X. M, 15. 3F10.b> 
C 

C DETERMINE WHICH PRIOR TO USE 

IF(NDOT.NE.NDS) GOTO 275 

C 

C CHOOSE PRIOR ACCORDING TO VERY FIRST REPETITION OF THE RUN 

IFC ( JINIT. NE. 0) . AND. (JFP. EQ. 1 ) ) GOTO 2S6 
IF( ( JINIT. NE. 0) . AND. ( JFP. EQ. 2) ) GOTO 259 
JINIT-1 


C CHOOSE PRIOR ACCORDING TO ESTIMATE FOR THE FIRST REPETITION 

c use THE EXP PRIOR ESTIMATE TO DECIDE 


PX-PSG( ND0T*2 ) 

IF(PX.LE.0.?1> GOTO 259 

C 

I QUAD PRIOR IS CHOSEN. NEED TO MOD I FY QUAD PRIOR PARAMETERS 

256 PX-PSG( (NDOT-1) *2+1 ) 

JFP- 1 

NMPR • ' QUAD ’ 



C 




A-6. 

IF (0.211.GT.PX) PX-0.211 
IF (PX.GT. 0.789) PX-0.789 
B • 12.*CPX-1. ) 

C - 5.-6.XPX 
CALL TH3PAR(A.B.C) 

IF( JFLAG.EQ. 1) WRITE(6,2S7)A.B.C 
257 FORMAT ( ’ RESETTING’. 

* ’ QUAD PRIOR PARAMETERS' A-’,F9.5. * B- .F9.S.’ C-’,F9.S> 

C 

GOTO 265 
C 

C EXP PRIOR IS CHOSEN, NEED TO SET EXP PRIOR PARAMETERS 

259 PX-PSG(NDOT*2) 

1FCPX.LE.0.01 ) PX-0.01 
JFP-2 

NMPR-’EXP ’ 

IF(PX.LE. 0.249) GOTO 262 
IF( JFl.AG.EQ. 1 ) WRITE(6»261 )PX 

261 FORMAT ( ’ ESTIMATE- ’.F9.6,’ TOO HIGH, BEING RESET TO AUOID’ , 

* ’ NEGATIVE ALPHA. ’ ) 

PX-0.249 

262 CALL SETEPC ( 1 . -4. *PX)/C 1 . -2. *PX) ) 

C 

c STORE THE CHOSEN PRIOR ESTIMATES, THESE WILL BE RETURNED 

26S DO 267 K-l.NDS 

267 PSG(K) -PSG( (K-1)*2+JFP ) 

C 

IFCJFLAG.EQ.l) WRITE(6,271)PX.NMPR.NMPR 
271 FORMAT ( ’ ESTIMATE- ’, F10. 6, ’ . NOW STORE THE ’,A4,’ PRIOR* 

* ’ ESTIMATES...’/’ SEO UAR AND DELTA UAR ARE’, 

* ’ BEING COMPUTED WITH ’,A4.* PRIOR.') 

C 

275 CONT I NUE 

C 

C 

I F < HDOT . LT . MAXDOT ) GOTO 191 
C 

2 77 DO 278 K-l.MX 
IFCNN(K) .LE.-l) NN(K) *0 

278 CONTINUE 
C 

C 

C END OF BAYES DOT ALLOCATION 

C 

c 

C WRITE CLUSTER INFORMATION 

280 CONTINUE 

IF(JFLAG.EQ.l) WRITE(6,282)MAXD0T,NPIXEL.MX 

282 FORMAT (//,5X,’ TOTAL NO. OF DOTS AUHILABLE -’,16, 

* / i 5X , ’ TOTAL NO. OF PIXEL IN THE CLUSTER MAP ’’,16. 

* /, 5X. ' TOTAL NO. OF CLUSTERS -’>16, 

* /• 5X, ’ Q.USTER OUSTER NO. OF NO. OF DOTS LABELLED’ 

* / > 5X » ’ NO. CODE PIXELS ASSIGNED SMALL GRAIN’) 
DO 283 K-l.MX 

IFCJFLAG.EQ.l) WRITE(6,284)K, t-LCK) ,MM(K) . NNCK) ,NX(K) 

283 CONTINUE 

284 F0RMATC8X. I2.5X, I4.3X, I6.4X, 13, 10X, 13) 

IFCMX.EQ.M) GOTO 288 

KX-0 

DO 285 K-l.M 
IFCNN(K) .NE.0) GOTO 285 
KX-KX+1 
NDARYCKXJ *K 


uyt\ 

/# 


OP THE 
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285 CONTINUE 

IF(JFLAG.ED.l) URITE(G,286)MX, (NDARY(K).K-l.KX) 

296 FORMAT O' THOSE SMALLER CLUSTERS GROUPED (AS CLUSTER' . 13. 

* ') ARE» ' » 30, 16(13.', ’)> ) 

288 CONTINUE 

C 

C 

C 

c 

GOTO 990 

C 

c 

901 HRITE(6,902)M 

902 FORMAT < ’ YOU HAUE’ , 13, ’ TOO MANY CLUSTERS (MAX- 50 ) ’ > 

GOTO 990 

911 WRITE(6.312>ND0T 

912 FORmTC YOU SPECIFIED’ , 14, ’ TOO MANY DOTS (MAX -150) ' ) 

C 

c 

i.: ENTRY TO DEFINE INTIAL NO. OF DOTS TO BE ASS I UNEP 

ENTRY INITS(NNIT) 

INIDOT-NNIT 
990 RETURN 
END 
C 
C 
C 

FUNCTION TH3(N,NX) 

C FUNCTION TO EVALUATE CLUSTER PROPORTION ESTIMATE 

c FOR the CASE (1) QUADRATIC PRIOR IS USED 

c (2) EXPONENTIAL PRIOR IS I SED 

C N«NO. OF DOTS ASSIGNED FOR A CLUSTER. 

C NX- NO. OF DOTS LABELLED AS SMAl L GRAIN 

C BOTH N AND NX ARE SIMPLE INTEGER (2 BYTES IN FDP) 

C TH3 (REAL NO.) CONTAINS THE ESTIMATE ON RETURN. 

C 

INTEGER N, NX 
REAL TH3 
C 

DATA A/6 . / , B/-7 . 877/ , C/2 . 9345/ 

C 

C CHOOSE THE DESIRED PRIOR 

IFCJFP.EU.2) GOTO 200 
C 

L THE FOLLOWING IS FOR QUADRATIC PRIOR 

120 X-NX 
TH3 • 

* ( A*(XU. )*(X+2. )*(X+3. ) + B*.:X+1. )#(X+2. )t(N+4. > 

* ¥ C4NX+1. )*(N+3. )*(N*4. ) ) / 

* ( A*tXU. )*(X+2. )*(N+4. ) + 8*<X+1. )*CN+3. )Y(N+4. 1 

* + C*(N+2. )*(N+3. )*(N+4. ) ) 

C 

RETURN 

r. 

c 

c THE FOLLOWING IS FOR EXPONENTIAL PRIOR 

200 KX-m+l 
C 

PI • ( 1 . -AL/KX ) 

211 IF(KX.GT.N) GOTO 221 
KX‘KX+1 

P1*P1*( 1.— AL/KX) 

GOTO 211 

r, 

221 TH3- ( (NX+l.-AL)/(N+2.-AL) - (NX+1. )/(H+2. )*P1 )/ ( 1. - PI ) 

r. 

RETURN 


reproducibility 
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~ ENTRY TO ALLOW CHANGING OF P AR AMET ER S 

ENTRY TH3PAR(AX.BX,CX) 

A»AX 
B-BX 

c-cx 

RETURN 

— ENTRY TO DEFINE ALPHA USED BY EXP PRIOR 

ENTRY 3ETEPC ALPHA) 

AL-ALPHA 
RETURN 

-- ENTRY TO DEFINE WHICH PRIOR TL CHOOSE 

ENTRY SETPRI ( JFPP) 

JFP-JFPP 
RETURN 
C 

END 
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12. A88t BAYESIAN DOT ALLOCATION, MAJORITY RULE LABELING 


12.1 DESCRIPTION 

This dot allocation and labeling schema is different from most 
other schemes in one important respect t the allocation procedure 
is fused with the labeling procedure. 


For each cluster, dots are allocated continuously until the 

cluster is labeled. The rule is Bayesian majority rule, stated 

as follows: 

a. Assign 2 dots; if their labels are identical, then give the 
cluster thrt label and stop. 

b. Assign 3 more dots; if the 5 labels yield a 1-4 split, then 
label the cluster by the majority and stop. 

c. Assign 2 more dots; if the 7 labels yield a 2-5 split, then 
label the cluster by the majority and stop. 

d. Assign 3 more dots; if the 10 labels yield a 3-7 split, then 
label the cluster by the majority and stop. 

e. Assign 3 more dots and label the cluster by the majority. 


It can be seen from the above that the number of dots assigned to 
the cluster may vary, depending on where the allocation and label- 
ing is stopped. Also, it should be noted that the labels above 
are either small grain or other. 

The equation for the proportion estimate is 



£ 

cluster i 
labeled as 
small grain 



( 12 - 1 ) 


where N i - number of pixele in cluster i 

N ■ total number of pixels in the entire scene 

This scheme, like other allocation and labeling ‘schemes, may be 
repeated with different pseudorandom requences. The main program 
resets the random number generator to a different starting point 
for each repetition. 

All computation of statistics such as mean and variance is per- 
formed in the main program. A grand summary is printed. 


I 2 * 2 LINKAGES 

The PDP 11/45 Fortran library (for functions and subroutines such 
as SQRT) , the PDP 11/45 system library (for routines such as the 
pseudorandom number generator) , and the utility package (described 
in section 14) are required. 


Main program 

or subprogram Subroutine or function required 


MAIN 

SAEBMR 

GETDOT 

BMR1 

(Bayesian majority 
rule 1) 

BMRSUP 

(Bayesian majority 
rule support) 


SAEBMR, GETDOT, CLMPCS, GTMPLB, LBLITP 

CLMPCS, BMR1 

RAN 

BMRSUP 


GETDOT, CLMPCS, GTMPLB, LBLITP 


12.3 INTERFACE 

Interface with other routines is through the common block PRTFLG 
(2 bytes), which is used to control the optional printing of dot 
files and other information. 


/*•> 


12.4 INPUTS 


Fortran- formatted input of the following parameters is needed: 
Card 

or Default Default 


line 


Parameters 

Format 


input 

value 

1 

Cluster map file name 

A13 

None 


2 

a. 

Number of repetitions 

13 

0 

or blank 

1 


b. 

Starting point of first 

15 

0 

or blank 

10 



pseudorandom sequence 






c. 

Number of repetition 

13 

0 

or blank 

5 



printings 





3 

Number of status messages on 

13 

0 

or blank 

NO MSG 


terminal 


Two disk files are needed, one containing the stripped cluster 
map, the other the stripped ground truth map. These files should 
have been created by conversion programs A81 and A86. 


12.5 OUTPUTS 

In addition to the normal echo of input specifications and status 
messages on the terminal, the following will be output on a line 
printer. 


Report 


Contents 


Individual repetition 1. Processor header 

2. Ground truth small-grain proportion 

3. A dot file showing dots chosen, their 
x-y position, and their ground truth 
labels 


4. A table showing number of clusters, 
cluster name codes, cluster sizes, 
number of dots assigned, and the split 


Grand summary 


1 


Processor header 


2. A table showing the estimate and bias 
of each repetition 

3. Bias, MSE, reduction in MSE, average 
estimate, variance, variance reduction, 
average and standard deviation of the 
total number of dots per repetition 

12.6 BRIEF ALGORITHM: FOR BAYESIAN DOT ALLOCATION, 

MAJORITY RULE LABELING 

1. Type messages at terminal and read from it run specifications. 

2. Compute ground truth small-grain proportion. 

3. Repeat a and b until all repetitions are finished: 

a. Set starting point of the pseudorandom number generator. 

b. For each cluster, allocate and label dots by Bayesian 
majority rule. 

4. Compute bias, MSE, reduction in MSE, average estimate, vari- 
ance, variance reduction, average and standard deviation of 
the total number of dots per repetition, and print a grand 
summary. 















Figure 12-1.- Concluded. 








t 


12.7 LISTING 


L Alb (HtJb. ISK) i PROPORTION ESTIMATE OF bhAL.L GRAIN 

C USING BAYES MAJ RULE (UNIFORM PRIOR) DOT ALLOCATION 

C LABELLING IS BAYES MAJORITY RULE. 

C 

BYTE NAME(15),N6T<i3) 

REAL PSGC3S0) 

INTEGER NDOTC350) 

C COMMON BLOCK IS A FLAG FOR PRINT ENABLE/D I SABLE 

COMMON /PRTFLG/JFLAG 

DATA NGT(S)/’ . ’/.NGT(S)/’S’/,NGT(7)/’T’/,NGT(8)/’P’/ 

C 

c 

C READING INPUT FROM TERMINAL 

C 

WRITE(S.l21) 

121 FORMAT ( ' FROGRAM« A1S (A88.TSK) . * 

* /* PROPORTIONAL ESTIMATION OF SMALL GRAIN* 

* /’ DOT ALLOCATION IS BAYES MAJ RULE (UNIFORM PRIOR) 

* /’ LABELLING IS ESAYES MAJORITY RULE.’ 

* /’ INPUT CLUSTER MAP FILENAME’/'’ AAfiAAAfi AAAAAA ’ ) 

READ(7i 122) iNAME(K) > K- 1 > 13) 

122 FORMAT C 13A1 > 

C 

URITE(G. 131 ) 

131 FORMAT ( ’ AN ASS I GWENT OF PSEUDORANDOM DOTS IS A REPETITION RUN’ 

* /’ HOW MANY REPETITION RUTS? SEED FDR FIRST PSEUDO SEQ?’ 

* /’ HOW MANY PRINTINGS OF INDIUIDUAL RUNS? 5 

* /’ II! Hill III’) 

READ(7, 1 32 ) JU > JSEED , JPAGE 

132 FORMAT ( 1 3 » IX. IS. IX, 13) 

1FCJU.LE.0) JU=1 
IFCJSEED.LE.0) JSEED=10 
IF(JPAGE.LE.0> JPAGE*5 
JSK IP* ( JU-1 ) /'JPAGE+1 

C 

NRITE(8, 141 ) 

141 FORMAT ( ’ NOTICE THAT NO. OF DOTS ASSGINED IN EACH’. 

* ’ REPETITION MAY UAEY ’) 

i" 

WRITEC8, 1SI ) 

.151 FORMAT (’ STATUS MESSAGES ON TERMINAL. HOW MANY?’/” Ill’) 

READC7, 152)NMES 
152 FORMAT (13) 

C 

WR I TE C B . 1 8 1 ) ( NAME < IO , K - 1 , 1 3 > , JU , JSEED 
181 FORMAT.://’ CLUSTER FILE GIUEN * ’.13A1, 

* /’ NO. OF REPETITION RUNS FOR EACH TOTDOT =’,I3. 

X /’ THE FIRST REPETITION RUN STARTS WITH SEED - ’.16) 

C 

C COMPUTE PROPORTIONAL ESTIMATE FOR GROUND TRUTH 

NLINE=117 
NC0L=196 
DO 21S K-i.4 
21S NGT ( K ) “NAME ( K ) 

OPEN ( UN I T «2 . NAME *NGT > TYPE * ’ 01 .D ’ , READONLY. FORM* ’UNFORMATTED’ . 

* ACCESS* ’ D I RECT ’ ) 

TO 231 L=1.NLINE 

DO 221 K*1.NC0L 
CALL GTMPLB(L,K. LABEL) 

CALL LBLITP(LABEL .LB1 i IP) 

221 CONTINUE 

231 CONTINUE 

P * FLOAT ( I P ) /NL 1 NE/ NCOL 
C 

OPEN ( UN IT * 1 , NAME “NAME , TYPE- ’ OLD ’ .FORM* ’UNFORMATTED’ , 

* READONLY, ACCESS* ’DIRECT’ ) 

CALL CLMPLC(HLINE.HCOL) 



c 

C START ESTIMATION FOR EACH TOTDOT AND REPETITIONS 

JMES-0 
JS-JSEED 
DO 361 J-l.JU 
JFLAG-0 

IF(MODCJ-1,JSKIP).EQ.0) JFLAG-1 
IFCJFLAG.EQ.l) WRITE(6,321)NAME.P 

321 FORMAT < ’ 1 ’ , 10X , ’ ACCURACY ACCESSEMENT SOFTWARE ( S-17-79 ) ’ , / , 2X , 

* ’PROGRAM A1S CA88.TSK>« EST I MAT 1 04 OF SMALL GRAINS PROPORTION* 

* /, 15X, ’ DOT ALLOCATION BAYES MAJ PULE (UNIFORM PRIOR) , ’ 

* /« 15Xi ’DOT LABELLING i EAVES MAJORITY RULE. ’ 

% //.10X. • INPUT CLUSTER MAP IS FILE * ,15A1, 

* /.10X*’ PROPORTIONAL ESTIMATE FOR GROUND TRUTH IS’ .F8.S) 
IF(JFLAG.EQ.l) WRITE<6. 322) J, JSEED 

322 FORMAT ( / , 2X , ’ REPETITION RUN-’ . 13. 

* ’ RANDOM DOT SEQUENCE STARTS WITH* , 18, * ’ ) 

CALL RANST(JSEED) 

CALL SAEBMRCNDOT(J).PSGCJ)) 

JMES-JMES+1 

IF(JMES.LE.NMES) WR I TE ( 8 , 323 ) MAXDOT , J , JSEED 

323 FORMAT ( ’ MAXDOT* ’. 13, ’ REPETITION RUN* ’ , 13, ’ SEED*’, 16) 

JSEED* JSEED+150 

361 CONTINUE 
C 


C PRINT GRAND SUMMARY FOR THIS JOB 

WR I TE ( 6 . 32 1 ) NAME > P 
WRITE<6. 371 ) JU, JS 

371 FORMAT (/.14X,* GRAND SUMMARY OF THIS JOB ’, 

* /> 14X, * NO. OF REPETITION RUNS -’,13, 

* /,14X, ’ RANDOM DOTS START WITH SEED- ’,16,/ 

* ,10X, ’REPETITION PSEDUO SEQ DOT SMALL GRAIN BIAS’,/, 

* 10X, ’ RUN SEED ASSIGNED ESTIMATE WRTG.T. ’) 

C COMPUTE BIAS AND M.S.E. 

PB=0. 

PM=0. 


AM*®. 

UM*0. 

DO 381 J-l.JU 
TEMP=PSG( J)-P 
PB-PB+TEMP 
PM-PM+TEMP**2 
AM-AM+NDQT(J) 

UM-UM+NDQT(J)*#2 

WR I TE ( 6 . 37S ) J , JS , NDOT ( J ) . PSG < J ) , TEMP 
375 FORMAT USX. 13. 6X. 16, 6X. 14, 5X,F8.S,3X,F9.6) 

381 JS=JS+1S0 
PB-PB/JU 
PM-PM/JU 
AM-AM/JU 
UM-UM/JU 

IFCJU.EQ.l) UM-UM - AM*X2 
IF(JU.GT.l) UM- CUM-AM*#2)*JU/(JU-1) 

UM- SORT CUM) 

AUERG-PB+P 

IFCJU.EQ.l ) UAR-- PM-PB**2 
IFCJU.GT.l) UAR- (PM-PB**2)*JU/<JU-1) 

C COMPUTE UARIANCE REDUCTION 

RR-FM/C PY ( 1 . -P ) /AM ) 

REDUAR-UAR/ C P*(1.-P)/AM ) 

WR I TE ( 6 . 38S ) PB , PM , RR , AUERG , UAR , RED JAR , AM , UM 
385 FORMAT </’ BIAS- ’ . F10.6, * M.S.E. •’ .F10.6, ’ REDUCTION-’ , FI®. 6, 
X /’ AUERAGE-' ,F10.6, ’ UARIANCE- F10. 6, ’ REDUCTION-’ , F10. 6, 

X /’ AUE DOT-’ ,F10.3, ’ DOT S. D. - ’ . F10. 3) 

C 
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C 

WRITE(6>401) 

401 FORMAT ( ' 1 END OF THIS JOB ’ ) 

C 

STOP 

END 

C 

SUBROUTINE SAEBMRCNDQT, PSG) 

C STRATIFIED AREAL ESTIMATION USING BAYES MAJORITY RULE DOT ALLOCATION. 

C FORMULA HERE ARE FOR BAYES MAJORITY RULE WITH UNIFORM PRIOR. 

C ON RETURN, THE PROPORTIONAL ESTIMATE IS IN PSG. AND THE 

C TOTAL NO. OF DOTS ASSIGNED IS IN NDOT. 

C 

c UR ITEM AND EDITED BY N.Y. CHU ON 5-17-79. 

C 

C 

INTEGER NDARY ( 350 ) « MM ( 51 ) , NN( 51 ) , L.L ( Si ) 

INTEGER NNSC2.51) 

INTEGER* 4 NSG 

C COMMON BLOCK IS A FLAG FOR PRINT ENABLE/D I SABLE 

COMMON /PRTFLG/JFLAG 
DATA NSTART/0/, MINSZ/13/ 

C GET PIXEL COUNT FOR EACH CLUSTER 

I F ( NSTART . NE . 0 ) GOTO 121 
NSTART ■ 1 

CALL CLMPCS(NPIXEL,M,MM,LL) 

IFCM. GE. 51 ) GOTO 901 
121 CONTINUE 
C 

C MAKE SURE EACH CLUSTER HAS A MINIMUM NO. OF PIXELS 

K1=0 

DO 151 K“1,M 
NN(K)>0 

I F ( MM ( K ) . LT . M I N3Z ) Kl-Kl+MMCK) 

IF(MMCK).LT.MINSZ) NN(K)=-1 
151 CONTINUE 
MX=M 

IF (Kl.LT.MINSZ) GOTO 161 
MX-M+1 
m(MX)=Kl 
161 CONTINUE 
C 
C 

C BEGIN BAYES MAJORITY RULE DOT ALLOCATION 

C 

IF( JFLAG.EQ. 1) WRITERS. 211) 

211 FORMAT t/V DOT FILE ’ 

* /’ CLUSTER DOT NUMBER POSITION GROUND TRUTH LABEL’ 

* /• NO. URT ITS CLUSTER LINE, COL RAM, CODE SMALL GRAIN’ ) 

C ASSIGN DOTS AND COMPUTE ESTIMATE 

NDOT-0 

PSG*0. 

ND*1 

DO 231 K* 1 , MX 
IF(NNCK).LE.-l) GOTO 231 

CALL BMR1 ( ISG, NNSC 1 , K> > NDARYCND) ,K,MM.NN.LL,M) 

PSG > PSG+FLOAT ( I SG ) *MM ( K > /NP I XEL 
NDOT-NDOT+NNCK) 

ND-ND+NNCK) 

231 CONTINUE 

r. 


c 

C END OF BAYES MAJORITY RULE DOT ALLOCATION 

C 

C WRITE CLUSTER INFORMATION 

IFCJFLAG.EQ. 1) WRITEC6,282)ND0T.NPIXEL,MX 

282 FORMAT C//.5X, * TOTAL NO. OF DOTS ASSIGNED -’,16. 

* /,S X,' TOTAL NO. OF PIXEL IN THE CLUSTER MAP -’.16, 

* / » 5X» * TOTAL NO. OF CLUSTERS -’.16, 

* /, 10X. ’ CLUSTER CLUSTER NO. OF NO. OF DOTS’ 

* ^.10X. ’ NO. CODE PIXELS ASSIGNED SPLIT’) 

1X3 283 K-l.MX 

IF(NNOO.LE.-l) NNCK)»0 
IF(NNS(l,lO.GE.NNS(2.K)) NSG-’S.G.’ 

IFCNNSC1,K).LT.NNSC2,K)) NSG- ’ 

IFCJFLAG.EQ. 1) MRITEC6,284)K,LLCK).MM<K).NN(K),NNSC1,K).NNS(2,K) 

* »NSG 

283 CONT I fSUE 

284 F0RMATC13X. I2.5X, I4.3X. I6.4X, 13, 7X, 12, ’ , 12, IX, A4) 

IFCMX.EQ.M) GOTO 2 88 

KX-0 

DO 28S K-l.M 
IF(NNCK) .NE.0) GOTO 285 
KX-KX+1 
NDARY CKX)-K 

285 CONTINUE 

IFCJFLAG.EQ. 1) WR I TE ( 6 , 286 ) MX , (NDARY(K) ,K-1,KX) 

286 FOt MAT ( / ’ L THOSE SMALLER CLUSTERS GROUPED CAS CLUSTER’ , 13, 

* ’) ARE*’. 3C/.16CI3,’,’)) ) 

288 CONTINUE 

IFCJFLAG.EQ. 1) WRITEC6, 351 )NDOT, PSG 
351 FORMAT C r TOTAL DOTS ASSIGNED-’ . 14, ’ PROPORTIONAL’ . 

X ’ ESTIMATE-' ,F10. 6) 

C 

GOTO 990 
C 

901 WRITEC6.9023M 

902 FORMAT C ’ YOU HAUE’ , 13, ’ TOO MANY CLUSTERSCMAX-50) ’ ) 

GOTO 990 

C 

C 

990 RETURN 
END 


C 

C 

C 

SUBROUTINE DMRICISG.NS, NDARY , K1 , MM , NN , LL , M ) 

C BAYES MAJORITY RULE CUNIFORM PRIOR) 

C TOGETHER WITH SUB. BMRSUP, THIS ROUTINE ASSIGN 

C DOTS ACCORDING A BAYES MAJORITY RULE. 

C ON RETURN, NNCK1 ) CONTAINS THE NO. OF DOTS ASSIGNED DURING THIS CALL 

C ISG-1 , ON RETURN, MEANS SMALL GRAIN LABEL IS ASSIGNED TO 

C THE CURRENT CK1) CLUSTER, IF ISG-0, NOT SMALL GRAIN 

C NDARY- DOT ARRAY FOR THE CURRENT CLUSTER 

C MM- ARRAY CONTAINING CLUSTER SIZES 

C NN- ARRAY CONTAINING CLUSTERS’ FUTURE DOT OF ASSIGNED DOTS 

C LL- ARRAY CONTAINING CLUSTERS’ CLASS NAME. 

C JUST FOR USE IN SUB. CLMPXY 

C 

INTEGER NDARY C 1 ) , MMC 1 ) , NNC 1 ) , LLC 1 ) 

INTEGER NSC 2) 

C COMMON BLOCK IS A FLAG FOR PRINT ENABLE/DISABLE 

COMMON /PRTFLG/JFLAG 



C 


NSCD-0 

NSC2)«C 


c 

C BAYES MAJ RULE WITH UNIFORM PRIOR 

CALL KiRSUPC2.NS.MR, NDARY, K1,I1M,NN,LL«M) 
IFCNSCMR) .EQ. 2) GOTO 401 
CALL K1RSUPC3.NS.MR. NDARY.Kl.MM.NN.LL.M) 
IFCNSCMR). EQ. 4) GOTO 401 
CALL K1RSUPC2.NS.MR. NDARY, Kl. MM, NN.t-L.M) 
IFCNSCMR). EQ. 5) GOTO 401 
CALL BMRSUPC3.NS.MR. NDARY. Kl, MM, NN,LL.,M) 
IFCNSCMR). EQ.?) GOTO 401 
CALL BMRSUPC3.NS.MR. NDARY.K1 .rtl.NN.LL.M) 
C 

401 IFCMR.EQ.l) ISG-1 
IFCMR.EQ.2) ISG-0 
RETURN 
END 
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c 

SUBROUTINE BMRSUPC JDOT. NS, MR. NDARY, KOLD. MM, NN.LL.M) 
C TOGETHER WITH SUBROUTINE KIRI, THIS ROUTINE 


C ASSIGN DOTS TO A CURRENT CLUSTER CSPECIFIED IN Ki), AND 

C DETERMINE THE SPLIT BETWEEN SMALL GRAIN AND NON SMALL GRAIN 

C LABELS. 

C JDOT - NO. OF DOTS DESIRED TO BE ASSIGNED IN THIS CALL 

C NS- CDIM 2) ARRAY CONTAINING THE SPLIT 

C FIRST MEMBER REFERS TO SMALL GRAIN, SECOND MEMEBER NON S.G. 

C NOTICE FOR NS, NEW KITS ARE ACCUMULATED TO PREUIOUS UALUES. 

C MR- 1 OR 2. IF -1, MEANS S.G. IS THE MAJORITY. OTHERWISE NON S.G. 

C NDARY- ARRAY NECESSARY FOR DOT GENERATION 

C 

INTEGER NSC2) , NDARY C l),ft1Cl),NNCl),LLCl) 

BYTE LABEL 
INTEGER44 NSG 

C COMMON BLOCK IS A FLAG FOR PRINT ENABLE/DISABLE 

COMMON /PRTFLG/JFLAG 


C ASSIGN JDOT DOTS 

DO 391 J-l.JDOT 
Kl-KOLD 

NNCK1 ) -NNCK1 )+l 

CALL GETDOTCMMCK1 ) .NNCK1 ) .NDARY, I DOT) 

C WRITECb, 99141 )K1. NNCK1 ). CNDARYCK),K«i. 15) 

■99141 FORMAT < ’ K1.NN«’,2I4. 4C/.2X.15I4) ) 

IFCK1.LE.M) GOTO 22? 

DO 220 Kl-l.M 
IFCNNCKD.GE.0) GOTO 225 
I F < ID0T.LT.MMCK1 ) ) GOTO 227 
I DOT- ID0T-MMCK1 ) 

£25 CONTINUE 

227 CALL CLMPXYCK1, IDOT.NL.NC.LL) 

CALL GTMPLBCNL.NC, LABEL) 

ISG-0 

CALL. LBLITPCLABEL.LB1. ISG) 

IFC ISG.EQ. 1) NSG- ’ YES’ 

IF 1 ' ISG. EO. 0) NSG-’ NO* 

IFC ISG.EQ. 1) NS(1)-NSC1)+1 
IFC ISG. EQ. 0) NS(2)»NSC2)+1 

I F C JFLAG . EQ . 1 ) WR I TE C 6 . 37 1 ) K 1 , I DOT . NL . NC , LABEL , LB 1 , NSG 
371 FQRMATC3X, I2.8X. I6.6X.2I4, iX, 14, 15. 4X.A4) 

391 CONTINUE 
C 

IFCNSC1) .GE.NSC2) ) MR-1 
IFCNSCn.LT.NSC2)) MR-2 
C WR1TEC6, 99301 )JDOT.NS, MR 
99301 KORUATC’ JDOT-’, 13 .’ NS"’ ,214, ’ 

C 

901 RETURN 

END 


MR-’, 13) 


13. A90: NON -BAYES IAN SEQUENTIAL DOT ALLOCATION, 

MAJORITY RULE LABELING 


13.1 DESCRIPTION 

This dot allocation and labeling scheme is different from most 
other schemes in one important respect: the allocation procedure 

is fused with the labeling procedure. 

For each cluster, dots are allocated continuously until the 
cluster is labeled or until the number of dots assigned exceeds 
a preselected value (currently, it is 35). The rule * «« non- 
Bayes ian sequential rule, stated as follows: 

a. Allocate INIDOT dots to the cluster. 1NIDOT was specified 
by the user during the program INPUT session. 

b. Allocate one dot, and compute the test interval (a, b) : 

x. 

a « -=• - 1.534 6, 
n i 1 

b - -i + 1.534 8< 
n i i 

e _/ X i (n i - *i> 

where n^ * number of dots allocated to cluster i 

x i ■ number of dots labeled as small grain for cluster i 
(notice < n^) 

c. If 0.5 e (a, b) and n^ does not exceed a preselected value 
(35) , repeat step b. 

d. Label the cluster by majority rule. 



//r 


It can be seen from tha abova that the number of dots assigned to 
the cluster may vary, depending on where the allocation and label- 
ing is stopped. Also, it should be noted that the labels above 
are either small grain or other. 


The equation for the proportion estimate is 


sg 


£ t 

cluster i 

labeled as 
small grain 


fr 


where N^ * number of pixels in cluster i 

N ■ total number of pixels in the entire scene 


(13-2) 


This scheme, like other allocation and labeling schemes, may be 
repeated with different pseudorandom sequences. The main program 
resets the random number generator to a different starting point 
for each repetition. 

All computation of statistics such as mean and variance is per- 
formed in the main program. A grand summary is printed. 


13.2 LINKAGES 

The PDP 11/45 Fortran library (for functions and subroutines such 
as SORT) , the PDP 11/45 system library (for routines such as the 
pseudorandom number generator) , and the utility package (described 
in section 14) are required. 

Main program 

or subprogram Subroutine or function required 
main saenbs, getdot, clmpcs, gtmplb, lblitp 

SAENBS CLMPCS , NBSMR 

getdot ran 



GETDOT, CLMPCS, GTMPLB, LBLITP 


NBSMR 

(non-Bayesian 
sequential 
majority rule) 

13.3 INTERFACE 

Interface with other routines is through the common block PRTGLG 
(2 bytea) , which is used to control the optional printing of dot 
files and other information. 

3.3.4 INPUTS 

Fortran-' ormatted input of the following parameters is needed: 
Card 


or 

line 

Parameters 

Format 

Default 

input 

Default 

value 

1 

Cluster map file name 

A13 

None 


2 

a. Number of repetitions 

13 

0 or blank 

1 


b. Starting point of first 
pseudorandom sequence 

15 

0 or blank 

10 


c. Number of repetition 
printings 

13 

0 or blank 

5 

J 

Number of initial dots 

13 

None 


4 

Number of status messages on 
terminal 

13 

0 or blank 

NO MSG 


Two disk files are needed, one containing the stripped cluster 
map, the other the stripped ground truth map. These files should 
have been created by conversion programs A81 and A86. 


tt? 


13.5 OUTPUTS 


In addition to the normal aoho of input specifications and status 
messages on the terminal, the following will he output on a line 
printer. 


Report 


Contents 


Individual repetition 1. Processor header 

2. Ground truth small-grain proportion 

3. A dot file showing dots chosen, their 
x-y position, their ground truth labels, 
and the test interval 


Grand summary 


4. A table showing number of clusters, 
cluster name codes, cluster sizes, 
number of dots assigned, and the split 

1. Processor header 

2. A table showing the estimate and bias 
of each repetition 

3. Bias, MSE, reduction in MSE, average 
estimate, variance, variance reduction, 
average and standard deviation of the 
total number of dots per repetition 


13.6 BRIEF ALGORITHM: FOR NON-BAYESTAN SEQUENTIAL DOT ALLOCATION. 

MAJORITY RULE LABELING 

1. Type message* at terminal and read from it run specifications. 

2. Compute ground truth small-grain proportion. 

3. Repeat a and b until all repetitions are finished s 

a. Set starting point of the pseudorandom number generator. 

b. For each cluster, allocate and label dots by non-Bayesian 
sequential majority rule. 



4 


. Compute bias, MSB, reduction in MSE, average estimate, vari- 
ance, variance reduction, average and standard deviation of 
the total number of dota per repetition, and print a grand 
summary. 






Op«i the ground truth date 
file end compute Its saell 
grain proportion 


Cell CIMKS to initialize 
cluster nap read 


Open the cluster i 



Call SAEhBS 

Coapute proportion estiaate 
(allocation: sequential) 
(labeling: majority rule) 


Store estimate a no 
number of dots used 


Update tpe seed to start 
pseudoranoen number 
generator 


Compute statistics of 
proportion estimate and 
number of dots 


Too 

many clusters? 


Check cluster sizes 
Group small clusters 
if necessary 


For each cluster 
call N8SMR to compute 
cluster small-grain 
estimate 


Print dot allocation 
information 




Print grand sumaary 


^Stop^ 



or more dots 
assigned? 


Assign and label a dot 


Figure 13-1.- Flow chart for non-Bayesian sequential 
majority rule labeling. 
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13.7 LISTING 

C — PROGRAM At? (A90.TSK) = PROPORTION ESTIMATE OF SMF1LL GRAIN 

C USING NON-BAYES I AN SEQUENTIAL DOT ALLOCATION 

C LABELLING IS NON-BAYES I AN SEQ MAJORITY RULE. 

C 

BYTE NAMEU5),NGT(13) 

REAL PSGC435) 

INTEGER NDOT(43S) 

C COMMON BLOCK IS A FLAG FOR PRINT ENABLE/DISABLE 

COfttON /PRTFLG/JFLAG 

DATA NGT(S)/’ . ’/,NGT(6)/’S’/.NGT(7)/’T’/,NGT(8)/’P’/ 

C 

c 

C READ I M3 INPUT FROM TERMINAL 

C 

WRITE(3,121) 

121 FORMAT ( * PROGRAMS fli? CA90.TSK) . ’ 

* /’ PROPORTIONAL ESTIMATION OF SMALL GRAIN’ 

* x' DOT ALLOCATION IS NON-BAYES I AN SEQUENTIAL 

* /’ LABELLING IS NON-BAYES I AN SEQ MAJORITY RULE.’ 

* /’ INPUT CLUSTER MAP FILENAME’/’ AAAAAAAAAAAAfi ' ) 

READC7. 122) (NAME(K) > K=1 • 13) 

122 FORMATC13A1) 

C 

WRITEC3, 131 ) 

131 FORMAT ( ’ AN ASSIGNMENT OF PSEUDORANDOM DOTS IS A REPETITION RUT 

* /’ HOW MANY REPETITION RUNS? SEED FOR FIRST PSEUDO SEQ?’ 

* /' HOW MANY PRINTINGS OF INDIUIDUAL RUNS’?’ 

* /’ III IIIII III’) 

READ (7,1 32 )JU, JSEED, JPAGE 

132 FORMAT (13; IX, IS, IX, 13) 

IF( JU.LE.0) JU=1 
IFCJSEED.LE.0) JSEED=10 
IF( JPAGE. LE.0) JPAGE =5 
JSK I P= ( JU-1 ) /JPAGE+1 

C 

WRITECS, 141) 

141 FORTtoT ( ’ NOTICE THAT NO. OF DOTS ASSGINED IN EACH’, 

* ’ REPETITION MAY UARY ’/ 

X ’ SPECIFY NO. OF DOTS INITIALLY ASSIGNED TO EACH CLUSTER.’/ 

X ’III’) 

READC7, 142)NNIX 

142 FORMAT (13) 

IF(NNIX.GE.l) CALL INIT(miX) 

C 

WRI TEC 0.151) 

151 FORMAT ( ’ STATUS MESSAGES ON TERMINAL, 1-04 MANY?’/’ Ill’) 

READ (7,152) MIES 

152 FORMAT (13) 

C 

WRITECS, 1S1 ) (NAME(K) , K=1 , 13) , JU, JSEED 
181 FORMAT (//’ CLUSTER FILE GIUEH = ’,13A1. 

X /’ NO. OF REPETITION RUNS FOR EACH TOTDOT -’,13, 

X /’ THE FIRST REPETITION RUN STARTS WITH SEED = ’ , 16) 

C 

C COMPUTE PROPORTIONAL ESTIMATE FOR GROUND TRUTH 

NLINE-117 
NC0L=196 
IX) 215 K=i ,4 
21S NGT(K) -NAMECK) 

0PENCUNIT=2, NAME-NGT, TYPE® ’OLD* , READONLY, FORM” ’ UNFORMATTED’ , 

X ACCESS® ’ DIRECT’ ) 

DO 231 L-l.NLINE 
DO 221 K = 1 , NCOL 
CALL GTMPLBCL.K, LABEL) 

CALL LBLITPCLABEL.LBl, IP) 

221 CONTINUE 

231 CONTINUE 

P= FLOAT ( I P ) /NL I NE/NCOL 
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C 

OPENCUNIT-l. NAME-NAME. TYPE- ’OLD’ . FORM- ’ UNFORMATTED ’ . 

X READONLY, ACCESS" ’DIRECT’ 1 
CALL CLMPLCCNLINE.NCOL) 

C 

c START ESTIMATION :‘OR EACH TOTDOT AND REPETITIONS 

JMES-0 
JS=JSEED 
DO 361 J-l.JU 
JFLAG-0 

IFCNODCJ-l.JSKIPI.EQ.^ tFLAS-1 
IF.CJFLAG.EQ.il WRITEC6.321 1NAME.P 

321 FORMAT ( ’IF , 10X. ’ ACCURACY ACCESSEMENT SOFTWARE C 5-29-79) ’ , /, 2X. 

* 'PROGRAM Al? CA90. TSKJ » ESTIMATION OF SMALL GRAINS PROPORTION’ 
X / . 1SX . ’ DOT allocation --non-bayesi an SEQUENTIAL .’ 

* /. 15X. ' DOT LABELL. I NG > NON-BAYES I AN SEQ MAJORITY RULE.’ 

X //.10X,’ INPUT CLUSTER MAP IS FILE ’.15A1, 

X /, 10X. ’ PROPORTIONAL ESTIMATE FOR GROUND TRUTH IS’.FB.S) 

IFCJFLAG.EQ.l) URITEC6.3221J, JSEED 

322 FORMATS. 2X, ’ REPETITION RUN- ’ , 13. 

X * RANDOM DOT SEQUENCE STARTS WITH’ > 18* ’ ’) 


CALL RANST ( JSEED 1 

CALL SAENBSCNDGTCJI.PSGCJI) 

JMES-JMES+1 

IFCJMES.LE.NMES1 WRITEC8.3231J, JSEED. NDOTCJI.PSGCJ) 

323 FORMAT ( ’ REPETITION RUN- 13. ’ SEED-’ , 16, ’ ND0T=’.I3, 

X ’ ESTIMATE-’ ,F9.5) 

JSEED-. TSEED+1S0 
361 CONTINUE 
C 

C PRINT GRAND SUMMARY FOR THIS JOB 

WRITE(6,321)NAME,P 
WRITEC6. 371 1 JU. JS 

371 FORMAT C/.14X. ’ GRAND SUMMARY OF THIS JOB 

X /, 14X, ’ NO. OF REPETITION RUNS =’.I3, 

X /, 14X, ’ RANDOM DOTS START WITH SEED-’, 16./ 

X , 10X. ’REPETITION PSEDUO SEQ DOT SMALL GRAIN BIAS’,/, 

X 10X, * RUN SEED ASSIGNED ESTIMATE WRTG.T. ’) 

C COMPUTE BIAS AND M.S.E. 

PB=0. 

PM-0. 


AM-0. 

UM-0. 

DO 381 J-l.JU 
TEMF-PSGC Jl-P 
PB-PB+TEMP 
PM- PM+TEMPXX2 
AM-AM+NDOTCJ) 

UM-UM+FLOAT C NDOT ( J) 1 XX2 
UR I TE ( 6 , 375 jJ. JS, NDOT ( J) , PSG ( J 1 , TEMP 
375 F0RMATC13X, I3.6X. I6,5X, 14, 5X, F8. 5, 3X. F9. 6) 

391 JS-JS+150 

FB-PB/.1U 
PM-PM/JU 
AM-AM/JU 
UM-UM/JU 

IFCJU.EQ.il UM-UM-AMXX2 

IFCJU.GT. 1) UM- C UM-AMXX2 1 X JU/ ( JU- i 1 

UM-SQRTcUMI 

AUERG-PB+P 

IFCJU.Ea.11 UAR<- PM-PBXX2 
IFCJU.GT.l) UAR- ( PM-PBXX2 1 X JU/ C JU- 1 ) 

C COMPUTE UARIANCE REDUCTION 

RR-PM/C PXC1.-P1/AM > 

REDUAR-UAR/C PX(1.-P)/AM 1 

WR I TE ( 6 . 395 1 PB , PM . RR , AUERS , UAR . REDUAR , AM , UM 
38S FORMAT C / ’ BIAS- ’ , F10. 6, ’ M.S.E. , F10.6, ’ REDUCTION- ’, F10 6, 
X /’ AUERAGE-’ -F10.6, ’ UARIANCE- ’ . F10.6, ’ REDUCTION-’ , FIO. 6. 

X /’ AUE DOT-’ .F10.3. ’ DOT S. D. - ’ . F10. 31 


A2J 





C ) 

C ■ J 

WRITEC6.401) 

401 FORNATC’l END OF THIS JOB * ) 

C 

STOP 

END 

C | 3 

C 

SUBROUTINE SAENBSCNDOT. PSG) 

C STRATIFIED AREAL ESTIMATION USING NON-BAYES I AN 

C SEQUENTIAL DOT ALLOCATION WITH MAJORITY RILE LABELLING. 

C ON RETURN, THE PROPORTIONAL ESTIMATE IS IN PSG, AND THE 

C TOTAL NO. OF DOTS ASSIGNED IS IN NDOT. 

C 

C UR I TEN AND EDITED BY N.Y. D-D ON 5-29-79. 

C 

c 

INTEGER NDARY(435).mi51),m(51),NX(51).LL(51) 

INTEGER NNS<2.51) 

INTEGERS NSG 

C COMMON BLOCK IS A FLAG FOR PRINT ENABLE/D I SABLE 

COMMON /PRTFLG/JFLAG 
DATA NST ART/0/, MINSZ/35/ 

C 

C GET PIXEL COUNT FOR EACH CLUSTER 

I F ( NSTART . NE . 0 ) GOTO 121 
NSTART = 1 

CALL CLMPCS ( NP I XEL , M , MI i , LL ) 

IFCM.6E.51) GOTO 901 
121 CONTINUE 
C 

C MAKE SURE EACH CLUSTER HAS A MINIMUM NO. OF PIXELS 

Kl-0 

DO 151 K-l.M 

NN(K>»0 

NX(K)-0 

I F ( MM ( K ' .LT.MINSZ) K1-K1+MMCIO 
IF(MMCK).LT.MINSZ) NN(K)=-1 
151 CONTINUE 
MX-M 

IFCK1.LT. MINSZ) GOTO 161 
MX-M+1 
MMCMX) =K1 
161 CONT I NUE 
C 
C 

C BEGIN NON-BAYES I AN SEQ MAJORITY RULE DOT ALLOCATION 

C 

C 

IFCJFl AG.EQ. 1) WRITECS, Oil. ) 

211 FORMAT <//’ DOT FILE ’ 

* /’ CLUSTER DOT NUMBER POSITION GROUND TRUTH 

* , ‘ TEST INTERNAL’ 

* /’ NO. MRT ITS CLUSTER LINE. COL RAH, CODE’ 

* ’ LOWER UPPER’ ) 

C 

C ASSIGN DOTS AND COMPUTE ESTIMATE 

NDOT-0 

P5G=0. 

ND-1 

DO 231 K-l.HX 

IFi NNCK) .LE.-l) GOTO 231 

I F C NDOT . GT . C -35+435 ) ) GOTO 251 

CALL NBSMRC ISG.NNS(l.K), NDARYCNP) . K. MM. NN, LL, NX, M) 

PSG-PSG+FLOAT CISC) *MM ( K > /NP I XEL 
NDOT « NDOT +NN C K ) 

ND=ND+NN(K) 

231 CONTINUE 
GOTO 255 
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c 

c 

251 

253 

255 

C 

C 

C 

c 

282 


283 

234 


28 S 


288 

C 

351 

C 

C 

301 

902 

C 

C 

930 

C 

C 

c 

c 

c 

c 

c 

C 

c 

c 

c 

c 

c 


— TERMINATE DOT ASSIGNING IF ARRAY OVERFLOW IS ANTICIPATED 
I F ( JFLAG . EQ . 1 ) WRITE(6> 253)NDCT 

FORMAT ( ’ DOTS ASSIGNED* ’, 14, ’ MAY CAUSE ARRAY QUERFLOW, ’ , 

* ’ ALLOCATION TERMINATED. ) 

CONTINUE 

— END OF NON-BAYES I AN SEQ MAJORITY RULE DOT ALLOCATION 

-- WRITE CLUSTER INFORMAT I OS 

IF( JFLAG.EQ. 1 ) WRITE(6,282)NMT,NPIXE1_,MX 

FORMAT'//, SX, ’ TOTAL NO. OF DOTS ASSIGNED , IS, 

* / , 5X > ’ TOTAL NO. OF PIXEL IN THE CLUSTER MAP =’,I£, 

* /.SX, ’ TOTAL NO. OF CLUSTERS -’,16, 

* /.10X,’ CLUSTER CLUSTER I'd. OF NO. OF MTS' 

* ,',10X.’ NO. CODE PIXELS ASSIGNED SPLIT’) 

DO 283 K-l.MX 

IF(NN(iO.LE.-l) NNOO-0 
IF(NNSCl.K) .GE.NNS(2,K) ) NSG-’S.G.’ 

I F ( NNS ( 1 , K ) . LT . NNS ( 2 , K ) ) NSG-’ 

IF(JFLAG.EQ.l) URITEC6. 284)K,LL(K) ,MMCE) ,NN(K) , NNS(l.K) ,NNS(2.K) 
t , NS3 

CONTINUE 

F0RMATU3X, 12, 5X, 14, 3X, 16. 4X, 13, 7X. 12. ’ , 12. IX, A4) 

IF(MX.EQ.M) GOTO 288 
KX=0 

DO 28S K-] ,M 

IFCNN(K) .NE.0) GOTO 28S 

KX-KX+1 

NDARYtKX) 

CONTINUE 

IFCJFLAG.PZQ.l) WRITE (6, 286) MX, (NDARY(K) , K«1 , KX) 

FORMAT </’L THOSE SMALLER CLUSTERS GROUPED (AS CLUSTER* , 13, 

* ’) ARE:’, 3(2,16(13, ’,’)) ) 

CONTINUE 

I F ( JR.AG . EQ . 1 ) WR 1 TE ( 6 , 35 1 ) NDOT . PSG 

FORMAT ( ’ TOTAL DOTS ASSIGNED* ’ > 14, ’ PROPORTIONAL’, 

* ’ ESTIMATE-’ ,F10. 6) 


GUTO 99© 

WR I TE C S . 302 ) M 

FORMAT ( ’ YOU HAUE’ , 13, ’ TOO MANY CLUSTERS (MAX-50) ’ ) 
GOTO 990 


RETURN 

END 


SUBROUTINE NBSMR( ISG.NS. ISDARY , K 1 , MM > NN > LL , NX > M ) 

— NON-BAYES I AN SEQ MAJORITY RULE 

— THIS ROUTINE ASSIGN A DOT OR MORE TO 

-- CLUSTER K1 ACCORDING TO A NON- BAYES I AN SEQ MAJORITY RULE. 

-- ON RETURN, NN(K1) CONTAINS THE NO. OF DOTS ASSIGNED 
DURING THIS CALL. 

— ISG-1, ON RETURN, MEANS SMALL GRAIN LABEL IS ASSIGNED TO 

— THE CURRENT (Kl) CLUSTER, IF ISG-0, NOT SMALL GRAIN 

— NDARY = DOT ARRAY FOR THE CURRENT CLUSTER 

— MM- ARRAY CONTAINING CLUSTER SIZES 

— NN- ARRAY CONTAINING CLUSTERS' NO. OF ASSIGNED DOTS 

— NX- ARRAY CONTAINING CLUSTERS’ NO. OF SMALL GRAIN DOTS 

— LL- ARRAY CONTAINING CLUSTERS’ CLASS NAME, 

JUST FOR USE IN SUB. CLMPXY 

— M- NO. OF CLUSTER (NOT COUNTING THE GROUPED CLUSTER, IF EXISTS; 



c 

INTEGER NDARY<1).MM<1).NN<1>.NX(1),LL<1) 

INTEGER NSC2) 

INTEGER** NSG 

C COMMON BLOCK IS A FLfiG FOR PRINT ENftBLE/D I SftBLE 

COmON /PRTFLG/JFLAG 
DftTft IN I DOT/ 1/ 

C 

C 

C ASSIGN ft DOT 

NSCD-0 

NS<2)=0 

KOLD-K1 

KCOUNT-0 

121 IF((KC0UNT.EQ.35).AND. C JFLAG. EQ. 1 ) ) WRITEC6, 122)K1 
I F ( KCOUNT . EQ . 35 ) GOTO 901 

122 FORMAT ( ’ FOR CLUSTER’ , 13. ’ DOT ASSIGNMENT TERMINATED’, 

* ' (ONLY 35 DOTS ALLOWED) . ’ ) 

KCOUNT-KCOUNT+1 

NNCK1)«NNCK1)+1 

CALL GETD0T(MMCK1).NN(K1),NDARY, IDOT) 

C WRITEC6, 99141 )Ki. NNCK1 ). CNDAFYCK) , K=1 . 15) 

99141 FORMAT (’ K1.NN-* .214. 4(/,2X.15I4) ) 

IF(Kl.LE.M) GOTO 227 
DO 225 K1=1.M 
IFCNNCKD.GE.0) GOTO 22S 
I FC IDOT.LT.MMCK1) > GOTO 227 
IDOT-IDOT-MMCK1) 

■ 325 COST I NUE 

227 CALI. CLMPXY ( K1 » I DOT , NL , NC . LL ) 

CALL GTMPLB(NL,NC. LABEL) 

Kl-KOLD 

ISG-0 

CALL LBLITPCLABEL.LB1, ISG) 

IF(ISG.EQ.l) NX<K1)-NXCK1)+1 
IFUSG.EQ.l) NSG-’S.G. ’ 

IF( ISG. EQ. 0) NSG- ’ 

I F C ISG. EQ. 1 ) NSC 1 ) -NSC 1 ’+1 
IFCISG.EQ.0) NSC2)-N3C2)+1 
C 

C NON-BAYES I AN SEQ CLUSTER LABELLING 

TX-NXCK1) 

TN-NN(Kl) 

IFCNNCKD.GT.l) SHEAD-SQRTC TX*<TN-TX)/( TN*TN*CTN-1. ) ) ) 
TL-TX/TN-1 . 534*SHEAD 
TU-TX/TN+1 . S34*SHEAD 
IF(NN(Kl).EQ.l) TL-0. 

IFCNNCKD.EQ.l) TU-0. 

IFC JFLAG. EQ. 1) WRITECS, 371)K1, I DOT . NL , NL , LABEL , LB1. NSG. TL.TU 
371 F0RMATC3X, 12, 8X. I6.6X.2I4, IX. 14, 15, IX, A4, 2F9.4) 

C CONT’NUE DOT ASSIGNMENT UNTIL ALL INITIAL DOTS ARE DONE 

IFCNNsKD.LT. INIDOT) GOTO 121 

C CONTINUE DOT ASSIGNMENT UNTIL THE CRITERIA IS NOT SATISFIED 

IFC C0.5.GT.TL) .AND. (0.5.LT.TU) ) GOTO 121 

C NOW TERN I NATE DOT ASSIGNMENT AND LABEL BY MAJORITY RULE 

IF(NS(1).GE.NS(2)> ISG-1 
IFCNSC 1 ) . l.T. NSC2) ) ISG-0 
C 

901 RETURN 
C 

C ENTRY TO SET INITIAL NO. OF DOTS ASSIGNED 

ENTRY INITCNNIT) 

I N I DOT -MS I T 
RETURN 

C END OF THIS ENTRY 

C 

END 
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14. A06: THE UTILITY PACKAGE 


The utility package consists of the subroutines CLMPCS, GETDOT, 
GTMPLB , MR, and LBLITP. Subroutine CLMPCS and its entries CLMPXY 
and CLMPLC read the cluster map from its disk file. Subroutine 
GETDOT chooses pixels from clusters and assigns them as dots. 
Subroutine GTMPLB reads the ground truth map from its disk file. 
Subroutine MR finds a majority label for a group of labels. Sub- 
routine LBLITP interprets a given label as small grain or other 
and increments a counter if the label is small grain. Each sub- 
routine is discussed in one of the subsections following. The 
discussion includes the purpose and entry points of the subroutine 
its linkages and interfaces, its inputs and outputs, the storage 
requirements, and a description of the operation of the subroutine 
The listing for the entire utility subroutines package is given 
as the last subsection. 

14 . 1 CLMPCS 

The subroutine CLMPCS and its entries read the cluster map. This 
subroutine has three entry points: 

CLMPCS — to read the cluster map and return cluster sizes and 
numbers 

CLMPXY - to return the x,y position of the jth pixel of the ith 
cluster 

CLMPLC — to set the size of the input cluster map 

14.1.1 LINKAGES 

This routine does not call any other subprogram. 

14.1.2 INTERFACES 

The CLMPCS subroutine interfaces with other routines through the 
calling arguments. 
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14.1.3 INPUTS 

Input to the subroutine consists of the cluster nap, named as 
*.STP, the file output by the A86 processor. 

Calling sequences CALL CLMPCS(NPIXEL,M,MM,LL) 


Parameter 

Type 

Dimension 

In/Out 

Definition 

NPIXEL 

1*2 

1 

Out 

Total number of pixels in 
this map 

M 

1*2 

1 

Out 

Total number of clusters 

MM 

1*2 

Variable 

out 

Array containing sizes of 
all clusters 

LL 

1*2 

Variable 

Out 

Array containing the clusters' 
numbers (the first cluster 
is 1, second 2, etc.) 

Entry calling sequences CALL CLMPXY (I, J,NL,NC,LL) 

Parameter 

Type 

Dimension 

In/ Out 

Definition 

I 

1*2 

1 

In 

The cluster ith in consideration 

J 

1*2 

1 

In 

The jth pixel of the ith 
cluster 

NL 

1*2 

1 

Out 

Line number of that pixel 

NC 

1*2 

1 

Out 

Column number of that pixel 

LL 

1*2 

Variable 

In 

Array containing the clusters' 
numbers 

Entry calling sequences CALL CLMPLC (LL1,LC1) 

Parameter 

Type 

Dimension 

In/ Out 

Definition 

LL1 

1*2 

1 

In 

Total line number of the 


cluster map 




LC1 1*2 1 In Total column number of the 

clutter map 

14.1.4 OUTPUTS 

Th3 results are returned for use by the calling routine. 

14.1.5 STORAGE REQUIREMENTS 

This subroutine requires a fair amount of storage. 

14.1.6 DESCRIPTION 

When CLMPCS Is called, It reads, line by line, the entire cluster 
map and determines the number of clusters and their sizes. Next, 
it sets up an internal table that indicates the number of pixels 
of cluster i in line 1, for all clusters i and lines l . This 
table will later be referenced to determine the x,y position of 
a pixel of a cluster. 

When CLMPXY is called, it searches through the internal table that 
was set up when CLMPCS was called. The table immediately gives 
the line number of pixel j of cluster i. Next CLMPXY reads in 
the line from the cluster map, and, by pixel-to-pixel comparison, 
it gives the column number of that pixel. 

When CLMPLC is called, the cluster map size is set to that speci- 
fied by the calling arguments, even when the actual cluster map 
is larger than that indicated by those arguments. 

14.2 GETDOT 

The subroutine GETDOT selects pseudorandomly with uniform proba- 
bility a pixel from a group of pixels (cluster) . This subroutine 
has two entry points t 

GETDOT - to get a dot assigned 
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RANST — to skip sons pseudorandom numbers so that the pseudoran- 
dom sequence will be different 

14.2.1 LINKAGES 

This routine calls the RAN subprogram. 

14.2.2 INTERFACES 

The GETDOT subroutine interfaces with other routines through the 
calling arguments. 

14.2.3 INPUTS 

Calling sequence: CALL GETDOT (NP,L,NDARY, IX) 


Parameter 

Type 

Dimension 

In/ Out 

Definition 

NP 

1*2 

1 

In 

Total number of dots in the 
cluster of interest 

L 

1*2 

1 

In 

Lth dot is requested (i.e., 
L - 1 dots exist in the dot 
airay) 

NDARY 

1*2 

Variable 

In/ Out 

Dot array (all L - 1 dots are 
in ascending order) 

IX 

1*2 

1 

Out 

The selected pixel (now 
called dot) 

Entry calling sequence: CALL RANST (ISTART) 

Parameter 

Type 

Dimension 

In/Out 

Definition 

ISTART 

1*2 

1 

In 

Number of pseudorandom num- 


bers to be skipped 


14.2.4 OUTPUTS 

The results are returned for use by the calling routine. 



14.2.5 STORAGE REQUIREMENTS 

This subroutine requires small storage. 

14.2.6 DESCRIPTION 

Depending on L, the subprogram GETDOT computes the number of free 
pixels available; it is (NP - L + 1). Using the system pseudoran- 
dom number generator (RAN) , a number between 0 and 1 is obtained 
and then scaled to between 1 and (NP • L + 1). If there are no 
previously selected dots (i.e., L ■ 1), then the scaled number is 
the pixel number of the selected dot. However, if L > 1, then 
the subroutine needs to adjust the scaled number by referring to 
dots in the dot array. Since the dots are present in ascending 
order, the insertion of the new dot can be done in an orderly 
fashion, starting with the first dot in the dot array. 

14.3 GTMPLB 

The subroutine GTMPLB returns the label of a pixel (given its 
x,y position) by reading the ground truth map. 

14.3.1 LINKAGES 

This routine does not call any other subprogram. 

14.3.2 INTERFACES 

The GTMPLB subroutine interfaces with other routines through the 
calling arguments. 

14.3.3 INPUTS 

Input to the subroutine consists of the ground truth map, named 
as *.STP, the file output by the A81 processor. 




Call ing sequence ! CALL QTMPLB (NL,NC, LABEL) 


Parameter 

Type 

Dimension 

In/Out 

Definition 

NL 

1*2 

1 

In 

Line number of the pixel 

NC 

1*2 

1 

In 

Column number of the pixel 

LABEL 

Byte 

1 

Out 

The label of the pixel 

14.3.4 OUTPUTS 





The results are returned for use by the calling routine. 

14.3.5 STORAGE REQUIREMENTS 

This subroutine requires small storage. 

14.3.6 DESCRIPTION 

Since the ground truth map stripped by processor A81 is directly 
accessible by Fortran , all that is needed to read in a line is 
just a simple READ statement. After a line has been read into 
a buffer, the label of the pixel can easily be picked up. 

14.4 MR 

The subroutine MR returns a majority rule label for a group of 
labels. 


14.4.1 LINKAGES 

This routine does not call any other subprogram. 

14.4.2 INTERFACES 

The MR subroutine interfaces with other routines through the 
calling arguments. 





14.4.3 INPUTS 

Calling sequences CALL MR (XG,NG, LABEL, Jl,JG,LGj 


Parameter 

ZZES. 

Dimension 

In/Out 

Definition 

IG 

Byte 

Variable 

In 

Array containing the group 
of labels 

NG 

1*2 

1 

In 

Dimension of IG 

LABEL 

Byte 

1 

Out 

The found majority rule label 

J1 

1*2 

1 

Out 

Number of different labels 
in IG 

JG 

Byte 

Variable 

Out 

Array storing those different 
labels (dimensioned as Jl) 

LG 

1*2 

Variable 

Out 

Numbers of the labels in JG 
(dimensioned as Jl) 

14.4.4 OUTPUTS 





The results are returned for use by the calling routine. 

14.4.5 STORAGE REQUIREMENTS 

This subroutine requires small storage. 

14.4.6 DESCRIPTION 

The subroutine MR scans the input a::r«-.y IG in ascending order, 
stores the different labels in JG, counts the number of each label 
and stores these counts in LG, and selects the label with the 
greatest count. If two lcbels have that count, the one appearing 
first will be selected. 

14.5 LBLITP 

The subroutine LBLITP interprets a raw ground truth label and 
increments a counter if the label is a small grain. 


14.5.1 LINKAGES 


This routine does not call any other subprogram. 


14.5.2 INTERFACES 

The LBLITP subroutine interfaces with other routines through the 
calling arguments. 


14.5.3 INPUTS 


Ca< ing sequence: CALL LBLITP (LABEL, LBl, IP) 

Parameter Type Dimension In/ Out Definition 

LABEL Byte 1 In The input raw ground truth 


LBl 

IP 


1*2 

1*2 


label 

1 Out The crop code of the label 

1 In/Out The counter to be incremented 

if the label is a small grain 


14.5.4 OUTPUTS 

The results are returned for use by the calling routine. 


14.5.5 STORAGE REQUIREMENTS 

This subroutine requires small storage. 


14.5.6 DESCRIPTION 

The crop code of the label is computed as LBl = LABEL + 128. The 
counter IP is incremented if any one of the following is true: 

a. 1 < LBl < 62 

b. 99 < LBl < 104 

c. LBl = 109 

d. 124 < LBl < 129 



e 


LBl = 134 


14.6 LISTING 


C IN PROGRAM A06.FTN* THE UTILITY PACKAGE 

C THIS CONTAINS ALL ESSENTIAL SUBROUTINES FOR AA SOFTWARE 

C- PREPARED DURING SPRING. 1979. FOR R.T.E USUAGE. 

C 

C DIRECT-ACCESS U.E. RANDOM ACCESS) TO CLUSTER MAP 

C — AND GROUND TRUTH MAP IS EMPLOYED. 

C THEY SHOULD BE THE STRIPPED UERSIONS CREATED BY 

C A DIFFERENT TASK. N CHU 5-3-79 

C 

C 


C 

SUBROUTINE CLMPCSCNPIXEL.M.MM.LL) 

C ROUTINES TO READ CLUSTER MAP 

C 

C ENTRY* 

C CLMPCS READ CLUSTER MAP AND RETURN CLUSTER SIZES AND NO. 

C CLMPXY GET JTH PIXEL OF I TH CLUSTER, RETURN X.Y POSITION 

C CLMPLC SET SIZE OF CLUSTER MAP, LINE & CCS. DIMENSION 

C 

C 

LOGICAL*! JLC5-4C, JD(468).J1 
INTEGER*2 J2. 1 1 (256) . 1 1 K2S6.S1) 

INTEGER MN< 1 ) >LL( 1 ) 

EQUIUAIENCE (JLC73). JD). <J1, J2) 

DATA LINEAUY0Y, NLINEY117Y.NC0L/196X 
C 


C ENTER I Na uLMPCS* MND CLUSTER SIZES AND Nu. 

C 

C 

C CLEAR THE ARRAYS 

DO 1-41 K* 1.256 
141 1 1 <K)=0 

DO 161 L- 1 . NLINE 
DO 1S1 K=l.St 
1S1 I I1(L>K) S 0 
161 CONTINUE 
C 

C CHECK FOR CLUSTERS 


Kl=NC0L+?2 
DO 221 L-l. NLINE 
READCl’L)<JD<K>,K«i,NCOL) 
IX) 211 K-i.NCOL 
J1=JD(K) 


C URITE(6,99002)JP 
99002 FORMAT ( ’ J2»’ . 16) 

211 I I ( J2) = I I C J2)+l 
C WRITEC6.99005)L, ( JD(K) . K-l.NCOL) 

99005 FORMAT (’ LINE= ’ , 13. /. 10C ’ ’.1614,/) ) 
221 CONTINUE 
C 

C WRITEC6, 99101) 1 1 
99101 FORMAT C ’ II = ’.13IS) 

C COMPACT CLUSTER SIZES IN ARRAY MM 


M-0 

NPIXEL=0 

DO 291 K= 1,256 

IF( 1 1 < K ) . EQ . 0 ) GOTO 291 

M-M+l 

MMCM) = 1 1 (K) 

LL(M) *K 
1 1 CK) =M 

C WRITEC6, 99104)K, M. LL(M) 

99104 FORMATS 99104 K,M,L.L=’ ,515) 

NP I XEL *NP I XEL+MM ( M ) 

291 CONTINUE 

L l*!R I TEC 6, 99121 ) CMMCK),K=1»M) 

99121 FORMAT C ’ MM=’,12I5) 





C URITE<6.99191)(LL(K). K-l.M) 

99131 FORMAT'.' LL-’.12I5) 

C 

C STORE PIXELNO FOR EACH CLUSTER ON EACH LINE 

DO 321 L’l.NLINE 

READ ( 1 * L H JD < K ) . K » 1 . NCOL ) 

C WRITE<6, 99135) (JDCK) .K-l.NCOL) 

99135 FORMATC AF*’,1SI4) 

DO 311 K-l.NCOL 
Jl-JD(K) 

MX- I I(J2) 

IIKL.MX)-IIi<L.MX)+l 
311 CONTINUE 

IF(L.EQ.l) GOTO 321 
DO 315 K-l.M 

315 II1(L.K)*II1(L»K)+I11( (L-l ) . K) 

C MRITE(6. 99141 ) ( I ll(L.K) .K-l.M) 

99141 FORMAT ( ’ I II- *,1215) 

321 CONTINUE 
C 

GOTO 990 
C 

c 

c END OF CLMPCS 

C - 
r 

C ENTERING CLMPXY< FIND LINE AND COL NO. OF PIXEL J IN CLUSTER I 

C 

ENTRY CLMPXYU.J.NL.NC.LL) 


C LOOK THROUGH TABLE II1CL.K) TO FIND LINE NO. 

DO 411 h#_=l >NLINE 
IFilUCNL, D.GE.J) GOTO 451 
411 CONTINUE 
451 CONT INLIE 

C 

C DIRECT-PCCESS TO THE CLUSTER MAP IN *.STP FILE 

IFfNL.EQ.LINEAU) GOTO 461 
LINEAU-NL 

READC l’NLHJD(K).K-l. NCOL) 

461 CONTINUE 
C 

C TEST THE LINE BY COUNTING TO THE COLUMN NO. 

4S1 CONTINUE 

IF(NL.EQ.l) JX-J 
IF(NL.NE.l) JX-J- I I1(NL-1» I ) 

MX-LL.C I ) 

C UR I TE ( 6 , 99741 ) JX . MX 

99741 FORMATS 99741**** JX,MX»’,SI5> 

DO 491 NC-l.NCOL 
IF( JD(NC) . EQ.MX) JX-JX-1 
IFC JX.EQ.0) GOTO 495 
491 CONT I IDE 

495 CONTINUE 

C WRIT£<6. 99751 )I,J.NL.NC 

99751 FORMAT ( ’ 99751*** I . J.NL.NC-’ ,615) 

GOTO 990 
C 

C ENTERING CLMPLC* TO SET SIZE OF CLUSTER MAP IN CONSIDERATION 

C 

ENTRY CLMPLC (LL1, LCl) 

NLINE-LL1 
NC0L-LC1 
GOTO 990 
C 

C END OF CLMPLC 

C 

c 

990 RETURN 
END 


1 


c 

c 

c 

SUBROUTINE GETDQT ( NP » L . NDARY . IX) 

C FROM ft POOL OF (NP-L+1) PIXELS. THIS ROUTINE 

C PROVIDES ft DOT BASED ON UNIFORM PROBABILITY DISTRIBUTION. 

C THE DOT HILL BE CORRECTLY INSERTED AMONG PREUIQUS GENERATED 

C DOTS. 

C NP- TOTAL NO. OF PIXEL IN THE CLUSTER OF INTEREST 

C L- L TH DOT IS REQUESTED 

C NDARY IS THE ARRAY CONTAINING ALL (L-l) DOTS t ASCENDING ORDER) 

C IX- THE PRESENT ASSIGNED DOT’S PIXEL NO. W.R.T. NP PIXELS 


C 

INTEGER NDARYC1) 
DATA Jl/0/. J2/0/ 
C 


C GENERATE RANDOMLY THE DOT’S PIXEL NUMBER 

C IX- THE NO W.R.T. (NP-L+1) PIXELS 

X-RANCJ1. J2) 

IX-X* (NP-L+1 )+i 
C WRITE(S.99201)L.X, IX 

99201 FORMAT ( ’ IN GETDOT. L.X. IX-M4.F7.3. 14) 

C 

C CHECK IF NO DOT PREVIOUSLY ASSIGNED 


IF(L.GT.l) GOTO 131 
NDftRY(L) * IX 
GOTO 901 

C 

C FIND THE DOT'S CORRECT PIXEL NO W.R.T. TO NP PIXELS 

131 Kl-L-1 

DO 161 K-1.K1 
IF( IX-NDARY(K) )181. 151. 1S1 
1S1 IX-IX+l 

C WRITEC6.9930DK1.K2. IX.NDARYCK) 

99301 FORMAT ( ’ K1 , K2, IX. NDARYCK) - ’ . 614) 

161 CONTINUE 
K-L 


C INSERT THE NEW EOT INTO THE ARRAY 

161 IF( (Kl+1 ) .LE.K) GOTO 191 
NDARY ( Kl+1 ) -NDARY (K1 ) 

Kl-Kl-1 

C WRITE(6. 99401 )K1.NDARY(K1).K 
99-YU FORMAT ( ’ 99401, Kl. NDARY (Kl) .K=’ .414) 
GOTO 181 

191 NDARYCK) -IX 
GOTO 901 


C 

C 

C— 
c — 

C 


401 

901 


ENTRY RANSTC I START) QJ 

— SETTING STARTING OF RANDOM NUMBER GENERATOR^ 

— INSTART MUST BE GREATER THAN ONE 

Jl-0 





J2-0 

DO 401 K-l. I START 
X-RANUJ1. J2) 
RETURN 
END 




c 

c 

SUBROUTINE GTMPLB ( NL . NC • LABEL ) 

c T0 find THE LABEL ON THE COMPRESSED GROUND TRUTH (IN *.STP) 

c THAT LOCATES AT NL LINE AND NC COL. 

C 

C 

BYTE J1 ( 360 ).JX, LABEL 
DATA L I NEAU/0/ , NCOL/ 1 96/ 

C 


/37 


C DIRECT-ACCESS TO GROUND TRUTH DISK FILE 

IFCNL.EQ.LINEAU) GOTO 201 
LINEAU-NL 

READ C 2 ’ NL) ( J 1 ( K > . K ■ 1 . NCOL ) 

201 LABEL- Jl(NC) 

990 RETURN 
END 
C 
C 

SUBROUTINE MR( IG.NG.LABEL, Jl, JG.LG) 

C IG, WHEN ENTERS. CONTAIN NG LABELS 

C ON RETURN. LABEL WILL CONTAIN THE FIRST LABEL 

C BY MAJORITY RULE 

C ALSO. ON RETURN. Jl WILL CONTAIN NO. OF DIFFERENT LABELS FOUND 

C JG CONTAIN THOSE LABELS, AND LG CONTAIN THE NO. OF THESE LABELS 

C 

BYTE I G C 1 J , JGC 1 ) .LABEL 
INTEGER LG( 1 ) 

C 

Jl-1 

JG(l)-IGCJl) 

LGCJD-l 

LABEL-IG(Jl) 

IFCNG.LE.l) GOTO 901 
C 

C COUNT SUBPIXELS HAUEING SAME LABELS 

DO 251 I-2.NG 
DO 221 J-l.Jl 

C URITE(6. 99108)1. J. IG( I) . JGCJ) . Jl 

99108 FORMAT (’ I , J, IG. JG. Jl- ’ . 1014) 

IFCIGCI).EQ.JG(J)) GOTO 231 
221 CONTINUE 
Jl-Jl+1 
JG( Jl )• IG( I ) 

LG ( J 1 ) * 1 
GOTO 251 

231 LG(J)»LG(J)+1 

251 COST I MJE 

C 

C CHOOSE THE MAJORITY LABEL 

MAX-0 

DO 281 J-l.Jl 

IFCLGC D.LE.MAX) GOTO 281 

MAX-LG(I) 

J-I 

281 CONTINUE 

LABEL-JG(J) 

901 RETURN 
END 
C 

C , ' 

SUBROUTINE LBLITPfLABEL.LBl. IP) 

C INTERPRETATE A RAW LABEL (ONE BYTE LONG) AND RETURN 

C CODE IN I.B1 (INTEGER). IP IS INCREMENTED BY ONE IF 

C- — THE LABEL SHOWS SMALL GRAIN. 

C 

BYTE LABEL 
LB1 -LABEL 
LB1-LB1+128 

IF( (LB1 . GE. 1 ) .AND. (LB1 . LE. 82) ) GOTO 401 
IF( (LB1.GE.99) .AND. (LB1.LE. 104) ) GOTO 401 
IFd.Bl . EQ. 109) GOTO 401 

IF( (LB1.GE. 124) .AND. (LB1.LE. 129) ) GOTO 401 
IFCLB1.EQ.134) GOTO 401 
GOTO 901 
401 IP- I P+1 
901 RETURN 
END 


15. EXAMPLE OF OPERATING PROCEDURE 


All programs are written for an interactive environment but execu- 
tion in batch is also possible. When a program begins execution , 
guide messages will appear on the terminal and inputs are expected 
to be keyed in. Most messages contain a format specification: 

"A" indicates alphanumeric characters are expected, "I" an integer 
expected, and N F N a floating point number. The format field is 
indicated by the length of a whain of characters on the terminal. 
Those keyed-in characters not underneath the format field will be 
ignored. To run any program, just type, after MCR>, "RUN (the 
program's name) " and hit the ESC key. 


Suppose an analyst desires to test a cluster map 100519101. DT2 
against the Bayesian dot allocation scheme with adaptive prior. 
Suppose the ground truth map is initially in 100577278. GT0. The 
analyst needs to execute three jobs. 

1. To strip the cluster map 

MCRsRUN A86S 

PPDGPAM All <t.A8b.TSK> . TO STRIP A CLUSTER MAR FILE. 

INPUT' CLUSTER MAP FILE NAME. 

AAAAAhAAAAAF.A 
1 0 05 1 A 1 i.i 1 . DT c‘ 


£0 

LINES 

WRITTEN 

UN 

1 i.i 0 

'MA1U1 . S 

TP 


40 

LINES 

WRITTEN 

□N 

10U 

1519101. S 

TP 


P M 

LINES 

WRITTEN 

ON 

1 1.1 IJ 

ib 191 Ul. S 

TP 


S 0 

LINES 

WRITTEN 

LIN 

1 0 U 

ib 191 i.i] . :s 

■ TP 


1 0 0 

LINES 

WRITTEN 

UN 

1 U U 

15191 ut.S 

TP 


STRIPPED MAP OF SIZE 

117 

LINES BY 

1 

. CLlCS; PRODUCE U 

NEST 

CLUSTER MAP PL 

..EA: 

SE . 

TO STOP* 

■JUS 

T TYPE STOP 

PROOF 

'’An All 

< ABC . T ST 

■ > . 

TD 

STRIP A 

CLUS 

TER MAP FILE. 


INPUT CLUSTER MAP FILE NAME. 

ahaahhahahaaa 

STOP 

ABC — STOP 


2. To strip a ground truth nap 

MCR> RUN RSI * 

PROGRAM : R1 0 < A8 1. TSfc> . TO REDUCE GROUND TRUTH 
RESOLUTION. OUTPUT IN STRIPPED PORMHT 
TYPE GROUND TRUTH FILE. NAME 
RAAAAAAAAARAA 
1 00577c!? 8. GTO 

GROUND TRUTH FILE FROM 100577S78.GT0 

LABEL BY MAJORITY RULE. . . EXECUT ION BEGINS. . .UNIT 


£0 

LINES 

WRITTEN 

ON 

1005. ST P 

4 0 

LINES 

WRITTEN 

ON 

1005.STP 

b 0 

LINES 

WRITTEN 

ON 

t 005. S TP 

S'.lj 

LINES 

WR 1 1 TEN 

ON 

1005. ST P 


100 LINES WRITTEN ON 1005.STF' 

STRIPPED MRP OF SIZE 11? LINES BY 19b COLS PRODUCED IN 1005. ST 

NEXT GROUND TRUTH MRP PLERSE . TO STOP* JUST TYPE STOP 
PROGRAM s A U' < HS 1 . TSK > . TO REDUCE GROUND T PUT H 
RESOLUTION. OUTPUT IN STRIPPED FDPMRT 
TYPE GROUND TRUTH FILE NAME 

aaaaaaaaarahh 

STDP 

A81 — STOP 

3. To run the Bayesian dot allocation 
There are two choices: 

a. Interactive — just type, after MCR>, "RUN A01.TSK$(ESC key)," 
and then answer the questions. 

b. Batch - submit a batch job with the following cards. 

$ JOB/NAME « SOMEONE /MCR/LIMIT = 360/ACCOUNT = 131 1 
$ CREATE POR007.DAT 

Data cards 

$EOD 

$MCR RUN A81.TSK 
$EOJ 


All program sources (Fortran) , object files, and task files can 
be accessed from UIC = [131,1] on the Image Processor (PDP 11/45) 
in JSC Building 17. 



APPENDIX A 


SUBPROGRAMS OF A81: TO STRIP A GROUND TRUTH MAP 
A.l MR6 

The subprogram MR6 determines a majority rule label for a group 
of six labels. 

A. 1.1 LINKAGES 

This routine does not call any other subprogram. It is called by 
the MAIN program, A81. 

A. 1.2 INTERFACES 

The MR6 subprogram interfaces with other routines through the 
calling arguments. 

A. 1.3 INPUTS 

Calling sequence: CALL MR6 (IG, LABEL) 


Parameter 

Type 

Dimension 

In/Out 

Definition 

IG 

Byte 

6 

In 

Array containing the six 
labels 

LABEL 

Byte 

1 

Out 

The chosen majority rule 
label 


A. 1.4 OUTPUTS 

The results are returned for use by the calling routine. 

A. 1.5 STORAGE REQUIREMENTS 

This subprogram requires little storage. 



A. 1.6 DESCRIPTION 


The subprogram MR6 scans the Input array IG In ascending order , 
counts the number of identical labels , and selects the label with 
the greatest count. If two labels have that count, the one 
appearing first will be selected. 

A. 1.7 LISTING 

The subprogram listing is provided in section 3.7. 



APPENDIX B 


SUBPROGRAMS OF A82: PROPORTIONAL DOT ALLOCATION 


B.l SAE 

The subprogram SAE computes a proportion estimate using propor- 
tional dot allocation. 

B.1.1 LINKAGES 

This routine calls the CLMPCS, GETDOT , GTMPLB , and LBLITP sub- 
programs. It is called by the MAIN program, A82. 

B.l. 2 INTERFACES 

The SAE subprogram interfaces with other routines through common 
block PRTFLG . 

B.l. 3 INPUTS 

Calling sequence: CALL SAE(NDOT,PSG) 


Parameter 

IZE®. 

Dimension 

In/Out 

Definition 

NDOT 

1*2 

1 

In 

Total number of dots to 
be allocated 

PSG 

R*4 

1 

Out 

The proportion estimate 
for small grain 


B.l. 4 OUTPUTS 

The results are returned for use by the calling routine. 
B . 1 . 5 STORAGE REQUIREMENTS 

This subprogram requires a fair amount of storage. 

B.l. 6 DESCRIPTION 

For a description of subprogram SAE, see section 5.'.. 


/^3 


B.1.7 FLOW CHART 

The subprogram flow chart Is provided in section 5.6. 
B.1.8 LISTING 

The subprogram listing is provided in section 5.7. 


APPENDIX C 


SUBPROGRAMS OF A83: PROPORTIONAL DOT ALLOCATION, 
MAJORITY RULE LABELING 


C.l SAEMR 

The subprogram SAEMR computes a proportion estimate using propor- 
tional dot allocation and majority rule labeling. 

C.1.1 LINKAGES 

This routine calls the GETDOT , CLMPCS, GTMPLB, LBLITP, and MR 
subprograms. It is called by the MAIN program, A83. 

C.l. 2 INTERFACES 

The SAEMR subprogram interfaces with other routines through common 
block PRTFLG . 

C.l. 3 INPUTS 

Calling sequence: CALL SAEMR (NDOT,PSG) . 


Parameter 

Type 

Dimension 

In/Out 

Definition 

NDOT 

1*2 

1 

In 

Total number of dots to 
be allocated 

PSG 

R*4 

1 

Out 

The proportion estimate 
for small grain 


C.l. 4 OUTPUTS 

The results are returned for use by the calling routine. 
C.l. 5 STORAGE REQUIREMENTS 

This subprogram requires a fair amount of storage. 


/KS" 



C.1.6 DESCRIPTION 

For a description of subprogram SAEMR, see section 6.1. 
C.1.7 FLOW CHART 

The subprogram flow chart is provided in section 6.6. 
C.1.8 LISTING 

The subprogram listing is provided in section 6.7. 




APPENDIX D 


SUBPROGRAMS OF A84; BAYESIAN DOT ALLOCATION 
(UNIFORM PRIOR) 


D.l SAEBl 

The subprogram SAEBl computes a proportion estimate using Bayesian 
dot allocation (uniform prior) . This subprogram has two entry 
points t 

SAEBl - to compute an estimate 

INITl - to set the number of initial dots 

D.1.1 LINKAGES 

This routine calls the CLMPCS, GETDOT, CLMPXY, GTMPLB , and LBLITP 
subprograms. It is called by the MAIN program, A84. 

D.l. 2 INTERFACES 

The SAEBl subprogram interfaces with other routines through common 
block PRTFLG. 

D.l. 3 INPUTS 

Calling sequence; CALL SAEBl (MAXDOT,PSG) 


Parameter 

T^ge 

Dimension 

In/Out 

Definition 

MAXDOT 

1*2 

1 

In 

Maximum number of dots to 
be assigned 

PSG 

R*2 

Variable 

Out 

The proportion estimate 
for small grain 

Entry calling sequence: CALL 

INITl (NNIT) 

Parameter 

Type 

Dimension 

In/Out 

Definition 

NNIT 

1*2 

1 

In 

Number of initial dots to 
be assigned to each cluster 




D.1.4 OUTPUTS 


The results are returned for use by the calling routine. 
D . 1 . 5 STORAGE REQUIREMENTS 

This subprogram requires a fair amount of storage. 


0.1.6 DESCRIPTION 

For a description of subprogram SAEBl , see section 7.1. 
D.1.7 FLOW CHART 

The subprogram flow chart is provided in section 7.6. 


0.1.8 LISTING 

The subprogram listing is provided in section 7.7. 




APPENDIX E 


SUBPROGRAMS OF A85: BAYESIAN DOT ALLOCATION 

(NO PRIOR) 


E.l SAEB2 

The subprogram SAEB2 computes a proportion estimate using Bayesian 
dot allocation (no prior) . This subroutine has two entry points: 
SAEB2 — to compute an estimate 
INIT2 — to set the number of initial dots 

E.1.1 LINKAGES 

This routine calls the CLMPCS, GETDOT, CLMPXY, GTMPLB , and LBLITP 
subprograms. It is called by the MAIN program, A85. 

E.l. 2 INTERFACES 

The 3AEB2 subprogram interfaces with other routines through common 
block PRTFLG . 

E.l. 3 INPUTS 

Calling sequence: CALL SAEB2 (MAXDOT , PSG) 


Parameter 

T yp e 

Dimension 

In/Out 

Definition 

MAXDOT 

1*2 

1 

In 

Maximum number of dots to 
be assigned 

PSG 

R*2 

Variable 

Out 

The proportion estimate 


for small grain 

Entry calling sequence: INIT2(NNIT) 


Parameter 

Type 

Dimension 

In/Out 

Definition 

NNIT 

1*2 

1 

In 

Number of initial dots to 


be assigned to each cluster 


/W 


E.1.4 OUTPUTS 


The results are returned for use by the calling routine. 
E.1.5 STORAGE REQUIREMENTS 

This subprogram requires a fair amount of storage. 


E.1.6 DESCRIPTION 

For a description of subprogram SAEB2, see section 8.1. 
E.1.7 FLOW CHART 

The subprogram flow chart is provided in section 8.6. 
E.1.8 LISTING 

The subprogram listing is provided in section 8.7. 




APPENDIX F 


SUBPROGRAMS OF A87: BAYESIAN DOT ALLOCATION 
(QUADRATIC PRIOR) 

F.l SAEB3 

The subprogram SAEB3 computes a proportion estimate using Bayesian 
dot allocation (quadratic prior) . This subprogram has two entry 
points: 

SAEB3 — to compute an estimate 

INIT3 — to set the number of initial dots 

F.1.1 LINKAGES 

This routine calls the CLMPCS, GETDOT, CLMPXY, GTMPLB , THl , and 
LBLITP subprograms. It is called by the MAIN program, A87. 

F.l. 2 INTERFACES 

The SAEB3 subprogram interfaces with other routines through common 
block PRTFLG . 

F . 1 . 3 INPUTS 

Calling sequence: CALL SAEB3 (MAXDOT,PSG) 


Parameter 

Type 

Dimension 

In/Out 

Definition 

MAXDOT 

1*2 

1 

In 

Maximum number of dots to 
be assigned 

PSG 

R*2 

Variable 

Out 

The proportion estimate 
for small grain 


Entry calling sequence: CALL INIT3 (NNIT) 

Parameter Type Dimension In/Out Definition 


NNIT 


1*2 


1 


In 


Number of initial dots to 
be assigned to each cluster 


P.1.4 OUTPUTS 

The results are returned for use by the calling routine. 

P.l.D STORAGE REQUIREMENTS 

This subprogram requires a fair amount of storage. 

P.1.6 DESCRIPTION 

For a description of subprogram SAEB3, see section 9.1. 

F . 1 . 7 FLOW CHART 

The subprogram flow chart is provided in section 9.6. 

F.1.8 LISTING 

The subprogram listing is provided in section 9.7. 

F.2 TH1 

A 

The subprogram THl computes an estimate 0 (n,x) for use in connec- 
tion with Bayesian dot allocation (quadratic prior) . 

F . 2 . 1 LINKAGES 

This routine does not call any other subprogram. It is called by 
subprogram SAEB3 . 

F.2. 2 INTERFACES 

The THl subprogram interfaces with other routines through the 
calling arguments. 

F.2. 3 INPUTS 

Function usage: THl (N, NX) 


Parameter 

Type 

Dimension 

Ir./Out 


Definition 

N 

1*2 

1 

In 

Number 

of dots assigned 

NX 

1*2 

1 

In 

Number 

of small-grain dots 


assigned 


F.2.4 OUTPUTS 

The function's value is returned. 

F.2.5 STORAGE REQUIREMENTS 

This subprogram requires small storage. 

F.2.6 DESCRIPTION 

For a description of subprogram THl, see section 9.1. 
F .2.1 FLOW CHART 

The subprogram flow chart is provided in section 9.6. 
F.2.8. LISTING 

The subprogram listing is provided in section 9.7. 


^J>3 — 

/S3 


APPENDIX G 


SUBPROGRAMS OF A89: BAYESIAN DOT ALLOCATION 
(MODIFIED QUADRATIC PRIOR) 


G.l SAEB4 

The subprogram SAEB4 computes a proportion estimate using Bayesian 
dot allocation (modified quadratic prior) . This subprogram has 
two entry points: 

SAEB4 - to compute an estimate 

INIT4 - to set the number of initial dots 


G.1.1 LINKAGES 

This routine calls the CLMPCS, GETDOT, CLMPXY, GTMPLB , TH2 , and 
LBLITP subprograms. It is called by the MAIN program, A89. 

G.l. 2 INTERFACES 

The SAEB4 subprogram interfaces with other routines through common 
block PRTFLG . 


G.l. 3 INPUTS 

Calling sequence: CALL SAEB4 (MAXDOT,PSG) 


Parameter 

Type 

Dimension 

In/Out 

Definition 


MAXDOT 

1*2 

1 

In 

Maximum number of 
be assigned 

dots to 

PSG 

R*2 

Variable 

Out 

The proportion estimate 
for small grain 

Entry calling sequence: CALL 

INIT4 (NNIT) 


Parameter 

1X21 

Dimension 

In/Out 

Definition 


NNIT 

1*2 

1 

In 

Number of initial 

dots to 


be assigned to each cluster 


/sv 


6.1.4 OUTPUTS 


The results are returned for use by the calling routine. 

G.1.5 STORAGE REQUIREMENTS 

This subprogram requires a fair amount of storage. 

G.1.6 DESCRIPTION 

For a description of subprogram SAEB4, see section 10.1. 

G.1.7 FLOW CHART 

The subprogram flow chart is provided in section 10.6. 

G.1.8 LISTING 

The subprogram listing is provided in section 10.7. 

G.2 TH2 

The subprogram TH2 computes an estimate 0 (n,x) for use in connec 
tion with Bayesian dot allocation (modified quadratic prior) . It 
has two entries: 

TH2 - return the value of the function 
TH2PAR — to set the parameters a, b, c 

G.2.1 LINKAGES 

This routine does not call any other subprogram. It is called by 
subprogram SAEB4. 

G.2. 2 INTERFACES 

The TH2 subprogram interfaces with other routines through the 
calling arguments. 



G.2.3 INPUTS 


Function usage: 

TH2 (N,NX) 



Parameter 

Type 

Dimension 

In/Out 

Definition 

N 

1*2 

1 

In 

Number of dots assigned 

NX 

1*2 

1 

In 

Number of small-grain do 
assigned 

Entry calling sequence: CALL 

TH2PAR(AX 

,BX,CX) 

Parameter 

T YE e 

Dimension 

In/Out 

Definition 

AX 

R*4 

1 

In 

Set parameter a to value 
in AX 

BX 

R*4 

1 

In 

Set parameter b to value 
in BX 

CX 

R*4 

1 

In 

Set parameter c to value 
in CX 

G.2.4 OUTPUTS 


i 



The function's value is returned. 

G . 2 . 5 STORAGE REQUIREMENTS 

This subprogram requires small storage. 

G.2.6 DESCRIPTION 

For a description of subprogram TH2, see section 10.1. The entry 
point TH2PAR provides a means of modifying the parameters a, b, 
c, which are initially set to a * 6, b ■ -7.877, and c » 2.9345. 

G.2.7 FLOW CHART 

The subprogram flow chart is provided in section 10.6. 

G.2.8 LISTING 

The subprogram listing is provided in section 10.7. 




APPENDIX H 


SUBPROGRAMS OF A91: BAYESIAN DOT ALLOCATION 
(ADAPTIVE PRIOR) 


H.l SAEB5 

The subprogram SAEB5 computes a proportion estimate using Bayesian 
dot allocation (adaptive prior) . This routine has two entry 
points: 

SAEB5 - to compute an estimate 

INIT5 - to set the number of initial dots 

H.1.1 LINKAGES 

This routine calls the CLMPCS, GETDOT , CLMPXY, GTMPLB , TH3, and 
LBLITP subprograms. It is called by the MAIN program, A91. 

H.l. 2 INTERFACES 

The SAEB5 subprogram interfaces with other routines through common 
block PRTFLG . 

H.l. 3 INPUTS 


Calling sequence: 

CALL SAEB5 (MAXDOT,: 

PSG) 


Parameter 

Type 

Dimension 

In/Out 

Definition 


MAXDOT 

1*2 

1 

In 

Maximum number of 

dots to 





be assigned 


PSG 

R*2 

Variable 

Out 

The proportion estimate 





for small grain 


Entry calling sequence: CALL 

INIT5 (NNIT) 


Parameter 

Type 

Dimension 

In/out 

Definition 


NNIT 

1*2 

1 

In 

Number of initial 

dots to 


be assigned to each cluster 


3 *^ 


ipffwr 


H.1.4 OUTPUTS 

Th« results are returned for use by the calling routine. 

H.1.5 STORAGE REQUIREMENTS 

This subprogram requires a fair amount of storage. 

H.1.6 DESCRIPTION 

For a description of subprogram SAEB5, see section 11.1. 

H.1.7. FLOW CHART 

The subprogram flow chart is provided in section 11.6. 

H.1.8 LISTING 

The subprogram listing is provided in section 11.7. 

H.2 TH3 

The subprogram TH3 computes an estimate 6 (n,x) for use in connec- 
tion with Bayesian dot allocation (adaptive prior) . It has four 
entries: 

TH3 - to return the value of the function 

TH3PAR - to set the parameters a, b, c, for the quadratic 
prior 

SETEP - to set the parameter a for the exponential prior 
SETPRI - to select the prior 

H.2.1 LINKAGES 

This routine does not call any other subprogram. It is called by 
subprogram SAEB5. 

H.2. 2 INTERFACES 

The TH3 subprogram interfaces with other routines through the 
calling arguments. 



H.2.3 INPUTS 

Function usage: TH3 (N,NX) 


Parameter 

Tyge 

Dimension 

In/Out 

Definition 

N 

1*2 

1 

In 

Number of dots assigned 

NX 

1*2 

1 

In 

Number of small-grain dots 
assigned 

Entry calling sequence: CALL 

TH3PAR(AX,BX,CX) 

Parameter 

Type 

Dimension 

In/Out 

Definition 

AX 

R*4 

1 

In 

Set parameter a to value 
in AX 

BX 

R*4 

1 

In 

Set parameter b to value 
in BX 

cx 

R*4 

1 

In 

Set parameter c to value 
in CX 

Entry calling sequence: CALL 

SETEP (ALPHA) 

Parameter 

Type 

Dimension 

In/Out 

Definition 

ALPHA 

R*4 

1 

In 

Set parameter a to value 
in ALPHA 

Entry calling sequence: CALL 

SETPRI (JFPR) 

Parameter 

T YP e 

Dimension 

In/Out 

Definition 

JFPR 

1*2 

1 

In 

Prior select: 


1 - quadratic prior 

2 - exponential prior 


H.2.4 OUTPUTS 

The function's value is returned. 

H.2.5 STORAGE REQUIREMENTS 

This subprogram requires small storage. 


H.2.6 DESCRIPTION 


For a description of subprogram TH3, see section 11.1. The abil- 
ity of the processor to select one prior (quadratic) or the other 
(exponential) at any instant cones from this subprogram, TH3. By 
calling SETPRI before using the function TH3, one can switch from 
one prior to the other. Notice that for quadratic prior, the 
parameters are initially set as follows t a ■ 6, b ■ -7.877, 
c ■ 2.9345. For the exponential prior, the parameter a is not 
initially set to any value. The calling program must set a by 
calling SETEP (ALPHA) before using the function. 

H.2.7 FLOW CHART 

The subprogram flow chart is provided in section 11.6. 

H.2.8 LISTING 

The subprogram listing is provided in section 11.7. 



APPENDIX I 


SUBPROGRAMS OF A88: BAYESIAN DOT ALLOCATION, 
MAJORITY RULE LABELING 


1.1 SAEBMR 

The subprogram SAEBMR computes a proportion estimate using 
Bayesian dot allocation and majority rule labeling. 

1.1.1 LINKAGES 

This routine calls the CLMPCS and BMR1 subprograms. It is called 
by the MAIN program, A88. 

1.1.2 INTERFACES 

The subprogram interfaces with other routines through common 
block PRTFLG . 

1.1.3 INPUTS 

Calling sequence: CALL SAEBMR (NDOT,PSG) 


Parameter 

Type 

Dimension 

In/Out 

Definition 

NDOT 

1*2 

1 

Out 

Total number of dots that 
are allocated 

PSG 

R*4 

1 

Out 

The proportion estimate 
for small grain 


1.1.4 OUTPUTS 

The results are returned for use by the calling routine. 

1.1.5 STORAGE REQUIREMENTS 

This subprogram requires a fair amount of storage. 


f(*f 


1.1.6 DESCRIPTION 


For a description of subprogram SAEBMR, see section 12.1. This 
subroutine uses the subroutine BMR1 (which in turn uses subroutine 
BMRSUP) to achieve the computation of the estimate. 

Subroutine SAEBMR is the executive routine to drive the rest of 
the routines. 

1.1.7 FLOW CHART 

The subprogram flow chart is provided in section 12.6. 

1.1.8 LISTING 

The subprogram listing is provided in section 12.7. 

1.2 BMRl 

The subprogram BMRl returns a small-grain or non- small -grain 
label for a cluster. 

1.2.1 LINKAGES 

This routine calls the BMRSUP subprogram. It is called by sub- 
program SAEBMR. 

1.2.2 INTERFACES 

The subprogram interfaces with other routines through common 
block PRTFLG . 

1.2.3 INPUTS 

Calling sequence: CALL BMRl (ISG,NS,NDARY,Kl f MM,NN,LL,M) 



Parameter 

Type 

Dimension 

In/Out 

Definition 

ISG 

1*2 

1 

Out 

Small-grain indicator: 
1 - yes, 0 - no 

NS 

1*2 

2 

Out 

Split of small-grain and 
non-small-grain dots 

NDARY 

1*2 

Variable 

In/Out 

Dot array 

Kl 

1*2 

1 

In 

The kth cluster the dots 
will be allocated to 

MM 

1*2 

Variable 

In 

Array containing sizes of 
all clusters 

NN 

1*2 

Variable 

In/Out 

Array containing numbers 
of dots assigned to all 
clusters 

LL 

1*2 

Variable 

In 

Array containing the clus- 
ters' numbers (the first 
cluster is 1, second 
cluster is 2, etc.) 

M 

1*2 

1 

In 

Total number of clusters 

1.2.4 OUTPUTS 




The results 

are returned for 

use by the 

calling routine. 


1.2.5 STORAGE REQUIREMENTS 

This subprogram requires small storage. 

1.2.6 DESCRIPTION 

For a description of subprogram BMRl, see section 12.1. This 
routine serves as a middleman between the executive subroutine 
SAEBMR and the workhorse routine BMRSUP. If the user wants to 
change the prior (from uniform prior) , only this routine, BMRl, 
needs reprogramming. 




1.2.7 FLOW CHART 


The subprogram flow chart is provided in section 12.6. 

1.2.8 LISTING 

The subprogram listing is provided in section 12.7. 

1.3 BMRSUP 

The subprogram BMRSUP performs the essential function of dot 
assignment and determination of the small-grain and non-small- 
grain split. 

1.3.1 LINKAGES 

This routine calls the GETDOT, GTMPLP , and LBLITP subprograms. 
It is called by subprogram BMRl. 

1.3.2 INTERFACES 

The subprogram interfaces with other routines through common 
block PRTFLG . 

1.3.3 INPUTS 


Calling 

sequence: 

CALL BMRSUP (JDOT, NS 

,MR , NDARY , KOLD ,MM , NN , LL , M) 

Parameter Type 

Dimension 

In/Out 

Definition 

JDOT 

1*2 

1 

In 

Number of dots requested 
to be assigned 

NS 

1*2 

2 

Out 

Split of small-grain and 
non-small-grain dots 

MR 

1*2 

1 

In/Out 

Majority indicator: 

1 — small grain 

2 — other 

NDARY 

1*2 

Variable 

In/Out 

Dot array 

HOLD 

1*2 

1 

In 

The kth cluster the dots 


will be allocated to 






Parameter 

Type 

Dimension 

In/Out 

Definition 

MM 

1*2 

Variable 

In 

Array containing sizes of 
all clusters 

NN 

1*2 

Variable 

In/Out 

Array containing numbers < 
dots assigned to all clus 
ters 

LL 

1*2 

Variable 

In 

Array containing the clus 
ters' numbers (the first 
cluster is 1, second 2, 
etc. ) 

M 

1*2 

1 

In 

Total number of clusters 


1.3.4 OUTPUTS 

The results are returned for use by the calling routine. 

1.3.5 STORAGE REQUIREMENTS 

This subprogram requires a fair amount of storage. 


1.3.6 DESCRIPTION 

For a description of subprogram BMRSUP, see section 12.1. The 
subroutine BMRSUP is the workhorse of the entire program A88; it 
assigns dots and determines labels. 


1.3.7 FLOW CHART 

The subprogram flow chart is provided in section 12.6. 


1.3.8 LISTING 

The subprogram listing is provided in section 12.7. 
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APPENDIX J 


SUBPROGRAMS OF A90t NON-BAYES IAN SEQUENTIAL DOT ALLOCATION, 

MAJORITY RULE LABELING 


J.l SAENBS 

The subprogram SAENBS computes a proportion estimate using non- 
Bayesian sequential dot allocation and majority rule labeling. 

J.1.1 LINKAGES 

This routine calls the CLMPCS and NBSMR subprograms. It is 
called by the MAIN program, A90. 

J.l. 2 INTERFACES 

The subprogram interfaces with other routines through common 
block PRTFLG. 

J.l. 3 INPUTS 

Calling sequence: CALL SAENBS (NDOT,PSG) 


Parameter 

T y£f. 

Dimension 

In/Out 

Definition 

NDOT 

1*2 

1 

Out 

Total number of dots that 
are allocated 

PSG 

R*4 

1 

Out 

The proportion estimate 
for small grain 


J . 1 . 4 OUTPUTS 

The results are returned for use by the calling routine. 
J.l. 5 STORAGE REQUIREMENTS 

This subprogram requires a fair amount of storage. 


J.1.6 DESCRIPTION 


For a description of subprogram SAENBS, see section 13.1. This 
subroutine uses NBSMR to achieve the computation of the estimate. 

J.1.7 FLOW CHART 

The subprogram flow chart is provided in section 13.6. 

J.1.8 LISTING 

The subprogram listing is provided in section 13.7. 

J.2 NBSMR 

The subprogram NBSMR assigns dots to a cluster according to non- 
Bayesian sequential majority rule. 

J.2.1 LINKAGES 

This routine calls the GETDOT, CLMPXY, GTMPLB, and LBLITP sub- 
programs. It is called by subprogram SAENBS. 

J . 2 . 2 INTERFACES 

The subprogram interfaces with other routines through common 
block PRTFLG . 

J.2. 3 INPUTS 

Calling sequence: CALL NBSMR ( ISG , NS , NDARY , K1 ,MM , NN , LL , NX , M) 


Parameter 

Type 

Dimension 

In/Out 

Definition 

ISG 

1*2 

1 

Out 

Small-grain indicator: 
1 — yes, 0 — no 

NS 

1*2 

2 

Out 

Split of small-grain and 
non-small-grain dots 

NDARY 

1*2 

Variable 

In/Out 

Dot array 


j** 


Parameter 

Type 

Dimension 

In/Out 

Definition 

K1 

1*2 

1 

In 

The kth cluster the dots 
will be assigned to 

MM 

1*2 

Variable 

In 

Array containing sizes of 
all clusters 

NN 

1*2 

Variable 

In/Out 

Array containing numbers 
of dots assigned to all 
clusters 

LL 

1*2 

Variable 

In 

Array containing the clus 
ters' numbers (the first 
cluster is 1, second 2, 
etc. ) 

NX 

1*2 

1 

In/Out 

Array containing numbers 
of small-grain dots 
assigned to the clusters 

M 

1*2 

1 

In 

Total number of clusters 

Entry calling sequence; CALL 

INIT(NNIT) 


Parameter 

Type 

Dimension 

In/Out 

Definition 

NNIT 

1*2 

1 

In 

Number of initial dots 
desired 


J . 2 . 4 OUTPUTS 

The results are returned for use by the calling routine. 
J . 2 . 5 STORAGE REQUIREMENTS 

This subprogram requires a fair amount of storage. 


J.2.6 DESCRIPTION 

For a description of subprogram NBSMR , see section 13.1. This 
subroutine is the workhorse used by subroutine SAENBS to compute 
the estimate. 


J.2.7 FLOW CHART 


The subprogram flow chart is provided in section 13.6. 
J.2.8 LISTING 

The subprogram listing is provided in section 13.7. 



